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ABSTRACT 


This  thesis  is  concerned  with  the  development  of  a  model  to  analyze  a  Digital  Im¬ 
age  Synthesizer  (DIS)  integrated  circuit  designed  to  create  false  target  images  to  deceive 
Inverse  Synthetic  Aperture  Radar  (ISAR).  The  DIS  is  able  to  recreate  the  scattering  ef¬ 
fect  of  a  moving  target  by  using  appropriate  phase  and  gain  modulations  on  an  inter¬ 
cepted  ISAR  chirp  signal  before  retransmitting  it  with  the  proper  time  delay.  The  DIS 
signal  processing  and  the  ISAR  compression  of  the  modulated  return  are  modeled  to  ex¬ 
amine  the  range-DoppIer  profile  of  a  synthesized  false  target  image.  The  image  is  repre¬ 
sentative  of  the  image  that  would  appear  on  an  ISAR  display.  ISAR  image  quality  is 
used  to  evaluate  different  DIS  architectures  and  bit  formats.  Evaluation  of  the  image 
quality  is  based  on  the  deviation  from  an  infinite  resolution  false  target  image.  The  re¬ 
sults  obtained  from  evaluating  different  DIS  architectures  indicate  that  the  design  is  toler¬ 
ant  of  significant  quantization  errors.  The  model  is  used  to  validate  the  architecture  of 
the  integrated  circuit  being  fabricated.  Finally,  various  different  ISAR  integration  times 
and  pulse  repetition  frequencies  are  used  to  confirm  the  integrity  of  the  model. 
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EXECUTIVE  SUMMARY 


This  thesis  is  concerned  with  the  development  of  a  model  to  analyze  the  perform¬ 
ance  of  a  Digital  Image  Synthesizer  (DIS)  integrated  circuit  designed  to  create  false  tar¬ 
get  images  to  deceive  an  Inverse  Synthetic  Aperture  Radar  (ISAR).  The  DIS  chip  is 
being  designed  at  the  Naval  Postgraduate  School,  Center  of  Joint  Services  Electronic 
Warfare.  The  circuit  is  able  to  recreate  the  scattering  effect  of  a  moving  target  by  using 
appropriate  phase  and  gain  modulations  on  an  intercepted  ISAR  chirp  signal  before  re¬ 
transmitting  it  with  the  proper  time  delays.  The  DIS  is  expected  to  be  integrated  with  the 
Advance  Multifunction  RF  Concept  (AMRFC)  test  bed.  The  AMRFC  is  radio  frequency 
system  that  will  integrate  radar,  communications  and  electronic  warfare  equipment,  and 
significantly  increase  the  effectiveness  of  the  shipboard  electronic  warfare  (EW)  system 
while  reducing  the  number  of  topside  antennas. 

The  model  that  is  developed  in  this  thesis  is  able  to  simulate  the  DIS  operation  by 
generating  the  required  input  signals  and  producing  the  chip  signal  processing  to  obtain 
the  output-jamming  signal.  A  second  part  of  the  model  uses  the  DIS  output  to  simulate 
the  compression  process  that  takes  place  in  the  ISAR  receiver.  The  output  of  the  ISAR 
range  and  azimuth  compression  is  an  image  with  the  range-Doppler  profile  of  the  false 
target.  This  image  is  representative  of  the  image  that  would  appear  on  an  ISAR  display. 

ISAR  image  quality  is  used  to  evaluate  different  DIS  architectures  and  bit  for¬ 
mats.  Evaluation  of  the  image  quality  is  based  on  the  deviation  from  an  infinite  resolu¬ 
tion  false  target  image.  The  results  obtained  from  evaluating  different  DIS  architectures 
indicate  that  the  design  is  tolerant  of  significant  quantization  errors.  The  model  is  used  to 
validate  the  architecture  of  the  integrated  circuit  being  fabricated. 

Additional  studies  are  performed  to  find  the  relation  between  the  dynamic  range 
of  the  DIS  output  and  the  image  contrast  at  the  ISAR  display.  Finally,  various  different 
ISAR  integration  times  and  pulse  repetition  frequencies  are  used  in  the  model  in  order  to 
study  the  image  effects  of  different  Doppler  resolutions.  These  studies  provide  helpful 
results  for  the  design  effort,  since  they  identify  critical  parameters  and  confirm  the  integ¬ 
rity  of  the  model. 
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I.  INTRODUCTION 

A.  DIGITAL  IMAGE  SYNTHESIZER  OPERATIONAL  CONCEPT 

The  Digital  Image  Synthesizer  (DIS)  is  an  application  specific  integrated  circuit 
(ASIC)  able  to  generate  false  target  images  to  deceive  an  Inverse  Synthetic  Aperture  Ra¬ 
dar  (ISAR).  ISAR  is  a  high-resolution  radar  technique  capable  of  developing  a  two- 
dimensional  image  of  a  target  based  on  the  relative  Doppler  of  the  reflecting  components. 
The  DIS  is  able  to  generate  a  set  of  scattering  effects  that  emulate  the  presence  of  multi¬ 
ple  false  targets,  thus  attempting  to  deceive  the  ISAR  target  recognition  process. 

The  DIS  design  takes  advantage  of  readily  available  digital  radio  frequency  mem¬ 
ory  (DRFM)  technology  to  sample,  store  and  replicate  an  intercepted  ISAR  signal.  The 
DRFM  is  used  to  feed  the  DIS  circuit  with  the  intercepted  ISAR  signal  in  the  form  of 
digitized  samples  of  the  chirp  pulses.  The  DIS  circuit  modulates  the  signal  samples  in 
amplitude  and  phase  in  order  to  recreate  the  scattering  effect  of  a  desired  false  target. 
The  DIS  output  (modulated)  signal  is  sent  through  a  digital-to-analog  converter  to  a  radio 
frequency  (RF)  module  for  frequency  up-conversion  and  retransmission.  Note  that  a 
DRFM  alone  cannot  perform  the  complex  modulations  required  to  synthesize  large  false 
targets. 

The  ISAR  platform  collects  the  returned  modulated  pulses  and  coherently  proc¬ 
esses  them  to  form  an  image.  The  RF  receiver  first  down-converts  the  signal  frequency, 
samples  the  signal,  and  performs  the  pulse  compression  digitally.  The  signal  is  then 
sampled  in  both  in-phase  (I)  and  quadrature  phase  (Q),  and  sent  to  the  ISAR  signal  proc¬ 
essor.  After  range  compression,  the  samples  within  each  of  the  pulses  are  integrated  co¬ 
herently  and  resolved  in  Doppler  (azimuth  compression  process)  to  display  a  two- 
dimensional  target  image.  A  block  diagram  of  this  process  is  presented  in  Figure  1. 

Technical  developments  in  naval  warfare  have  improved  the  capabilities  of  anti¬ 
ship  cruise  missile  (ASCM)  strike  systems,  which  are  still  considered  the  foremost  threat 
to  naval  vessels.  The  range  of  these  missiles  has  been  significantly  increased  and  capa¬ 
bilities  such  as  supersonic  speeds,  high  maneuverability,  data-link  control,  GPS  aided 
guidance,  re-attacking  options  and  active/passive  seekers  have  been  incorporated. 
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Launching  platforms  with  their  increased  missile  ranges  have  forced  an  increase  in  their 
sensor  capabilities  in  order  to  properly  evaluate  potential  hostile  targets  at  greater  dis¬ 
tances.  This  problem  has  also  resulted  in  maritime  patrol  and  strike  aircrafts,  such  as  the 
U.S.  Navy  P-3C  and  S-3B,  being  fitted  with  ISAR  in  order  to  identify  surface  targets 
based  on  their  profile  before  assigning  a  strike  mission. 


ISAR  RADAR 


PROTECTED  TARGET 


Figure  1  Block  Diagram  of  the  DIS  Operational  Concept. 

The  ASCM  strike  can  be  divided  into  six  phases:  surveillance,  identification,  targeting, 
missile  launch,  guidance  and  the  terminal  homing  phase.  Naval  vessels  attempt  to 
counter  this  threat  at  all  phases.  Numerous  tactics  have  been  developed  to  counter  each 
phase,  but  are  not  within  the  scope  of  this  thesis.  However,  it  must  be  mentioned  that  the 
DIS  concept  can  be  used  for  a  number  of  situations  including  counter-identification, 
counter-targeting  and  counter-terminal  (ISAR  seekers). 

A  common  counter-identification  scenario  to  illustrate  the  DIS  operational  con¬ 
cept  would  be  a  destroyer  being  detected  by  a  maritime  patrol  aircraft  (MPA)  at  a  range 
beyond  100  nautical  miles.  The  aircraft  knows  the  ship’s  position  and  has  tracked  its 
course  and  speed.  The  aircraft’s  next  step  is  to  attempt  to  evaluate  the  unknown  target. 
Electronic  support  (ES)  could  provide  radar  interceptions  suitable  for  identification,  but 
the  ship,  once  aware  of  the  presence  of  an  MPA,  will  have  restricted  its  emissions.  The 
range  is  too  great  for  infrared  (IR)  imaging  and  is  aggravated  by  drizzle  that  severely  af¬ 
fects  IR  propagation.  The  MPA  decides  to  switch  the  radar  to  ISAR  mode  in  order  to 
identify  the  silhouette  of  the  ship’s  hull.  To  obtain  the  best  image,  the  MPA  positions 
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itself  45  degrees  off  the  ship’s  bow.  The  ship  intercepts  the  ISAR  chirp  signal  and  acti¬ 
vates  the  DIS  countermeasure.  Based  on  the  aircraft’s  relative  location  obtained  from  the 
ES,  the  DIS  is  then  able  to  generate  a  ship  silhouette  similar  to  a  merchant  vessel  that  was 
supposed  to  have  transited  the  area  that  same  day.  In  addition,  the  ship  activates  a  com¬ 
mercial  navigation  radar  to  make  the  presence  of  a  merchant  ship  more  credible.  The 
ISAR  display  in  the  MPA  shows  the  silhouette  of  a  merchant  ship,  in  agreement  with 
what  was  expected  to  transit  that  area,  which  is  also  supported  by  the  emission  of  a  com¬ 
mercial  radar.  The  MPA  classifies  the  target  as  neutral  and  leaves  the  scene. 

A  second  scenario  for  employing  a  DIS  is  on  an  aircraft  countering  Non- 
Cooperative  Target  Recognition  (NCTR)  techniques  based  on  scattering  characteristics. 
Radar  NCTR  is  the  classification  or  identification  of  a  type  of  aircraft  using  radar  based 
techniques,  including  range  profile,  2D-scattering  centers,  time-frequency  analysis,  natu¬ 
ral  frequencies  and  feature  extraction.  A  properly  programmed  DIS  mounted  on  a  threat 
aircraft  would  be  able  to  replicate  the  scattering  characteristics  of  a  friendly  aircraft. 

Further,  DIS  applications  can  also  include  the  generation  of  false  land  targets  to 
protect  an  area  against  a  reconnaissance  aircraft  mapping  with  a  Synthetic  Aperture  Ra¬ 
dar  (SAR)  or  any  defensive  effort  against  platforms  using  high  resolution  Doppler  proc¬ 
essing  radars. 

The  DIS  is  currently  being  designed  at  the  Center  for  Joint  Services  Electronic 
Warfare  at  the  Naval  Postgraduate  School  and  is  expected  to  be  integrated  with  the  Ad¬ 
vance  Multifunction  RF  Concept  (AMRFC)  test  bed.  The  AMRFC  is  an  RF  system  that 
will  integrate  radar,  communications  and  electronic  warfare  equipment,  and  will  signifi¬ 
cantly  increase  the  effectiveness  of  the  shipboard  electronic  warfare  (EW)  system  while 
reducing  the  number  of  topside  antennas. 

B.  PRINCIPAL  CONTRIBUTIONS 

The  objective  of  this  thesis  is  to  study  the  DIS  design  and  develop  an  updated  ver¬ 
sion  of  the  model.  The  principal  contributions  in  this  thesis  can  be  identified  as  follows: 

•  Decrease  the  DIS  simulation  processing  time  by  increasing  the  efficiency 
of  the  mathematical  models  and  reducing  their  complexity 
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•  Incorporate  a  model  to  simulate,  measure  and  manage  the  overflow  effect 
within  the  DIS  output  adders 

•  Develop  a  model  to  simulate  a  DIS  with  infinite  bit  resolution 

•  Create  a  model  to  compare  and  correlate  the  ISAR  image  with  the  infinite 
resolution  image 

•  Change  the  DIS  and  ISAR  models  to  enable  the  evaluation  of  chip  proc¬ 
essing  without  the  preceding  constant  phase  sequence  that  is  currently 
used 

•  Evaluate  the  impact  of  the  DIS  architecture  and  bit  formats  on  the  output 
image  quality 

•  Evaluate  the  impact  of  the  different  ISAR  parameters  and  sampling  crite¬ 
ria  on  the  output  image  quality 

•  Incorporate  the  capability  to  model  a  512-range-bin  target 

•  Produce  input  and  output  vectors  to  support  the  very  high-speed  integrated 
circuit  hardware  description  language  (VHDL)  simulations.  (VHDL  is  the 
hardware  description  language  used  for  documentation,  verification  and 
synthesis  of  large  digital  designs  currently  used  in  the  DIS  project.) 

•  Create  a  common  software  version  usable  for  different  types  of  targets  and 
radar  signal  parameters 

In  summary,  this  thesis  provides  the  DIS  design  team  with  an  efficient  modeling  tool  to 
support  the  DIS  chip  design  efforts,  by  being  able  to  accurately  predict  the  performance 
of  various  DIS  architecture  configurations. 

C.  THESIS  OUTLINE 

This  thesis  first  introduces  the  concepts  of  the  ISAR  and  DIS.  A  description  of 
the  DIS  software  model  and  the  performance  studies  obtained  with  its  simulation  follows. 
Chapter  II  presents  the  SAR  concept  and  its  angular  resolution  enhancement.  Next, 
ISAR  is  introduced  with  a  description  of  its  resolution,  the  output  target  image  and  the 
basics  of  its  compression  signal  processing.  This  analysis  provides  a  background  for  un- 
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derstanding  the  math  models  used  to  simulate  ISAR  reception  and  image  generation. 
Chapter  III  mathematically  describes  the  DIS  process  from  the  generation  of  the  target 
modulation  parameters  to  the  complex  pulse  modulation.  It  also  includes  a  mathematical 
description  of  the  ISAR  range  and  azimuth  compression  processes  used  in  the  model  to 
generate  the  expected  target  image.  Chapter  IV  describes  the  DIS  model  developed  in 
MATLAB.  The  highlights  of  the  software  are  explained  using  selected  code  lines,  plots 
with  signal  characteristics  and  tables  with  values  obtained  from  simulation  runs.  Finally, 
this  chapter  describes  the  modifications  required  to  obtain  an  ideal  model  with  infinite  bit 
resolution.  Chapter  V  studies  the  performance  of  the  DIS  design  using  the  simulation 
model  presented  in  Chapter  IV.  The  image  obtained  from  the  current  DIS  design  is  com¬ 
pared  to  an  infinite  bit  resolution  image  to  quantify  quantization  and  resolution  errors. 
The  effects  from  individual  stages  of  the  DIS  process  are  isolated  to  numerically  evaluate 
their  contribution  to  the  image  quality.  Chapter  VI  examines  the  effectiveness  of  the  DIS 
against  an  ISAR  using  different  pulse  repetition  frequencies  (PRFs)  and  integration  times. 
The  relation  between  ISAR  signal  parameters  and  image  Doppler  resolution  is  investi¬ 
gated  in  theory  and  demonstrated  with  simulation  outputs.  Chapter  VII  studies  the  dy¬ 
namic  range  of  the  DIS  signal  and  its  effect  on  ISAR  target  image  contrast.  Theoretical 
dynamic  range  values  are  developed  and  compared  with  results  obtained  from  simula¬ 
tions.  Chapter  VIII  presents  conclusions  to  this  thesis  study  and  recommendations  for 
future  efforts  in  designing  and  testing  the  DIS. 

This  thesis  also  includes  three  appendices  to  support  the  information  presented. 
Appendix  A  contains  the  MATLAB  m-files  that  comprise  the  DIS  simulation  model. 
Appendix  B  contains  a  brief  description  of  the  Digital  RF  Memory  (DRFM),  which  is  an 
integral  part  of  the  DIS  concept.  Finally,  Appendix  C  contains  a  brief  description  of  the 
AN/APS-I37,  an  ISAR  radar  operated  by  the  U.S.  Navy  onboard  the  P3-C  and  the  S-3B 
aircrafts. 
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II.  SAR/ISAR  FUNDAMENTALS 


ISAR  is  a  high-resolution  radar  usually  operated  from  airborne  platforms.  The 
high  resolution  obtained  by  this  radar  in  both  range  and  azimuth  makes  it  possible  to  ob¬ 
tain  the  image  of  a  moving  target  in  two  dimensions.  Among  the  most  popular  opera¬ 
tional  applications  of  ISAR  is  that  of  a  maritime  patrol  aircraft  tracking  a  surface  ship.  In 
this  application,  ISAR  exploits  the  ship’s  movement  to  obtain  a  two-dimensional  image 
of  its  length  and  height.  This  broadside  aspect  is  very  suitable  for  target  recognition,  but 
as  will  be  demonstrated  in  later  chapters,  the  use  of  electronic  attack  (EA)  makes  it  vul¬ 
nerable  to  deception.  This  chapter  describes  the  basic  principles  of  SAR  and,  more  spe¬ 
cifically,  those  of  ISAR.  This  background  information  helps  the  reader  understand  the 
fundamentals  of  the  DIS  EA  in  the  chapters  to  follow. 

A.  RADAR  AZIMUTH  RESOLUTION 

For  a  radar  antenna,  the  angular  resolution  is  approximately  equivalent  to  the  an¬ 
gular  width  of  the  antenna  beam.  Using  the  definition  of  arc-length,  the  radar  azimuth 
resolution  can  be  approximated  as: 

‘>a„,.=  R0y,.  (2-1) 

where  R  is  the  range  to  the  target  and  is  the  half-power  beamwidth  in  radians.  A  uni¬ 
formly  illuminated  sin(x)/x  linear  aperture  antenna  has  a  half  power  beamwidth  given  by: 

=  0-88  ^ID  (2.2) 

where  T  is  the  wavelength  of  the  carrier  frequency  and  D  is  the  length  of  the  physical 
antenna  [1].  The  beamwidth  expression  in  (2.2)  can  be  simplified  using  a  4dB  width 
equal  to  XI D  .  The  resulting  azimuth  resolution  for  the  conventional  antenna  becomes: 

d.„.,=  R^/D.  (2.3) 

From  Equation  (2.3),  it  can  be  concluded  that  an  extremely  narrow  beam  is  re¬ 
quired  to  obtain  a  fine  azimuth  resolution  at  a  long  range.  Such  a  beam  can  be  obtained 
either  with  a  big  antenna  or  a  very  small  wavelength.  Big  antennas  are  impractical  when 
operated  from  an  aircraft  because  of  size,  weight  and  aerodynamic  considerations.  On 
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the  other  hand,  a  radar  signal  with  a  very  small  wavelength  is  severely  affeeted  by  the 
strong  atmospherie  attenuation  oeeurring  at  very  high  frequeneies. 

B.  SAR  CONCEPT 

SAR  takes  advantage  of  the  motion  of  an  aireraft  to  simulate  a  large  antenna  in 
order  to  obtain  a  fine  azimuth  resolution.  An  aireraft  flying  a  straight  path  with  a  eoher- 
ent  side  looking  pulsed  radar  is  used  to  ereate  an  extremely  large  synthetie  aperture.  The 
equivalent  aperture  size  of  this  antenna  is  the  flight  path  length  L. 

SAR  pulses  are  transmitted  at  equally  spaeed  locations  along  the  flight  path.  The 
complex  samples  for  each  pulse  of  the  returning  signal,  amplitude  and  phase,  are  stored 
in  memory  and  added  coherently.  This  summation  results  in  a  similar  effect  as  a  physical 
linear  array  antenna  of  aperture  length  L. 

The  main  difference  between  the  synthetic  aperture  and  a  physical  array  antenna, 
is  that  the  SAR  beam  and,  thus  its  resolution  is  finer  by  a  factor  of  two  [2].  This  distinc¬ 
tion  occurs  because  physical  array  antennas  resolve  angle  selectivity  using  the  phase  dif¬ 
ference  at  the  reception  process  (one  way  trip).  Conversely,  SAR  resolves  azimuth  using 
phase  differences  resulting  from  round  trip  distance  differences  between  aperture  ele¬ 
ments.  By  inserting  a  factor  of  two  and  using  the  aircraft  path  length  L  as  the  aperture 
size  in  Equation  (2.2),  the  SAR  effective  beamwidth  becomes: 

6,jr,„=  0.44  1 /L.  (2.4) 

From  Equation  (2.1)  and  the  SAR  effective  beamwidth  in  (2.4),  it  follows  that  the  equiva¬ 
lent  azimuth  resolution  can  be  expressed  as: 

0.44  Rl/L.  (2.5) 

To  simplify  the  expression  in  (2.5),  a  4-dB  SAR  effective  beamwidth  of  A/D  can  be 
used,  yielding  an  azimuth  resolution  of: 

-R -1/(22.).  (2.6) 

In  summary,  a  side  looking  SAR,  using  a  physically  wide  beam,  is  able  to  synthe¬ 
size  a  narrow  beam,  whose  width  and  azimuth  resolution  depends  on  the  range  R  and  the 
aircraft  path  length  L.  A  simple  analysis  may  state  that  maximizing  the  distance  of  path 
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length  by  extending  the  integration  time  would  minimize  the  resolution.  However,  there 
is  a  point  at  which  a  further  increase  in  the  synthetic  aperture  length  degrades  resolution. 

1.  Unfocused  SAR 

Assuming  that  the  SAR  aircraft  path  follows  a  straight  line,  the  range  from  each 
point  in  the  path  to  a  distant  target  differs.  As  the  aircraft  path  length  is  increased,  the 
range  difference  changes  quadratically.  Due  to  the  fairly  short  wavelength  of  the  signal, 
this  range  difference  results  in  considerable  phase  rotation.  As  the  path  length  is  in¬ 
creased  on  an  unfocused  SAR,  there  is  a  limit  at  which  the  returns  received  at  the  end  of 
the  aperture  become  subtractive,  thus  degrading  the  overall  performance. 

SAR  applications  compensate  this  phase  degradation  by  focusing  the  signals  to  a 
specific  range.  An  unfocused  SAR  does  not  perform  this  correction  and  integrates  the 
return  signals  within  limits  where  degradation  is  negligible.  Stimson  [2]  calculates  the 
maximum  synthetic  aperture  length  for  the  unfocused  case  as: 

4^^=1.27Ir.  (2-7) 

By  inserting  (2.7)  in  (2.6),  the  azimuth  resolution  of  the  unfocused  SAR  can  be  repre¬ 
sented  as: 

‘I.s,,=0a4Ir.  (2.8) 


Equation  (2.8)  shows  that  for  the  unfocused  SAR  the  cross-range  resolution  is  range  de¬ 
pendent. 


2.  Focused  SAR 

The  limitations  on  the  SAR  aperture  size  presented  in  (2.7)  can  be  eliminated  by 
applying  an  appropriate  phase  correction  for  the  return  pulse  samples  according  to  the 
pulse  location  in  the  synthetic  aperture.  This  phase  correction  at  each  array  element  is 
given  by  [2]: 


Z  =  - 


^  2;^  ^ 


AR 


(2.9) 


where  rj'  is  the  distance  of  array  element  n  from  the  array  center.  The  use  of  these  phase 
corrections  allows  a  focused  SAR  to  integrate  samples  through  a  longer  synthetic  aper- 
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ture  length  L.  The  integration  time  can  be  selected  so  that  the  synthetic  aperture  length  L 
matches  the  azimuth  resolution  for  a  conventional  antenna  depicted  in  Equation  (2.3), 
such  that: 


L  =  RX/D. 


(2.10) 


The  resulting  synthetic  aperture  length  of  Equation  (2.10)  can  in  turn  be  used  in  (2.6), 
yielding  a  focused  SAR  azimuth  resolution  of: 


RA  D 


SAR 

2 

/D 

(2.11) 


where  D  is  the  aperture  size  of  the  physical  antenna. 

The  azimuth  resolution  achievable  with  a  focused  SAR  seems  contradictory  in  the 
case  of  conventional  antennas.  Fine  resolution  in  focused  SAR  requires  a  small  physical 
aperture  (large  beamwidth)  rather  than  a  large  one.  Furthermore,  the  resolution  for  the 
focused  case  is  independent  of  the  range  and  signal  wavelength. 

3.  SAR  Modes 

Several  operating  modes  can  be  identified  depending  on  the  type  of  target  and  the 
orientation  of  the  side  looking  antenna.  They  all  require  different  signal  processing  tech¬ 
niques  for  compensation,  integration  and  azimuth  compression.  These  modes  are  strip- 
map,  squint  mode,  spotlight  and  inverse  SAR. 

Strip-map  SAR  (Figure  2a)  is  designed  to  map  the  ground  within  a  strip  perpen¬ 
dicular  to  the  aircraft  velocity  vector.  The  antenna  is  fixed  to  one  side  of  the  aircraft  at  a 
right  angle  and  squinted  down  in  elevation.  The  aircraft  flies  over  a  ground  track  parallel 
to  the  strip  and  sweeps  a  broad  swath.  The  return  signal  is  processed  to  select  the  sam¬ 
ples  corresponding  to  the  range  where  the  strip  of  interest  is  located.  Both  negative  and 
positive  Doppler  ground  patches  are  processed. 

Squint  mode  SAR  (Figure  2b)  trains  the  antenna  to  a  fixed  position  forward  or 
backwards  from  its  strip-map  right  angle  position.  Here  the  ground  patches  have  either 
positive  or  negative  Doppler.  By  training  the  antenna  forward,  SAR  does  not  need  to 
wait  for  the  leading  edge  of  the  map  to  pass  by  a  full  array  length  to  obtain  a  map.  A 
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real-time  mapping  can  thus  be  made  over  territory  not  yet  flown  over.  A  drawback  of 
this  mode  is  that  the  azimuth  resolution  is  degraded,  such  that: 


d  = 


RX 


2-L-cos(^) 


(2.12) 


where  q)  is  the  azimuth  angle  of  the  antenna  measured  from  the  broadside  normal  [2]. 

Spotlight  SAR  (Figure  2c)  trains  the  antenna  in  azimuth  and  elevation  in  order  to 
continuously  point  to  a  specific  ground  patch  where  a  target  of  interest  is  located.  Due  to 
the  constant  steering  of  the  antenna,  its  beamwidth  does  not  limit  synthetic  aperture 
length  [3].  Therefore,  this  mode  offers  a  finer  azimuth  resolution  than  the  Strip-map  and 
Squint  modes. 


a.  b.  c. 


Figure  2  SAR  Modes:  a)  Strip-Map,  b)  Squint  and  c)  Spotlight  (After:[3]). 


Inverse  SAR  (ISAR)  is  designed  to  obtain  the  image  of  a  moving  target.  In  con¬ 
trast  to  SAR,  the  ISAR  synthetic  aperture  is  created  using  the  rotational  motion  of  the 
target  instead  of  using  the  motion  of  the  radar.  Simply  put,  ISAR  can  be  understood  as  a 
spotlight  SAR  mode,  where  the  target  moves  and  the  radar  is  stationary  as  shown  in 
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Figure  3.  The  end  result  of  ISAR  is  a  two-dimensional  image  in  range  and  Doppler  suit¬ 
able  for  target  recognition. 

The  ISAR  technique  has  probably  been  best  applied  to  ship  target  recognition  by 
airborne  radars,  where  the  pitch,  roll  and  yaw  movements  provide  an  excellent  way  to 
recognize  hull  shapes,  masts  and  other  superstructures.  Other  applications  range  from 
recognition  of  moving  ground  vehicles  to  the  imaging  of  satellites  from  ground  stations. 


Cross-range 


t 


Figure  3  ISAR  Concept  for  Two-Dimensional  Imaging  Based  on  Target  Rotation. 


C.  ISAR  PRINCIPLES 

ISAR  target  imagery  is  obtained  from  the  signal  processing  of  the  return  pulses 
collected  while  the  target  is  rotating.  The  Doppler  shift  produced  by  the  radial  rotational 
components  of  the  moving  target  is  extracted  through  the  cancellation  of  the  contribu¬ 
tions  from  the  velocities  of  both  the  radar  platform  and  the  target.  In  addition,  ISAR  re¬ 
lies  on  a  tracking  system  to  keep  the  target  inside  the  antenna  beam  during  the  data 
collection  period. 

The  main  data  processing  differences  between  SAR  and  ISAR  are  in  the  azimuth 
dimension.  In  the  range  dimension,  ISAR  pulse  compression  of  the  received  chirp  is 
similar  to  that  of  SAR.  In  the  azimuth  dimension,  however,  SAR  resolution  is  a  function 
of  the  physical  antenna  size  as  shown  in  (2.1 1).  In  ISAR,  it  is  a  function  of  integration 
time  (dwell  time)  [3]. 

ISAR  systems  often  use  a  target-based  processing  technique  known  as  prominent 
point  processing,  in  which  a  single  large  scatter  from  the  target  is  tracked  and  used  as  a 
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reference  signal  to  cancel  phase  errors  common  to  all  targets  [4],  The  calculation  of 
phase  corrections  relies  on  external  data  from  the  radar  platform,  such  as  flight  plan,  on¬ 
board  INS  or  a  ground  tracker. 

The  compensated  signals  are  resolved  in  range  (time  of  arrival)  and  in  cross-range 
by  their  Doppler  component  using  a  bank  of  Doppler  fdters.  In  order  to  achieve  a  high 
resolution  in  both  range  and  cross-range  (Doppler),  a  short-duration  pulse  with  a  wide¬ 
band  characteristic  is  needed.  In  other  words,  pulse  compression  techniques  are  required. 
The  chirp  and  the  stepped  frequency  pulse  waveforms  satisfy  these  characteristics. 

1.  ISAR  Resolution 

The  slant  range  resolution  for  a  chirp  waveform  is  given  by: 


(2.13) 


where  c  is  the  velocity  of  propagation  and  A  is  the  bandwidth  of  the  chirp  pulse.  Actual 
values  for  this  resolution  are  close  to  one  meter.  Appendix  C  includes  a  brief  description 
of  the  APS- 137  radar  with  some  of  its  resolution  parameters. 


If  a  target  rotates  at  a  constant  angular  rate  co ,  about  an  axis  perpendicular  to  the 
line  of  sight  (LOS),  a  scatterer  at  a  distance  h  from  the  axis  has  a  tangential  velocity 
ct)-h.  The  instantaneous  Doppler  frequency  shift  produced  by  the  rotational  component 
of  that  scatterer  is  given  by: 


Ieoh 


(2.14) 


where  fc  is  the  (center)  carrier  frequency  of  the  signal.  Solving  for  h  in  (2.14)  the  azi¬ 
muth  or  cross-range  resolution  can  be  expressed  as: 


^ a  ISAR  ~  ^  ^  ^  (2-15) 

2cd 

where  A/  is  width  of  the  Doppler  filters  or  Doppler  resolution,  which  in  turn  is  related  to 
the  coherent  integration  time  T/by: 
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(2.16) 


Therefore,  the  cross-range  resolution  can  now  be  expressed  in  terms  of  the  viewing  angle 
ij/  (aspect  angle  of  Figure  3)  as: 


X  _  X 

2  CO  Tj  2  y/ 


(2.17) 


In  conclusion,  ISAR  range  resolution,  as  with  most  classical  radars,  improves  by 
increasing  the  pulse  bandwidth.  Conversely,  ISAR  azimuth  resolution  improves  when 
the  target  has  a  fast  angular  rotation  co  or  when  the  integration  time  T/  is  increased. 

2.  ISAR  Range-Doppler  Image 

The  ISAR  image  is  the  result  of  a  collection  of  echo  signals  during  a  small  time 
interval  while  the  target  rotates  with  a  small  aspect  angle.  Each  scatter  returned  from  the 
target  has  a  reflectivity  value  that  is  plotted  on  a  range  versus  Doppler  frequency  rectan¬ 
gular  plot.  The  range  of  the  image  lies  along  the  radar  LOS  and  is  often  referred  to  as 
slant  range.  The  cross-range  dimension  of  the  image  lies  in  a  plane  perpendicular  to  the 
radar  LOS  and  can  be  referred  to  as  the  Doppler  dimension. 

Due  to  the  direct  proportionality  between  the  cross-range  location  of  a  scatterer 
and  its  radial  Doppler,  a  Range-Doppler  map  is  equivalent  to  a  two-dimensional  image  of 
the  target.  Although  the  image  cross-range  dimension  does  not  represent  true  distance, 
its  proportionality  makes  it  possible  to  use  it  for  target  recognition  purposes. 

Figure  4  shows  a  picture  of  a  ship  and  its  equivalent  ISAR  image  although  not 
taken  at  the  same  time.  Even  though  the  website  publicizing  the  image  does  not  clarify 
how  it  was  obtained,  it  can  be  estimated  that  the  ISAR  aircraft  was  directed  near  the  bow 
of  the  ship.  The  horizontal  distribution  of  scatterers  in  the  image  is  obtained  from  their 
time  delay  (range  imaging),  while  the  vertical  distribution  is  obtained  from  the  Doppler 
caused  by  the  pitch  movement  of  the  ship  in  the  direction  of  the  aircraft.  Thus,  the  masts 
of  the  ship  are  expected  to  generate  a  high  Doppler  due  to  their  longer  distance  to  the 
ship’s  rotation  center. 
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D.  RANGE  IMAGING 


1.  Pulse  Compression 

High  range  resolution  on  conventional  radars  can  be  achieved  using  a  very  short 
duration  pulse.  However,  the  bandwidth  of  a  short  pulse  is  large,  which  increases  the 
complexity  of  the  radar.  Among  the  various  limitations  to  implementing  such  a  pulse, 
Skolnik  [4]  mentions  greater  demands  on  the  signal  processing,  higher  probability  of  mu¬ 
tual  interference  with  other  radars,  less  accurate  Doppler  measurements  and  energy  limi¬ 
tations  to  achieve  long  ranges,  requiring  higher  peak  power. 


Figure  4  a)  Photograph  of  the  USS  Crocket,  b)  ISAR  Image  of  the  Same  Ship  (From: 

http://radar-www.nrl.navv.mil/Areas/ISAR/') 

The  benefits  of  the  long  duration  and  the  wide  band  spectrum  for  a  pulse  can  be 
combined  using  pulse  compression.  This  process  uses  a  long  pulse  of  width  T  and  ap¬ 
plies  a  frequency  or  phase  modulation  to  widen  the  spectrum  to  a  bandwidth  A .  The  re- 
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turning  pulse  at  the  receiver  is  compressed  to  a  shorter  pulsewidth  r  « 1  /  A ,  thus  improv¬ 
ing  range  resolution  by  a  factor  equal  to  the  pulse  compression  ratio  of  T !  r . 

Various  pulse  compression  techniques  have  been  developed  for  frequency  and 
phase  modulation.  Polarization  modulation  techniques  have  been  employed  on  an  ex¬ 
perimental  basis  [5].  Frequency  modulation  techniques  include  linear  frequency  modula¬ 
tion  (LFM  or  chirp)  and  stepped  frequency  waveforms.  Phase  modulation  techniques 
include  binary  phase  codes,  such  as  Barker  and  pseudo  random  codes,  and  poly-phase 
codes,  such  as  Frank,  PI,  P3  and  P4  techniques  [5]. 

The  DIS  countermeasure  was  designed  to  deceive  an  ISAR  using  LFM  pulse 
compression  [6].  The  LFM  chirp  is  produced  by  sweeping  the  carrier  frequency  over  the 
length  of  the  pulse  T  by  a  frequency  amount  of  A .  Its  waveform  is  described  in  detail  in 
Chapter  III. 

2.  Range  Imaging  Model 

Considering  the  one-dimensional  scenario  of  Figure  5  where  an  imaging  radar  il¬ 
luminates  a  one-dimensional  target,  the  range  and  reflectivity  of  its  Nr  target  scatterers 
can  be  represented  by  the  range  profile  function: 


Cr{x)  =  ^CT^S{x-x^)  (2.18) 

m=\ 

where  is  the  radar  cross-section  of  the  scatterer  located  at  the  range  x^.  Assuming 

that  a  chirp  signal  s(t)  is  used  as  the  transmitted  signal  in  a  lossless  media,  the  expected 
received  signal  after  sampling  has  the  form: 


IimTJ  = 


f 


2  X. 


A 


mT^  - 
V  c  J 


(2.19) 


where  Ts  is  the  analog-to-digital  converter  (ADC)  sampling  period  and  2x^/c  is  the 

round  trip  time  delay.  Using  the  ideal  target  function  in  (2.18),  the  received  signal  in 
(2.19)  can  be  written  as: 


I{mT^)  =  CXx)*s{mT^)  =  C^ 


^  c 

V  2  J 


(2.20) 
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where  the  *  denotes  convolution.  The  target  profile  function  Cr  is  also  linearly  trans¬ 
formed  from  the  range  domain  to  the  discrete  time  domain.  From  the  convolution  theo¬ 
rem,  Equation  (2.20)  can  be  expressed  as  the  product  of  discrete  Fourier  transforms 
(DFT)  so  that  the  range  profile  function  can  be  found  as: 


C 


c  m  T 


=  F 


-1 


F{5(m7;)}J  1S(A:)J 


(2.21) 


where  F  and  F“'  denote  respectively  a  DFT  and  an  inverse  DFT,  normally  implemented 
using  fast  Fourier  transforms  (FFT).  The  k  variable  is  an  index  for  frequency  resolution 
cells. 


3.  Range  Profile  from  Match  Filtering 

The  inverse  DFT  of  I(k)  /  S(k)  is  known  as  source  deconvolution  and  assumes  that 
S(k)  is  nonzero  at  all  frequencies.  In  practice,  radar  signals  are  band-pass  signals  and 
thus  range  imaging  is  based  on  matched  filtering.  Based  on  matched  filter  theory,  a  range 

I 

profile  ,  slightly  different  from  in  (2.21),  can  be  obtained  from  the  cross-correlation 
received  signal  with  the  transmitted  signal: 

C'^{mT^)  =  I{mT^)®S{mT^)  (2.22) 
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where  0  denotes  cross-correlation.  From  the  correlation  theorem,  the  DFT  of  Equation 
(2.22)  can  be  expressed  as  [8]: 

C[{k)  =  ¥{l{mT^)®s{mT^)]  =  l{k)  S\k)  (2.23) 


where  the  Fourier  transform  of  the  cross-correlation  is  equivalent  to  the  product  of  Fou¬ 
rier  transforms,  in  which  one  transformed  signal  is  complex  conjugated.  Replacing  the 
Fourier  domain  signal  I(k)  in  (2.23)  by  the  DFT  of  (2.19),  the  result  is: 


m-\ 


r 

mT^- 

V 


(2.24) 


From  the  inverse  DFT  of  (2.24)  the  target  range  profile  can  be  given  as: 


\^m=l 

JV,  / 


=y(T  s 

/  j  m 


n=\ 


2x. 


(2.25) 


mT^- 
V  c  J 


F-‘  Sik)\ 


The  target  profile  obtained  from  cross-correlation  in  (2.25)  when  transformed  to 
the  range  domain  has  a  similar  form  as  the  range  profile  in  (2.18)  except  for  the  term 
F“'  \S(k)f ,  which  Soumekh  [7]  defines  as  a  point-spread  function  (PSF).  Using  (2.25), 

the  target  range  profile  obtained  from  cross-correlation  can  be  related  to  the  ideal  decon¬ 
volution  range  profile  as: 

C:{mT,)  =  CjmT,)*¥-'  |5(i)f .  (2.26) 


From  the  previous  relation,  it  can  be  seen  that  the  PSF  term  determines  the  time 
(range)  resolvability  of  the  targets.  The  PSF  in  the  case  where  the  band-pass  spectrum 
possesses  an  amplitude  equal  to  one  has  the  form  of  the  sine  function  [7]  within  an  ap¬ 
proximate  time  interval  of  1/  A .  The  range  resolution  can  therefore  be  obtained  as: 


d 


r 


c  \  _c  T 
2 


(2.27) 


which  agrees  with  the  result  presented  in  (2.13). 
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The  previously  described  matched  filtering  technique  when  applied  with  chirp 
signals  is  essentially  a  pulse  compression.  In  the  following  sections,  it  will  be  stated  as 
range  compression  to  differentiate  it  from  the  azimuth  or  Doppler  compression. 


E.  ISAR  SIGNAL  PROCESSING  EOR  A  CHIRP  PULSE  RADAR 

1.  Range  Compression 

The  range  compression  process  usually  occurs  after  the  received  signal  is  detected 
in  quadrature,  sampled  and  digitized.  The  resulting  digital  signal  must  retain  both  magni¬ 
tude  and  phase  of  each  sample  for  subsequent  coherent  processing.  A  target  range  profile 
can  be  obtained  from  each  pulse  of  the  returned  signal  using  the  expression  in  (2.22), 
where  the  multiplication  of  Fourier  transforms  is  equivalent  to  a  signal  cross-correlation. 

Figure  6  outlines  a  block  diagram  of  the  fast  correlation  that  implements  the  range 
compression  process.  The  ISAR  received  signal  is  down-converted  in  frequency  and 
processed  with  a  Fast  Fourier  Transform  (FFT).  A  coherent  reference  signal  from  the 
transmitted  chirp  is  also  processed  with  an  FFT  and  then  converted  into  its  complex  con¬ 
jugate  equivalent.  Both  signals  are  multiplied  and  processed  with  an  inverse  Fourier 
transform  to  obtain  the  target  range  profile. 

2.  Azimuth  Compression 

The  range  compression  process  produces  one  range  profile  for  each  of  the  Np  re¬ 
turned  pulses  that  are  collected.  These  Np  range  profiles  can  be  arranged  in  a  matrix 
where  each  row  corresponds  to  a  different  returned  pulse  and  each  column  is  comprised 
of  profile  elements  belonging  to  the  same  range  bin.  This  matrix  can  be  expressed  as: 


C'r{mT^,n) 


for  pulse 

Cvp-l  forpulseAfp-1 

C'rlim  ^  for  pulse  2 
for  pulse  1 


(2.28) 


where  (t)  are  range  profiles  for  different  pulses  as  defined  in  (2.22). 
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Figure  6  Range  Compression  Using  Fast  Correlation. 

The  process  of  azimuth  compression  is  basically  the  FFT  of  the  time  history  of 
range  profiles.  At  each  range  cell  location  (columns),  the  time  dimension  expressed  in 
terms  of  pulse  indices  is  transformed  to  the  frequency  dimension  in  terms  of  Doppler 
cells.  Essentially,  the  FFT  of  each  range  profile  matrix  in  the  column  dimension  yields 
the  ISAR  range-Doppler  image  [8].  Thus,  the  ISAR  image  can  be  expressed  as: 

C(m,A:)  =  FFT{c:(in7;,n)}=  — (2.29) 

The  creation  of  the  matrix  C^{mT^,n)  by  stacking  range  profiles  from  different 
pulses  relies  on  a  precise  range  alignment.  Thus,  range  tracking  systems  are  used  to  con¬ 
tinually  align  the  range  profiles  («T^)  by  adjusting  time  delays.  These  systems  lock 

on  a  prominent  scatterer  point  to  set  a  range  trigger  delay  reference  to  compensate  the 
rest  of  the  scatterers  of  the  range  profile. 

In  the  azimuth  dimension,  the  oscillatory  movement  of  ship  targets  varies  the  ap¬ 
parent  cross-range  dimensions  depending  on  the  rotational  velocity  during  the  data  inte¬ 
gration  time.  Figure  7  shows  different  ISAR  images  obtained  from  the  same  target  at 
different  times.  Image  1  shows  a  ship  pitching  away  from  the  radar  with  its  masts  pro¬ 
ducing  a  negative  Doppler.  Images  2  and  3  show  that  the  data  was  collected  during  a  pe- 
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riod  of  slow  ship  rotation.  Thus,  in  accordance  with  Equation 
Error!  Reference  source  not  found.,  their  azimuth  resolution  is  poor.  Images  4,  5  and  6 
were  generated  at  an  integration  period  of  fast  ship  rotation.  Therefore,  the  Doppler  shift 
produced  by  the  masts  stands  out  and  gives  a  ship  silhouette  less  difficult  to  recognize. 


Image  3  Image  4 


Image  5 


Image  6 


Figure  7  Sequence  of  Six  ISAR  Images  (From:  http://airbome.nrl.naw.mi1/I 


F.  CHAPTER  SUMMARY 

The  concept  of  SAR  and  ISAR  were  presented  and  explained  in  terms  of  their 
resolution  advantage.  Next,  the  range  imaging  concept  was  presented  in  a  one¬ 
dimensional  analysis,  where  range  compression  was  described  in  terms  of  match  filtering 
by  implementing  fast  correlation.  Finally,  the  ISAR  processing  technique  of  azimuth 
compression  and  the  generation  of  the  Range  Doppler  profile  were  described.  The  next 
chapter  presents  the  DIS  concept,  for  which  a  primary  level  of  range  profiling  and  ISAR 
processing  is  assumed. 
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III.  THE  DIS  PROCESS 


This  chapter  presents  the  DIS  eoncept  and  its  signal  proeessing  theory  as  first  de¬ 
scribed  in  [10].  Further  work  on  this  concept  has  been  discussed  in  referenees  [6],  [9] 
and  [11-18],  Following  the  theoretical  description  of  the  proeess,  the  DIS  hardware  is 
deseribed  at  a  system  level. 

A.  OVERVIEW 

The  DIS  was  designed  for  eleetronic  attaek  against  linear  frequeney  modulation 
(LFM)  ISARs  by  emulating  the  seattering  of  a  false  target.  In  order  to  effectively  deceive 
an  ISAR,  the  DIS  must  be  able  to  generate  the  same  seattering  effect  that  a  real  target 
would  produce  on  the  returning  radar  signal.  The  false  target  should  be  of  a  realistic  size 
and  its  Doppler  effects  should  reflect  a  feasible  target  rotational  motion.  The  transmitted 
power  of  the  DIS  deeeiving  signal  should  be  high  enough  to  eover  the  true  scattering  re¬ 
turn  produced  by  the  DIS  platform  and  low  enough  to  avoid  the  saturation  of  the  ISAR 
reeeiver,  and  thus  revealing  the  presenee  of  jamming.  Figure  1  of  Chapter  I  shows  a 
block  diagram  summarizing  the  DIS  operational  eoncept  to  deeeive  an  ISAR  platform. 

To  produee  the  scattering  effeet  of  a  false  target,  the  DIS  is  provided  with  a  set  of 
pulses  that  are  intercepted  and  down-eonverted  in  frequeney  from  an  ISAR.  The  DIS 
signal  proeessing  ineludes  amplitude  and  phase  modulation  to  recreate,  on  a  pulse-to- 
pulse  basis,  the  amplitude  ehange  and  Doppler  profile  caused  by  the  reflecting  surfaces  of 
the  target  at  each  location  in  range  and  eross-range.  The  DIS  output  is  upconverted  to  the 
earrier  frequency  and  retransmitted  back  to  the  radar.  A  bloek  diagram  summarizing  the 
signal  proeessing  oecurring  in  the  DIS  eircuit  is  shown  in  Figure  8. 
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Figure  8  Bloek  Diagram  of  the  DIS. 
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The  ISAR  being  targeted  receives  the  DIS  modulated  signal  and  passes  it  through 
a  matched  filter  that  compensates  for  platform  motion.  The  range  profiles  are  obtained 
for  each  pulse  and  collected  in  a  matrix  to  build  a  time  history  of  the  target  rotational 
movement.  An  TFT  of  this  time  history  provides  a  Doppler  mapping  of  the  range  pro¬ 
files,  thus  creating  a  two-dimensional  profile  or  image.  The  image  is  enhanced,  com¬ 
pared  with  other  images,  analyzed  for  features  that  can  be  extracted,  and  processed  for 
automatic  target  recognition.  Figure  9  shows  a  block  diagram  summarizing  the  ISAR 
processing  of  the  DIS  signal  and  the  generation  of  the  false  target  range-Doppler  image. 
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Figure  9  Block  Diagram  of  the  ISAR  Reception  and  Image  Generation. 

This  chapter  describes  the  signal  processing  occurring  in  the  DIS  and  its  peripher¬ 
als.  The  description  is  divided  into  the  generation  of  modulation  coefficients,  the  sam¬ 
pling  of  the  ISAR  signal  and  DIS  signal  processing.  It  also  includes  a  simplified  model 
of  the  ISAR  receiver  signal  processing  to  describe  the  basics  of  the  target  range  Doppler 
image  formation.  Finally,  a  brief  description  of  the  DIS  architecture  explains  how  the 
signal  processing  is  implemented  in  hardware.  The  equations  presented  in  this  chapter 
form  the  basis  of  the  simulation  software  used  to  evaluate  the  effectiveness  of  the  DIS 
design  as  an  ISAR  countermeasure.  These  simulation  results  are  presented  in  Chapter  V. 

B.  GENERATION  OF  THE  DIS  MODULATION  COEFFICIENTS 

In  order  to  synthesize  echoes  corresponding  to  a  false  target,  the  proper  phase  and 
amplitude  modulation  must  be  performed  for  each  range  bin  to  synthesize  the  Doppler 
effect  due  to  target  rotation  and  the  expected  radar  cross-section  value  at  each  scatterer 
location  in  the  range-Doppler  plane. 
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A  false  target  range-Doppler  profile  is  built  by  creating  a  two-dimensional  matrix 
of  complex  scattering  elements.  In  the  case  of  a  ship  target,  the  silhouette  can  be  laid  out 
in  an  Nr  by  Nd  matrix,  where  Nr  is  the  number  of  the  resolution  cells  in  the  range  axis  and 
Nd  is  the  number  of  Doppler  resolution  cells  in  the  cross-range  axis.  For  simplicity,  the 
range-Doppler  profde  is  usually  set  with  the  ship’s  longitudinal  axis  oriented  along  the 
range  dimension  and  the  ship’s  vertical  axis  oriented  along  the  cross-range  axis  as  shown 
in  Figure  10. 


Range  resolution  cells  (r) 


Figure  10  Range-Doppler  Profile  Template  for  a  False  Target  (Cell  Doppler  f(r,d)  in  Hz). 


The  amplitude  of  each  scatterer  in  the  range-Doppler  profde  depends  on  the  radar 
cross-section  at  that  location.  The  phase  of  each  scatterer  depends  on  the  Doppler  shift 
due  to  the  velocity  component  along  the  LOS  for  each  target  location.  The  Doppler  shift 
produces  a  time  dependent  phase  shift  and,  thus,  is  different  for  each  pulse  arrival.  The 
mathematical  representation  for  each  target  return  as  a  function  of  the  transmitted  pulse 
index  n  and  the  scattering  location  in  the  false  target  profile  can  be  expressed  in  a  tri¬ 
dimensional  matrix  as: 

T{r,d,n)  =  A{r,d) 

where  r  and  d  are  target  location  indices  is  terms  of  range  and  Doppler  resolution  cells,  A 
is  the  amplitude  of  the  scatterer  due  to  the  radar  cross-section  at  the  (r,d)  location  and 
f(r,d)  is  the  Doppler  frequency  shift  that  the  scatterer  exerts  on  the  ISAR  signal. 
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A  simplified  expression  for  the  return  signal  can  be  obtained  by  adding  the  com¬ 
plex  scatterers  along  the  azimuth  (Doppler)  dimension.  The  resulting  superposition  ex¬ 
pression  is  given  by: 

T'[r,n)  =  ^T[r,d,n).  (3.2) 

d=\ 

Equation  (3.2)  represents  a  return  signal  for  each  pulse  n  at  each  target  range  resolution 
bin  r.  This  signal  can  be  used  to  obtain  both  the  phase  and  the  magnitude  modulation  co¬ 
efficients  required  by  the  DIS. 

The  phase  of  the  return  signal  is  calculated  for  each  range  cell  and  pulse  index 
combination  using  the  angle  function,  which  is  the  arctan  function  of  the  division  be¬ 
tween  the  imaginary  and  the  real  components,  using  their  signs  to  solve  the  ±7:  ambigu¬ 
ity: 

(jdj  (r,«)  =  angle{real\T'{r,n^^,  imag{T'[r,n)]y  (3.3) 

The  phase  modulation  effect  produced  by  a  target  reflection  must  be  understood  as  a 
phase  increment  on  a  pulse-by-pulse  basis.  For  that  reason,  the  phase  of  the  returned 
pulse  is  subtracted  from  the  phase  of  the  previous  pulse  at  each  range  bin  r  and  added  to 
the  previous  increments  to  obtain  a  cumulative  sum  of  increments: 

The  cumulative  phase  increment  describes  the  ideal  phase  modulation  coefficient  needed 
for  image  synthesis.  The  phase  increment  from  Equation  (3.4)  is  scaled  and  truncated  to 

integer  values  from  0  to  2*^  - 1 .  The  final  expression  for  the  kp-h\i  phase  coefficient  fed 
to  the  DIS  at  each  range  bin  r  and  pulse  n  is  given  by: 


where  [  J  is  the  floor  integer  rounding  function,  defined  as  the  greatest  integer  less  than 
or  equal  to  the  argument. 
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The  amplitude  modulation  effect  in  the  return  signal  can  be  synthesized  by  multi¬ 
plying  the  intercepted  signal  by  a  gain  value  equal  to  the  magnitude  of  the  sum  of  com¬ 
plex  scatterers  at  the  same  range-bin/pulse  combination.  Thus,  the  first  step  to  obtain  the 
gain  coefficients  is  to  normalize  the  magnitude  of  T\r,n)  in  Equation  (3.2)  so  that  all 
magnitudes  have  a  range  of  values  between  0  and  1 : 


TN{r,n) 


T'{r,n) 
max  T'  {r,n) 


(3.6) 


The  DIS  amplitude  modulation  is  implemented  by  shifting  the  I  and  Q  bit  formats 
by  a  specific  number  of  binary  places.  This  simple  procedure  avoids  the  use  of  multipli¬ 
ers  and  reduces  processing  time.  The  amount  of  bit  shift  g(r,n)  is  set  by  a  control  code  of 

kg  bits.  Although  the  available  control  codes  are  2*^^ ,  as  it  will  be  explained  in  the  next 
chapter,  some  of  them  may  produce  the  same  amount  of  shift.  Since  the  available  gains 
are  power  of  two  values,  the  quantization  of  the  normalized  magnitudes  in  (3.6)  to  assign 
these  gains,  should  also  be  distributed  exponentially.  For  a  gain  shifter  controlled  with  a 
4-bit  format,  the  16  available  codes  only  produce  11  unique  shifts.  One  possible  gain 
quantization  scheme  for  this  case  can  be  implemented  as  shown  in  Table  1: 


Magnitude 

Ranges 

Number  of 
Shifts  ^(r,n) 

Effective  Gain 

0.8  to  1.0 

10 

1024 

0.4  to  0.8 

9 

512 

0.2  to  0.4 

8 

256 

0.1  to  0.2 

7 

128 

0.05  to  0.1 

6 

64 

0.025  to  0.05 

5 

32 

0.0125  to  0.025 

4 

16 

0.0063  to  0.0125 

3 

8 

0.0032  to  0.0063 

2 

4 

0.0016  to  0.0032 

1 

2 

Oto  0.0016 

0 

1 

Table  1  Gain  Modulation  Quantization  Scheme. 


C.  INTERCEPTION  AND  SAMPLING  OF  THE  ISAR  SIGNAL 

An  LFM  ISAR  pulse  with  constant  amplitude  and  upward  frequency  sweep 
(chirp)  can  be  expressed  as: 


r{t)  =  rect 


\T^J 


(3.7) 
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where  t  is  the  uncompressed  pulsewidth  in  seconds,  A  is  the  modulation  bandwidth  in 
Hertz  and  fc  is  the  transmitted  carrier  frequency,  also  in  Hertz.  The  envelope  of  the  chirp 
pulse  is  given  by  the  rectangular  function  rect(t/r),  where  there  is  a  constant  amplitude 

of  one  for  values  of  the  argument  tjr  between  -0.5  and  +0.5.  Figure  1 1  shows  the  pulse 
envelope,  the  modulating  frequency  used  at  each  time  and  the  resulting  LFM  pulse. 


Instantaneous 

Frequency 


Signal 

Amplitude 


Figure  1 1  Construction  of  an  LFM  Chirp  Pulse  (From:  [6]). 


For  an  ISAR  mounted  on  a  moving  platform,  the  signal  intercepted  by  a  moving 
target  has  an  instantaneous  phase  advance  at  time  t  given  by: 

0{t)  =  ^^(R  +  vt  +  ae  ll)  (3.8) 

where  R,  v  and  a  are  the  initial  range,  initial  relative  velocity  and  relative  acceleration 
between  the  ISAR  platform  and  the  target.  The  corresponding  frequency  variation  due  to 
this  phase  change  is  the  first  derivative  of  the  previous  equation  given  by: 

/(0  =  ^(v+a()-  (3-9) 

c 
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Therefore,  the  LFM  ISAR  signal  in  Equation  (3.7),  modified  for  the  case  of  a  moving 
platform  and  a  moving  target,  can  be  expressed  as: 


r{t) 


=  rect 


^  t  \  '/2r 


(3.10) 


The  Doppler  frequency  component  f(t),  considering  realistic  platform  velocities 
and  ISAR  frequencies  in  the  X  band  (8-12  GHz),  is  on  the  order  of  tens  of  Hertz.  Since 
the  modulation  bandwidth  is  usually  on  the  order  of  hundreds  of  MHz,  it  can  be  assumed 
that  the  intercepted  chirp  pulse  has  a  negligible  Doppler  (f~0)  or  constant  phase  change. 
From  (3.10),  the  corresponding  complex  envelope  obtained  after  down  converting  the 
intercepted  chirp  is  given  by: 


s(0 


=  rect 


f-1 

UJ 


I T 


(3.11) 


This  complex  envelope  is  stored  using  a  digital  RF  memory  (DRFM).  This 
equipment  permits  the  indefinite  and  non-degraded  storage  of  radio  frequency  signals  in 
digital  form.  The  stored  signals  can  be  transcribed  back  into  an  exact  replication  of  the 
original  signal,  which  allows  coherent  transponding  of  radar  signals  that  are  very  suitable 
for  EA. 


A  DRFM  is  fed  with  the  intercepted  ISAR  signal  and  replicates  it  in  terms  of 
phase  samples  or  complex  envelope  amplitude  samples  (I  and  Q),  at  sample  times  given 
by: 

t  =  —  +  nPRl  (3.12) 

fs 


where  m  is  the  index  of  the  intrapulse  samples,  n  is  the  index  of  each  pulse,  ^  is  the  ADC 
sampling  frequency  and  PRl  is  the  ISAR  pulse  repetition  interval  in  seconds.  Although 
many  imaging  radars  use  variable  PRI,  the  PRI  stays  constant  within  one  integration 
frame  time.  Substituting  the  sample  time  from  Equation  (3.12)  into  (3.11),  the  sampled 
phase  of  the  complex  envelope  can  be  expressed  as: 


K  A 

T 


{  m 


Kfs 


Y 

+  nPRI 

J 


(3.13) 


29 


The  phase  signal  in  Equation  (3.13),  as  generated  by  a  kp-hit phase  sampling  DRFM,  ean 
be  expressed  as: 


In 


(3.14) 


which  is  obtained  by  rescaling  the  complex  envelope  phase  and  truncating  it  to  the  inte¬ 
ger  value.  The  kp-hit  DRFM  phase  expression  of  (3.14)  gives  2*^  integer  values,  ranging 

from  0  to  2^^  - 1 .  Therefore,  the  DRFM  phase  samples  have  the  same  kp-hit  format  as 
the  phase  modulation  coefficients  from  Equation  (3.5).  The  phase  resolution  for  this 
DRFM  output  is  given  by  the  number  of  phase  quantization  levels  {in^jl’"^  . 

For  the  case  of  a  complex  amplitude-sampling  DRFM,  the  intercepted  signal  is 
stored  in  terms  of  the  orthogonal  components  of  the  complex  envelope: 


Do^(m,«)  +  7  Dog(/n,n)  =  cos(<zio(/n,«))  +  7  sin((t>^[m,n)).  (3.15) 

The  complex  envelope,  as  generated  by  the  two  output  channels  of  a  A:^-bit  amplitude¬ 
sampling  DRFM,  can  be  expressed  as: 


Dj{m,n)=  D^j{m,n)  2 


D^{m,n)=  D^Q{m,n)  2 


't„-i 


(3.16) 


Either  the  phase  or  amplitude  DRFM  signals  can  be  fed  into  the  DIS  chip  in  order  to  syn¬ 
thesize  the  effects  of  a  reflecting  target. 


D.  THE  DIS  SIGNAL  PROCESSING 

The  DIS  is  a  parallel  array  of  identical  complex  digital  modulators  and  each  one 
synthesizes  the  scattering  effects  for  one  range  bin.  The  input  coefficients  necessary  to 
synthesize  an  image  from  the  DRFM  samples  consist  of  the  target  dependent  phase  and 
gain  coefficients  and  the  desired  target  extent  provided  by  a  microprocessor.  Figure  12 
presents  a  block  diagram  of  the  DIS  along  with  the  external  components.  The  wideband 
chirp  signal  that  feeds  the  block  diagram  is  described  by  Equation  (3.10).  After  down- 
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conversion,  the  complex  envelope  components  are  represented  by  Equation  (3.11).  The 
phase  or  amplitude  (1  and  Q)  DRFM  samples  are  described  by  Equations  (3.14)  and 
(3.16),  respectively.  For  the  case  of  amplitude  DRFM  samples,  they  are  processed  by  a 
phase  conversion  circuit. 


l(jn,ri) 


Figure  12  DIS  Block  Diagram  (From:  [6]). 

For  the  cases  where  an  amplitude-sampling  DRFM  is  used  to  feed  the  DIS,  it  is 
necessary  to  convert  the  A:^-bit  I  and  Q  samples  to  ^^-bit  phase  samples.  This  arctan  op¬ 
eration  is  carried  out  using  the  Coordinate  Rotation  Digital  Computer  (CORDIC)  algo¬ 
rithm.  CORDIC  rotates  the  phase  of  /  +  jQ  in  the  complex  plane  by  multiplying  powers 
of  two  values  using  shift-add  algorithms.  These  incremental  phase  rotations  make  it  pos¬ 
sible  to  calculate  trigonometric  functions  such  as  arctan  without  the  need  of  multipliers 
[15]. 
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Each  phase  sample  strobed  from  the  DRFM  <j){m,n)  is  fed  in  parallel  to  all  the 

DIS  complex  modulators  (range  bins)  and  are  first  added  with  the  false  target  phase  coef¬ 
ficients  corresponding  to  each  range  bin.  The  phase  coefficient  present  at  each  range  bin 
changes  for  every  new  ISAR  intercepted  pulse.  That  is,  at  the  r*  array  branch,  the  w* 
sample  of  the  ISAR  pulse  is  added  with  the  phase  coefficient  for  that  range-bin/pulse 
combination,  giving  a  phase  output  of: 

^{r,m,n)  =  (^{m,n)  +  (p.^^{r,n).  (3.17) 

The  modulated  phase  obtained  in  (3.17)  at  each  range  bin  is  used  to  construct  a 
complex  signal  of  unit  amplitude  given  by: 

L[r,m,n)  =  cos{^^{r,m,n  ))  +  7  sin(<^(  r,m,n  )).  (3.18) 

The  sine  and  cosine  operations  for  the  2*^^  different  phase  values  are  implemented  using  a 
look-up  table  (LUT)  with  a  depth  of  '2^’’  and  a  resolution  of  ki  bits.  The  LUT  output  pre¬ 
sented  in  (3.18)  has  two  output  channels:  I  (cosine)  and  Q  (sine),  each  with  amplitude 
values  ranging  from  - 1  to  1 . 

The  I  and  Q  output  channels  of  the  LUT  are  amplitude  modulated  with  the  gain 
coefficients  for  each  range-bin/pulse  combination.  The  multiplication  by  the  power  of 
two  gain  coefficients  is  implemented  by  shifting  bit  positions  to  the  left.  The  gain  modu¬ 
lated  channels  can  be  expressed  as: 

S{r,m,n)  =  L{r,m,n)  =  (3.19) 

where  g(r,n)  is  the  gain  coefficient. 

Following  the  gain  modulation,  all  further  range-bin  processors  are  summed  with 
the  current  processor  and  passed  forward.  The  final  DIS  output  is  taken  from  range  bin 
one  with  complex  samples  per  pulse,  where  r  is  the  pulsewidth.  The  DIS 

output  expression  for  the  w*  sample  of  the  w*  pulse  is  given  by: 

I{m,n)  =^S{r,m,n)='^  (3.20) 

r=0  r=0 
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The  DIS  output  samples  are  fed  to  a  digital-to-analog  converter  (DAC),  coherently  modu¬ 
lated  with  the  same  local  oscillator  signal  used  for  reception  and  transmitted  back  to  the 
victim  ISAR. 


E.  THE  ISAR  RECEPTION  AND  IMAGE  GENERATION 

In  order  to  evaluate  the  proper  generation  of  a  false  target,  it  is  necessary  to  model 
the  process  occurring  in  the  ISAR  receiver.  ISAR  signal  processing  is  described  in  Chap¬ 
ter  II,  but  the  model  used  in  the  simulation  is  simplified  to  a  range  compression  stage  fol¬ 
lowed  by  a  Doppler  (cross-range)  compression.  This  section  describes  the  range  and 
azimuth  compression  processes  in  terms  of  the  mathematical  equations  applied  in  the 
simulation  software. 


The  ISAR  range  compression  is  implemented  using  fast  correlation  between  the 
received  DIS  signal  containing  false  target  echoes  and  a  coherent  reference  signal.  The 
correlation  process  is  implemented  as  a  product  of  FFT  signals  as  shown  in  Figure  6. 
From  equations  (3.11)  and  (3.12),  the  ISAR  reference  function  is  described  by  a  time- 
periodic  chirp-pulse  train  with  samples  given  by: 


'{m,n)  =  rect 


j^\  +n  PRI 


(3.21) 


The  FFT  of  (3.21)  determines  the  transfer  function  of  the  pulse  compression  network  for 
each  pulse  n  at  the  Doppler  fdter  t. 

=  (3.22) 

A  m=0 


where  N  represents  the  number  of  FFT  points.  For  better  FFT  resolution  N  can  be  in¬ 
creased  by  padding  zeroes  at  the  end  of  the  reference  signal  vector  in  (3.21). 

The  DIS  signal  received  by  the  ISAR  after  down-conversion  to  base  band  corre¬ 
sponds  to  the  expression  in  (3.20).  The  FFT  is  given  by: 

/(/t,n)  =  FFT{/(m,n)}  =  — (3.23) 

A  „=o 


33 


Multiplying  the  complex  conjugate  of  the  transfer  function  in  (3.22)  by  the  frequency 
signal  in  (3.23)  and  applying  an  inverse  FFT,  yields  the  complex  range  profde  of  the  false 
target  echoes  corresponding  to  the  n*  transmitted  pulse: 

{m,n)  =  FFT"'  ^l{k,n)S*  =  —'^l{k,n)S*  (3.24) 

^  k=Q 

The  range  profde  obtained  for  each  transmitted  pulse  after  range  compression  is  a  collec¬ 
tion  of  complex  scatterers  grouped  at  each  range  bin. 

The  Doppler  profile  at  each  range  bin  is  obtained  from  the  magnitude  of  the  FFT 
of  the  complex  scatterers  in  the  time  dimension  (pulse  indices).  Using  Equation  (2.29), 
the  result  of  this  azimuth  compression  for  all  range  bins  gives  the  range-Doppler  profde 
of  the  false  target  as: 

I  I  1 

C(w,A:)  =  |FFT{Q(m,n)}|=  —  (3.25) 

N  „=Q 

The  range-Doppler  profde  obtained  in  Equation  (3.25)  is  the  end  result  of  the  DIS  EA 
effort.  This  image  profde  should  deceive  the  ISAR  operator  and  the  supporting  target 
recognition  algorithms. 

F.  ARCHITECTURE  OF  THE  DIS  MODULE 

The  DIS  module,  as  presented  in  Figure  12,  is  a  parallel  array  of  Nr  identical 
modulators  controlled  by  an  off-chip  microprocessor.  Each  modulator  represents  the  ma¬ 
jor  reflective  surfaces  at  each  range  bin  of  the  target  profde  and  thus  will  be  referred  to  as 
a  range  bin.  The  modulation  at  each  range  bin  is  performed  through  a  pipelined  process 
composed  of  a  phase  modulation,  a  polar-to-rectangular  conversion  of  complex  samples 
with  a  sine/cosine  table,  an  amplitude  modulation  and  a  cumulative  sum.  The  hardware 
implementation  of  a  range  bin  depicted  in  Figure  13  is  briefly  described  in  this  section. 
The  current  DIS  chip  design  has  Nr  =  5\2  range  bins.  Additional  chips  can  be  cascaded 
to  synthesize  larger  targets. 

1.  DIS  Phase  Modulation 

The  DIS  input  signal  is  a  series  of  DRFM  generated  phase  samples  (or  I  and  Q 
samples  converted  to  phase  by  a  CORDIC  circuit)  from  Np  chirp  pulses.  Each  phase 
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sample  is  fed  simultaneously  to  all  Nr  range  bins.  These  samples  will  be  modulated 
within  each  range  bin  with  a  different  phase  modulation  coefficient  updated  for  every 
new  pulse.  This  pulse-to-pulse  phase  modulation  (or  phase  rotation)  synthesizes  the  de¬ 
sired  Doppler  frequency  profile  for  each  range  bin. 


1.  Phase  increment 

2.  Phase  data  from  DRFM 

3.  Pre-load  register 

4.  Pipeline  register 

5.  5-bit  adder 

6.  Look  Up  Table 

7.  16-bit  adder 

8.  Gain  shifter 

9.  Input  from  previous  bin 

10.  Control  Logic 

1 1 .  Output  to  next  bin 


1  2 


Figure  13  Hardware  Diagram  for  a  DIS  Range  Bin  (After:[14]). 

As  described  in  Equations  (3.5)  and  (3.14),  for  kq=5,  both  phase  signals  have  a  5- 
bit  format  from  00000  to  11111  representing  32  phase  levels  with  increments  of  1 1.25°. 
Table  2  lists  the  32  possible  values  of  this  phase  format.  The  phase  summation  is  imple- 


35 


merited  with  a  5-bit  phase  adder  as  shown  in  Figure  14.  A  eontrol  microprocessor  pro¬ 
grams  the  phase  coefficients  which  are  latched  separately  into  the  corresponding  range 
bins  at  each  new  PRI. 


5-bit 

format 

Phase 

value 

5-bit 

format 

Phase 

value 

5-bit 

format 

Phase 

value 

5-bit 

format 

Phase 

value 

00000 

0° 

01000 

90° 

10000 

180° 

11000 

270° 

00001 

11.25° 

01001 

102.5° 

10001 

191.25° 

11001 

281.25° 

00010 

22.5° 

01010 

112.5° 

10010 

202.5° 

11010 

292.5° 

00011 

33.75 

01011 

123.75° 

10011 

213.75° 

non 

303.75° 

00100 

45° 

01100 

135° 

10100 

225° 

11100 

315° 

00101 

56.25° 

01101 

146.25° 

10101 

236.25° 

11101 

326.25° 

00110 

67.5° 

OHIO 

157.5° 

10110 

247.5° 

lino 

337.5° 

00111 

78.75° 

01111 

168.75° 

10111 

258.75° 

11111 

348.75° 

Table  2  Equivalence  Table  for  the  DIS  5-Bit  Phase  Format. 


B4  B3  B2  B1  BO  A4  A3  A2  A1  AO 

5-bit  Adder 

S4  S3  S2  SI  SO 


00000 

Figure  14  Phase  Adder  Diagram. 


2.  DIS  Cosine  and  Sine  LUT 

The  output  of  the  phase  modulation  adder  is  a  5-bit  format  representing  the  phase 
of  a  complex  sample  with  amplitude  one.  The  next  step  is  to  convert  this  complex  ex¬ 
pression  from  polar  to  rectangular  format  as  described  in  Equation  (3.18).  The  conver¬ 
sion  is  implemented  using  the  5-bit  phase  adder  outputs  as  addresses  of  a  cosine/  sine 
LUT  with  a  dual  channel  8-bit  output  format.  An  input/output  diagram  of  the  LUT  is 
shown  in  Figure  15. 

The  sine/cosine  LUT  has  32  (2^)  possible  pairs  of  output  amplitudes.  Each  output 
pair  is  composed  of  the  1  and  Q  rectangular  component  of  a  unitary  amplitude  complex 
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sample.  The  I  and  Q  output  channel  amplitudes  are  each  represented  by  an  8-bit  two’s 
complement  format.  The  format  has  one  sign  bit  and  seven  decimal  bits  from  2"'  to  2"^ 
as  shown  in  Table  3.  Chapter  IV  uses  additional  formats  to  study  the  impact  of  replacing 
this  8-bit  resolution  LUT  with  tables  of  resolution  less  than  8-bits. 
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^3‘> - 

Ei> - 
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COSO 

COSl 


OO 

COSJ 

C0S4 

C0S5 


C0S2  ] - |cOsf)> 

--{cosL)> 
--{cosr> 

-(cosL)> 


C0S6 

C0S7 

SIHO  j-- 
smi 


CO  SIH4 
SINS 


-{cosr)> 


-{siNr)> 
SIH2  j - (siNf  > 


— (sraf)> 
— (sraL> 


SIH7  j — 
_ I 


— {sra^) 
-^sraL) 


Figure  15  LUT  Diagram  (From:  [14]). 


Bit  Position 

8 

7 

6 

5 

4 

3 

2 

1 

Value 

Sign 

2'^ 

2'^ 

2'^ 

2-' 

2'^ 

2'® 

2-’’ 

Table  3  Bit  Format  of  the  LUT  Output  Amplitudes. 

The  absence  of  an  integer  bit  in  the  LUT  format  of  Table  3  does  not  allow  the  rep¬ 
resentation  of  a  value  of  ‘1.0’,  needed  for  the  phase  inputs  of  0  and  90  degrees.  The  am¬ 
plitude  values  represented  with  this  format  are  limited  to  a  range  from  - 1  to  1  -  2"’ , 
where  the  upper  limit  can  be  obtained  as: 

^  =  2"‘+2""+2"'+2-^  +  2"'+2"®+2“"  =1-2“"  =0.9921875  .  (3.26) 

In  order  to  provide  a  solution  for  the  cosine  and  sine  values  of  0  and  90  degrees,  the  sine 
and  cosine  LUT  values  are  scaled  down  by  a  factor  of  s,  equivalent  to  the  upper  limit 
found  in  (3.26).  Using  the  scaled  LUT  values,  the  sine  output  for  the  phase  angle  equiva¬ 
lent  to  90  degrees  is  given  by: 

sin  (8-1 1.25°)  =  1 «  4,^  (OlOOO)  =  .s  (3.27) 
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where  Lsin  is  the  LUT  sine  output  and  01000  is  the  5 -bit  phase  representation  for  90  de¬ 
grees  as  given  in  Table  2. 

3.  DIS  Amplitude  Modulation 

The  signals  in  the  I  and  Q  channels  are  amplitude  modulated  by  multiplying  each 
binary  signal  by  the  gain  modulation  coefficients  Similar  to  the  phase  modulation, 

the  gain  coefficients  are  fed  separately  for  each  range  bin  and  updated  for  every  new 
ISAR  pulse  (each  PRI).  The  use  of  a  multiplier  at  this  stage  is  avoided  due  to  fast  speed 
requirements.  Instead,  the  gain  modulation  is  implemented  by  shifting  bit  locations.  A 
gain  shifter  module  shifts  the  binary  digits  by  g(r,n)  places  (from  0  to  10)  using  four  ar¬ 
rays  of  2-bit  to  1-bit  multiplexers.  An  input/output  diagram  of  the  gain  shifter  is  shown 
in  Figure  16.  Details  of  the  multiplexer  implementation  can  be  found  in  Fonts  et  al.  [9]. 


la - 

ja - 


Gain  Shifter 


aiE  ais 


ai3  ai! 


■n - ^ - ^ - ^ - ^ - ^ - ^ - ^ - ^ - ^ - ^ - ^ - T" 


Figure  16  Gain  Shifter  Diagram  (From:  [14]). 


The  left  side  of  the  gain  shifter  diagram  in  Figure  16  shows  four  control  inputs 
and  the  negated  equivalents.  These  inputs  represent  a  set  of  16  control  codes  to  select 
gains  according  to  the  list  included  in  Table  4.  The  missing  control  codes  of  4,  8,  9,  10 
and  12  are  discarded  because  they  produce  a  shift  also  produced  by  a  lower  control  code. 
For  instance,  the  control  codes  0011  (3)  and  0100  (4),  when  applied  to  the  multiplexer 
array,  both  produce  a  shift  of  three  binary  digits. 

The  outputs  of  the  I  and  Q  channels  after  gain  modulation  are  amplitudes  from 
-1028  to  1028-2“’.  They  are  represented  in  an  18-bit  two’s  complement  format  with 
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one  sign  bit,  ten  integer  bits  from  2*^  to  2®  and  seven  decimal  bits  from  2“' to  2“’.  The 
format  described  is  then  truncated  to  the  five  least  significant  bits  (LSB).  Thus,  each 
channel  is  composed  of  a  13 -bit  format  with  a  sign  bit,  ten  integer  bits  and  two  decimal 
bits,  as  listed  in  Table  5. 


Control 

Code 

Size  of  Shift 
g(r,n) 

Multiplication 
Factor  2^^''’"^ 

0 

0 

1 

1 

1 

2 

2 

2 

4 

3 

3 

8 

5 

4 

16 

6 

5 

32 

7 

6 

64 

11 

7 

128 

13 

8 

256 

14 

9 

512 

15 

10 

1024 

Table  4  Gain  Shifter  Control  Codes. 


Bit  Position 

13 

12 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

Value 

Sign 

2*^ 

2^ 

2’ 

2^ 

2' 

2^ 

2^ 

2^ 

2^ 

2^ 

2'^ 

2'^ 

Table  5  Bit  Format  of  the  Gain  Modulation  Output  after  Truncation  of  5  LSBs. 

The  use  of  bit  shifting  from  0  to  10  binary  places,  instead  of  a  multiplier,  restricts 
the  gains  to  a  binary  exponential  distribution  from  1  to  1028.  This  gain  distribution  was 
the  main  design  consideration  in  choosing  the  gain  coefficient  quantization  scheme  given 
in  Table  1. 

4.  DIS  Cumulative  Sum 

The  complex  outputs  for  all  range  bins  are  added  to  form  a  consolidated  output 
(return  signal)  where  the  scattering  effects  of  the  different  range  locations  of  the  target 
are  superimposed.  The  cumulative  sum  of  the  modulated  outputs  start  from  the  output  of 
the  last  range  bin  being  added  with  the  adjacent  range-bin  output.  On  the  next  clock 
pulse,  the  result  of  this  sum  is  added  with  the  next  adjacent  range  bin  until  all  range-bin 
outputs  are  added.  The  implementation  of  this  cumulative  sum  is  shown  in  Figure  17. 
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Range  Bin  0 


Range  Bin  Nj.-3 


Range  Bin  Nj.-2  Range  Bin  N^.-! 


IQ  IQ  I  Q  I  Q 


The  hardware  implementation  of  the  cumulative  sum  uses  a  16-bit  adder  for  each 
channel  (I  and  Q)  in  each  range  bin.  A  diagram  of  this  adder  with  the  input  and  output 
connections  is  shown  in  Figure  18.  The  13-bit  truncated  signal  (gain  modulation  output) 
with  the  format  described  in  Table  5  is  fed  into  one  of  the  adder  inputs.  To  complete  the 
three  missing  bit  locations,  the  sign  bit  at  location  13  is  also  fed  at  locations  14,  15  and 
16.  The  replication  of  the  sign  bit  in  the  three  most  significant  bits  is  equivalent  to  a 
‘zero’  padding  for  a  positive  number  and  a  ‘one’  padding  for  a  negative  number.  Both 
cases  do  not  alter  the  value  of  the  data  in  two’s  complement  form.  The  second  adder  in¬ 
put  is  fed  with  the  16-bit  output  from  the  cumulative  sum  of  previous  range  bins. 


Figure  18  16-bit  Adder  Diagram. 
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The  DIS  I  and  Q  output  channels  have  the  16-bit  format  shown  in  Table  6.  The 
DIS  sends  the  outputs  channels  to  the  digital-to-analog  converters  (DAC).  The  analog 
signals  are  then  up-converted  in  frequency  and  transmitted  to  the  victim  radar. 


Bit  Position 

16 

15 

14 

13 

12 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

Value 

Sign 

212 

2^^ 

2io 

2*^ 

2^ 

2^ 

2^ 

2= 

2^ 

2^ 

2^ 

2^ 

2^ 

2'^ 

2'^ 

Table  6  Bit  Format  of  the  DIS  Output  (per  Channel). 


5.  DIS  Control  Signals 

The  DIS  is  controlled  by  an  off-chip  microprocessor  that  generates  appropriate 
control  signals.  Description  of  these  DIS  control  signals  follow: 

a.  Use  Range  Bin  (URB) 

This  signal  is  used  to  enable  or  disable  a  range  bin  as  appropriate,  accord¬ 
ing  to  the  target  extent  or  the  existence  of  low  radar  cross-section  segments.  If  the  URB 
signal  is  a  logic  one,  the  range  bin  is  included  in  the  cumulative  sum  chain  of  Figure  17. 
When  a  logic  zero  is  applied  in  the  URB  signal,  the  range  bin  is  closed.  This  means  that 
it  feeds  a  zero  to  the  output  adder.  Thus,  a  closed  range  bin  propagates  the  cumulative 
sum  to  the  next  range  bin  without  incrementing  it. 

b.  Phase  Sample  Valid  (PSV) 

This  signal  indicates  when  data  coming  from  the  DRFM  is  valid.  This 
signal  is  only  enabled  during  the  pulse  duration  and  disabled  during  the  interpulse  period. 

c.  Output  Data  Valid  Out  (ODVout) 

This  signal  is  used  to  cascade  additional  DIS  chips  in  order  to  synthesize 
very  large  targets.  The  signal  indicates  that  the  output  of  the  last  range-bin  processor  in 
the  summation  chain  of  the  previous  DIS  chip  is  valid. 

d.  Output  Data  Valid  In  (ODViJ 

This  signal  is  used  in  conjunction  with  the  ODVout  signal  to  cascade  DIS 
chips.  The  ODVout  signal  from  a  cascaded  chip  is  fed  to  the  next  chip  as  the  ODVin  sig¬ 
nal.  For  the  case  of  the  first  chip,  the  ODVin  signal  is  always  enabled. 
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e.  Program  Range  Bin  (PRB) 

This  signal  causes  the  phase  and  gain  modulation  coefficients  to  be  loaded 
into  the  range  bins.  This  signal  clocks  a  preload  buffer  register  to  update  these  values 
every  PRI. 

f.  Use  New  Programming  (UNP) 

This  signal  causes  the  phase  and  gain  rotation  registers  to  load  the  respec¬ 
tive  modulation  coefficients  from  the  pre-load  registers  every  PRI. 

6.  DIS  Clock 

The  DIS  chip  is  a  synchronous  digital  integrated  circuit.  An  off-chip  clock  signal 
is  required  to  control  the  flow  of  data  through  the  internal  logic  and  to  synchronize  the 
DIS  chip  with  the  DRFM  providing  the  input  data  samples,  and  the  DAC  receiving  the 
output  of  the  DIS. 

G.  CHAPTER  SUMMARY 

The  DIS  concept  was  described  mathematically  from  the  generation  of  the  input 
signals,  including  the  signal  processing  at  the  complex  range-bin  modulators,  and  finish¬ 
ing  with  the  generation  of  the  output-jamming  signal.  The  process  that  takes  place  at  the 
ISAR  processor  (range  and  azimuth  compression)  generates  a  range  Doppler  image  of  the 
intended  false  target.  A  brief  description  of  the  DIS  hardware  was  included  to  understand 
the  pipelined  process  that  occurs  at  each  range-bin  modulator.  In  the  next  chapter,  this 
process  is  represented  with  a  software  model  to  simulate  the  DIS  input  signals,  the  signal 
processing,  and  the  follow-on  ISAR  processing  that  results  in  the  target  image  generation. 
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IV.  DIS  SIMULATION  ALGORITHM 


This  chapter  describes  a  software  model  created  to  simulate  the  performance  of 
the  DIS  system  design  against  chirp-pulse  ISARs.  The  software  model  is  implemented 
with  MATLAB  version  6. 1 .  Rather  than  analyzing  the  software  code  in  a  per  line  basis, 
only  the  major  signal  processing  steps  will  be  addressed,  presenting  the  code  lines  in¬ 
volved  for  each  case.  This  chapter  must  be  read  in  conjunction  with  the  complete  soft¬ 
ware  codes  included  in  Appendix  A  in  order  to  understand  the  origin  of  variables  that  are 
defined  before  the  presented  code  lines. 

The  MATLAB  code  is  composed  of  three  main  files  and  a  number  of  auxiliary 
files.  The  three  main  files  must  be  run  in  sequence.  First,  the  extract_v5.m  file  generates 
the  modulation  coefficients  for  a  false  target.  Next,  the  mathost_v5.m  file  simulates  the 
DRFM  samples  used  by  the  DIS  and  the  reference  signal  used  at  the  ISAR  receiver.  Fi¬ 
nally,  the  simhwchk_v5.m  file  simulates  the  DIS  signal  processing,  generates  the  deceiv¬ 
ing  signal  and  simulates  the  ISAR  compression  process  to  produce  a  target  range-Doppler 
image. 

A.  GENERATION  OF  FALSE  TARGET  MODULATION  COEFFICIENTS 

1.  Loading  of  a  Target  Doppler  Template 

The  extract_v5.m  file  loads  the  range-Doppler  template  of  a  target  as  presented  in 
Figure  10.  The  template  contains  the  main  target  scatterers  arranged  according  to  their 
range  location  r  and  Doppler  shift  value  f(r,d).  The  r  coordinate  is  the  range  position 
within  the  ship  structure  in  terms  of  range  cell  units,  usually  in  the  horizontal  dimension. 
The  d  coordinate  is  the  cross-range  or  azimuth  location  in  terms  of  the  Doppler  resolution 
units  or  Doppler  filters,  usually  in  the  vertical  dimension. 

The  size  of  the  DIS  hardware  array  (number  of  range-bin  processors)  is  deter¬ 
mined  by  the  target  extent  in  range  cell  units.  The  current  software  model  has  been  set  to 
simulate  DIS  array  sizes  for  2,  32,  128  and  512  range-bin  processors.  This  is  achieved  by 
loading  false  target  templates  from  the  data  files  targetOl.mat,  targets 2. mat,  tar- 
getl28.mat  and  target512.mat  that  are  created  from  the  files  doppler_02_target.m,  dop- 
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plerjOS.m,  doppler_128.m,  doppler_32.m  and  doppler_512.m.  The  creation  of  new  files 
with  target  templates  allows  the  simulation  of  additional  DIS  array  sizes. 

The  extract _y 5. m  file  starts  with  the  definition  of  the  numtaps  variable  as  the 
number  of  range-bin  processors  in  the  DIS  hardware  array  or  target  size.  The  target 
Doppler  template  file  is  chosen  and  loaded  depending  on  the  value  of  the  numtaps  vari¬ 
able.  The  Doppler  frequency  values  are  saved  into  a  matrix  labeled  as  freq.  A  plot  of 
this  matrix  is  displayed  by  the  routine  in  order  to  corroborate  the  target  template  that  will 
be  synthesized.  Figure  19  shows  the  display  of  the  freq  matrix,  indicating  the  range 
Doppler  arrangement  of  the  scatterers  for  a  32-range  cell  target. 


Figure  19  Scatterer  Distribution  for  a  32-Range-bin  Target  Using  128  Doppler  Cells. 

2.  Generation  of  a  Complex  Scatterer  Expression 

The  Doppler  frequency  values  from  the  freq  matrix  can  be  used  in  conjunction 
with  time  and  magnitude  values  to  produce  a  complex  scatterer  expression  in  polar  form. 
Hence,  an  array  of  target  scatterers  T(r,d)  for  each  pulse  is  obtained  using  a  common 
magnitude  matrix  A(r,d),  equivalent  to  the  RCS  of  the  scatterers,  and  phase  matrices  ob¬ 
tained  from  the  product  of  2k,  the  Doppler  shift  /  contained  in  the  freq  matrix  and  pulse 
dependant  time  values: 

T{r,d,t)  =  A{r,d)e-^^’'^^'''^^\  (4.1) 
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The  time  used  in  the  phase  produet  of  (4.1)  must  be  the  leading  edge  time  of  the  respec¬ 
tive  pulse,  thus  it  can  be  expressed  as  the  product  n  PRI ,  where  n  is  the  pulse  index. 

Replacing  this  product  in  the  time  expression  of  (4.1),  the  array  of  scatterers  can  be  ex¬ 
pressed  as  a  function  of  the  pulse  index  n  resulting  in  the  T(r,d,n)  three-dimensional  scat- 
terer  array  previously  defined  in  Equation  (3.1). 

The  ISAR  total  integration  time  is  obtained  from  the  product  of  the  PRI  and  the 
number  of  pulses  integrated  at  the  radar  receiver  Np,  equivalent  to  the  number  of  Doppler 
resolution  cells.  The  value  for  Np  is  defined  by  the  dp _pts  variable.  The  integration  time 
is  calculated  as  PRI  and  saved  under  the  variable  name  of  tstop  at  the  beginning  of 

the  extract_v5.m  file.  Knowing  the  integration  time  period  and  the  number  of  pulses  in¬ 
volved,  it  is  possible  to  define  a  vector  with  the  leading  edge  time  for  every  pulse  as: 

t  =  nPRl  fov  n  =  \..N^.  (4.2) 

This  time  vector  is  defined  in  the  simulation  model  as  ti  in  the  following  code  lines: 

tstop  =  pri*dp_pts;  %  Pulse  train  duration  (integration  time) 

ti  =  linspace(0, tstop, dp_pts):  %  Vector  with  pulse  leading  edge  times 

The  T(r,d,t)  expression  in  (4.1)  can  be  defined  as  T(r,d,n)  using  the  ti  vector.  The 
T(r,d,n)  expression  contains  a  complex  scattering  matrix  for  each  of  the  Np  pulses  and  is 
saved  as  target.  This  three-dimensional  expression  can  be  reduced  to  a  two-dimensional 
matrix  T\r,n) ,  by  adding  the  scatterers  coherently  in  the  cross-range  or  azimuth  dimen¬ 
sion.  The  two-dimensional  result  T\r,n)  is  saved  as  the  targetSum  matrix  and  represents 
a  collection  of  range  profiles  (rows)  for  each  pulse  n.  This  matrix  is  equivalent  to  the  ma¬ 
trix  of  range  profiles  formed  after  the  ISAR  range  compression  output.  Thus,  the  process 
of  generating  T\r,n)is  inverse  to  the  ISAR  azimuth  compression.  The  calculation  of 
T(r,d,n)  {target)  and  T\r,n)  {targetSum)  is  performed  using  a  three  loop  structure  as 

shown  in  the  following  code  lines: 

%  GENERATION  OF  COMPLEX  SCATTERER  EXPRESSIONS 
A  =  ones(numtaps,dp_pts):  %  Scatterer  Amplitudes  (uniform) 

targetSum  =  zeros(numtaps,dp_pts); 

for  i1  =  1  :numtaps,  %  At  each  range  bin 

target  =  zeros(dp_pts,length(ti)): 

for  i2  =  1  :numfreq(i1 ),  %  At  each  scatterer 

for  i3  =  1  :length(ti),  %  At  each  integrated  pulse 

target(i2,i3)  =  A(i1,i2)  *  exp(-j*2*pi*freq(i1,i2)*ti(1,i3)): 
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%  Creates  complex  scatterer 


end 
end 

for  14  =  1  :numfreq(i1 )  %  Combined  signal  p/rbin  p/pulse 

targetSum(i1,:)  =  targetSum(i1,:)  +  target(i4,:): 
end 
end 

where  il,  i2  and  i3  are  loop  variables  and  numfreq(il)  is  the  number  of  scatterers  co¬ 
existing  within  the  same  range  bin  il. 

3.  Magnitude  Coefficients 

The  magnitudes  of  the  array  of  range  profiles  T\r,n)  are  normalized  and  quan¬ 
tized  into  1 1  levels.  Each  level  is  assigned  to  a  gain  coefficient.  The  quantization  levels 
are  exponentially  distributed  between  threshold  values  given  by: 

(4.3) 

where  5  goes  from  1  to  10,  and  fq  is  an  arbitrary  factor  between  1  and  2.  The  resulting 
quantization  scheme  is  shown  in  Table  7  for  an  fq  value  of  1.5.  The  threshold  values 
slightly  differ  from  those  in  the  quantization  example  in  Table  1,  where  they  were  gener¬ 
ated  with  an  fq  value  of  1.6.  The  choice  of  the  fq  value  has  a  direct  effect  on  the  number 
of  times  that  the  gain  coefficient  has  the  maximum  value  of  1024.  Therefore,  this  factor 
can  be  used  to  control  the  number  of  overflow  occurrences  in  the  DIS  output  adder.  For 
example,  if  an  of  1.8  is  used,  the  gain  coefficient  of  1024  will  be  only  assigned  to  the 
scatterer  locations  where  the  normalized  magnitude  of  T\r,n)  has  a  value  from  0.9  to 
1.0.  Thus,  the  probability  of  obtaining  an  overflow  is  greatly  reduced. 


Magnitude  Range 

Number  of  Shifts  g(r,n) 

Effective  Gain 

0.75  to  1.0 

10 

1024 

0.375  to  0.75 

9 

512 

0.1875  to  0.375 

8 

256 

0.0938  to  0.1875 

7 

128 

0.0469  to  0.0938 

6 

64 

0.0234  to  0.0469 

5 

32 

0.01 17  to  0.0234 

4 

16 

0.0059  to  0.0117 

3 

8 

0.0029  to  0.0059 

2 

4 

0.0015  to  0.0029 

1 

2 

Oto  0.0015 

0 

1 

Table  7  Quantization  Scheme  for  DIS  Magnitude  Coefficients. 
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The  gain  quantization  scheme  described  in  Table  7  is  implemented  in  the  simula¬ 


tion  program  using  the  following  code  lines: 


%  GENERATION  OF  GAIN  MODULATION  COEFFICIENTS 


amp  =  abs(targetSum): 
gain  =  ones(numtaps,dp_pts); 
fq  =  1.5; 

for  i1  =  Tnumtaps, 
for  12  =  1  :dp_pts, 

normAmp  =  amp(l1,l2)/max(max(amp)); 
for  s=1;10, 

if  normAmp  >  (2'^-(1 1-s))*fq, 
gain(i1,i2)  =  2^s: 
end 
end 
end 
end 


%  Magnitude  of  combined  signal. 
%  Default  gain  coefficient  is  1 . 

%  Arbitrary  threshold  factor 


%  Magnitudes  are  normalized. 

%  Thresh;  0.001 5...., 0.375, 0.75 

%Gain:  2,4,8,16, . 512,1024 

%  Assigns  1 1  power-of-two  gains 


Figure  20a  shows  the  resulting  amplitude  distribution  of  gain  coefficients  for  a  test  target 
of  32  range  cells  with  an  integration  of  128  pulses.  To  prove  that  the  exponential  as¬ 
signment  of  gain  coefficients  is  evenly  generated,  a  logarithmic  plot  of  the  gain  coeffi¬ 
cient  distribution  is  shown  in  Figure  20b. 


Figure  20  Quantization  Scheme  for  Gain  Coefficients  of  a  32-Range-bin  Target. 


The  plot  of  gain  coefficients  can  also  be  rearranged  with  respect  to  their  range  cell/pulse 
location.  The  resulting  plot  is  shown  in  Figure  21.  This  plot  shows  occurrences  of  the 
highest  gain  coefficient  of  1024  at  the  combination  of  pulses  1,  2,  63,  64,  65  and  66,  and 
range  cells  9  and  25  corresponding  to  the  target  mast  locations  in  Figure  19.  These  peak 
gain  locations  are  potential  sources  of  overflow  in  the  output  adder  bit  format.  In  addi- 


47 


tion,  the  figure  shows  that  the  distribution  of  gain  coefficients  is  symmetric  between 
pulses  1  through  64  and  65  through  128. 


5  10  15  20  25  30 


Down  Range  Cells. 

Figure  21  Gain  Coefficient  Distribution  on  a  Range  Cell-Pulse  Index  Map. 

4.  Phase  Coefficients 

The  phase  of  the  array  of  range  profiles  is  extracted  by  applying  the  angleQ  func¬ 
tion  to  the  targetSum  array.  The  resulting  phase  values  are  then  converted  into  a  scale 
from  0  to  31.99.  For  each  scatterer  phase,  a  phase  increment  is  calculated  with  respect  to 
the  scatterer  of  the  previous  Doppler  location.  The  first  Doppler  location  is  assigned  a 
phase  value  of  zero.  At  each  scatterer  location,  the  cumulative  sum  of  phase  increments 
results  in  a  relative  phase  value  with  respect  to  the  first  Doppler  location.  These  relative 
phase  values  are  stored  in  a  matrix  labeled  as  phaseinc.  The  DIS  phase  modulation  coef¬ 
ficients  are  obtained  from  the  quantization  of  the  values  of  the  phaseinc  matrix  into  32 
integer  phase  values  from  0  to  31.  These  coefficients  are  stored  in  a  matrix  labeled 
phasecoeff.  The  calculation  of  the  phasecoeff  phase  coefficient  matrix  is  implemented  is 
the  following  code  lines: 

%  GENERATION  OF  PHASE  MODULATION  COEFFICIENTS 
vphasel  =  angle(targetSum):  %  Phase  of  combined  signal 

vphase2  =  vphasel  *2^nbitsdop/(2*pi):  %  Adjusts  values  from  0  to  31 .99 

phaseinc  =  zeros(numtaps,dp_pts): 
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%  Turn  phases  into  phase  increments 


for  i6  =  1:numtaps, 
for  i7  =  1:dp_pts, 
if  i7==1, 

phaseinc(i6,i7)  =  phaseinc(i6,i7)+vphase2(i6,i7): 
else  %  Finds  phase  increments  bet/pulses 

phaseinc(i6,i7)  =  phaseinc(i6,i7-1  )+vphase2(i6,i7-1  )-vphase2(i6,i7): 
end 
end 
end 

for  i8  =  1  :dp_pts, 

for  i9  =  1  :numtaps,  %  5-bit  quantization  in  32  levels. 

phasecoeff(i9,i8)  =  flx(rem(phaseinc(l9,l8)+32,32)): 
end 
end 

Similar  to  Figure  21,  the  phase  modulation  coefficients  can  be  plotted  with  re¬ 
spect  to  their  range  cell/pulse  location.  The  resulting  plot  is  shown  in  Figure  22.  This 
distribution  shows  that  the  phase  variation  becomes  faster  for  higher  range  bins.  In  addi¬ 
tion,  there  is  a  pseudo  symmetry  between  the  phase  coefficients  from  pulses  1  through  64 
and  pulses  65  through  128. 


5.  Storage  of  Coefficients 

The  gain  and  phase  coefficient  values  are  stored  in  a  text  file  labeled  as  para- 
MULTIq4Ship.txt  after  a  header  of  three  identification  values.  The  first  three  stored 
quantities  are  the  number  of  intrapulse  samples  (500  for  the  default  values  of  0.5  us 


pulsewidth  and  1  GHz  sampling  frequency),  the  number  of  Doppler  cells  (default  is  128) 
and  the  target  extent  in  range  cells  (default  is  32).  For  the  default  case  of  32  range  bins 


and  128  Doppler  cells,  the  total  number  of  gain  coefficients  is  4096.  For  the  case  of  these 
default  settings,  the  gain  modulation  coefficients  are  stored  from  the  4*  location  to  the 
4099*  location.  The  phase  modulation  coefficients  are  stored  from  the  4100*  location  to 
the  8195*  location.  The  storage  of  the  gain  and  phase  modulation  coefficients  in  the  pa- 


raMULTIq4Ship.txt  file  is  implemented  throughout  the  following  code  lines: 


%  STORAGE  OF  MODULATION  COEFFICIENTS 


f4  =  fopen('paraMULTIq4Shlp.txt','w'): 
fprlntf(f4,'%d\r\n',rag_pts): 
fprlntf(f4,'%d\r\n',dp_pts): 
fprintf(f4,'%d\r\n',numtaps): 
for  aa=1:dp_pts, 
for  bb=1  :numtaps, 

fprintf(f4,'%d\r\n',gain(bb,aa)); 

end 

end 

for  aa  =  1  :dp_pts, 
for  bb  =  fnumtaps, 


%  Open  file  to  store  coefficients 
%  Stores  #  of  chirp  pulse  samples 
%  Stores  #  of  Doppler  bins 
%  Stores  #  of  tap  lines  (target  extent) 


%  Stores  gain  coefficients 
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fprintf(f4,'%d\r\n',phasecoeff(bb,aa));  %  Stores  phase  coefficients 
end 
end 

fclose(f4): 


Figure  22  Phase  Coefficient  Distribution  on  a  Range  Cell-Pulse  Index  Map. 


B.  SIMULATION  OF  DIS  INPUT  AND  ISAR  PULSE  COMPRESSION 

REFERENCE  SIGNALS 

I.  DIS  Input  Signal 

The  DIS  input  signal  is  a  digitized  version  of  the  ISAR  chirp  pulse  that  is  used  as 
the  main  input  for  the  DIS.  Each  pulse  is  fed  into  the  DIS  as  a  batch  of  DRFM  phase 
samples.  Due  to  the  squared  time  term  in  the  exponent  of  the  chirp  phase  in  (3.11),  as 
time  increases  the  phase  values  become  very  large  rotated  values.  If  the  Doppler  effect  is 
ignored  as  stated  in  that  equation,  the  chirp  phase  as  a  function  of  time  can  be  represented 
as: 

=  Ink  1^/2  j  =  (4.4) 

where  k  is  the  compression  factor  equivalent  to  the  ratio  of  the  bandwidth  A  over  the 
pulsewidth  r .  The  resulting  phase-value  samples  for  each  pulse  are  stored  in  a  vector 


50 


labeled  as  oldphase.  These  rotated  values  converted  to  equivalent  phase  values  within 
the  0  to  360  degrees  interval  and  then  quantized  into  32  integer  phase  levels  from  0  to  31. 
The  quantized  values  correspond  to  the  5 -bit  representation  of  the  phase  sampling  DRFM 
output.  The  quantized  phase  samples  for  all  the  integrated  pulses  are  stored  in  a  matrix 
labeled  as  pri_rg_phaseq.  These  DIS  phase  inputs  are  arranged  by  rows  for  each  differ¬ 
ent  radar  pulse.  A  baseband  digital  representation  of  a  chirp  pulse  can  be  obtained  by 
plotting  one  row  of  the  pri_rg_phaseq  matrix.  The  resulting  chirp  waveform  is  plotted  in 
Figure  23. 


The  reference  signal  that  is  used  at  the  receiver  to  perform  the  range  compression 
of  each  returned  pulse  can  be  obtained  for  the  oldphase  vectors.  These  vectors  are  stored 
as  columns  of  a  matrix  named  cref.  The  process  used  to  generate  the  pri_rgj}haseq  DIS 
input  matrix  and  the  cref  range  compression  reference  matrix  is  implemented  with  the 


following  code  lines: 


%  GENERATION  OF  DRFM  PHASE  SAMPLES 
pri_rg_phaseq  =  zeros(dp_pts,rag_pts): 
p  =  2*pi/(2''nbitsdop): 
pri  =  1/prf: 
for  idx1  =  1  :dp_pts, 
t1  =  to  +  idx1*pri: 
oldphase  =  2*pi*(k*t1.*t1)/2: 
oldphase  =  renn(oldphase,2*pi): 
prLi'g_phaseq(iclx1 ,:)  =fix(oldphase/p): 
crefc(:,idx1)  =  exp(j*oldphase.'); 
end 


%  Initialize  DRFM  phase  matrix 
%  Phase  quantization  factor 

%  Loop  per  pulse 
%  PRI  time  shift  for  all  samples 
%  Chirp  pulse  phase  samples 
%  limit  phase  from  0  to  2pi 
%  Phase  quantized  in  2"  levels 
%  ISAR  reference 


where  nbitsdop  has  a  value  of  5  to  obtain  2^  phase  levels. 


2.  Storing  of  DRFM  Samples  and  Reference  Signal 

The  phase  samples  contained  in  the  pri_rg_phaseq  matrix  are  stored  row  by  row 
in  the  rawint.txt  file.  Each  row  of  this  text  file  corresponds  to  the  phase  samples  of  one 


radar  pulse.  The  last  element  of  the  matrix  is  stored  followed  by  a  carrier  return  character 


to  skip  to  the  next  line  of  the  data  file.  The  code  lines  involved  are  listed  below: 


%  STORAGE  OF  DRFM  PHASE  MATRIX  : 
f2  =  fopen('rawint.txt','w'): 
for  i  =  1  :dp_pts, 

int_raw  =  pri_rg_phaseq(i,1:rag_pts-1): 
fprintf(f2,'%d,',int_raw): 
int_raw  =  fpri_rg_phaseq(i,rag_pts): 
fprintf(f2,'%d\r\n',int_raw); 
end; 

fclose(f2): 


%  Opens  file  for  phase  storage 
%  Loop  for  all  pulses. 

%  Stores  pulse  phase  samples 

%  Store  last  sample  of  the  pulse 
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Figure  23  Baseband  Plot  of  the  Intercepted  ISAR  Pulse  as  Generated  by  a  DRFM. 


The  coherent  reference  for  each  pulse  needs  to  be  stored  in  order  to  be  available 
for  the  simhwchk_v5.m.  This  file  performs,  among  other  things,  the  cross-correlation  to 
simulate  the  ISAR  range  compression.  Since  the  cross-correlation  uses  the  complex  con¬ 
jugate  of  the  Fourier  transform  of  the  reference  (Figure  6),  the  complex  conjugation  and 
the  transform  process  is  performed  before  storing  the  signals.  To  approximate  the  ideal 
infinite  signal  condition,  the  FFT  is  performed  using  2rg _pts-\  points,  where  rg_pts  is 
set  arbitrarily  as  the  maximum  between  the  number  of  samples  within  a  radar  pulse  and 
ten  times  the  target  extent  in  range  bins.  The  FFT  of  the  reference  is  then  complex  con¬ 
jugated  and  stored  in  a  matrix  labeled  as  crefc,  where  each  column  corresponds  to  the 
conjugated  FFT  vector  of  a  pulse  reference. 

%  GENERATION  OF  ISAR  REF.  FOR  RANGE  COMPRESSION  -  FAST  CORRELATION 
cref  =  conj(fft(crefc,2*rg_pts-1 ));  %  FFT  of  ISAR  reference 

save  pc_ref  cref  t1  %  Save  FFT  of  ISAR  reference 
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3.  Execution  of  Follow-on  Routines 

At  the  end  of  the  mathost_v5.m  file,  there  are  instructions  to  execute  the 
simhwchk_v5.m  and  the  Image _Analy sis. m.  The  first  routine  emulates  the  DIS  and  ISAR 
compression.  The  second  routine  is  an  image  quality  evaluation  performed  by  compari¬ 
son  between  the  current  output  image  and  an  ideal  infinite  resolution  quantization  image. 
This  analysis  is  only  executed  for  cases  where  the  target  has  a  size  of  32  range  bins  and 
the  number  of  Doppler  resolution  cells  is  128,  as  the  ideal  image  is  made  available  only 
for  that  case.  The  code  lines  involved  in  this  procedure  are  listed  below: 

%  DIS  SIMULATION  and  IMAGE  GENERATION  AN  ANALYSIS 
sinnhwchk_v5:  %  DIS  &  ISAR  processing 

if  numtaps==32  &  nDopplerCell==128  &  prf==200, 

lnnage_Analysis:  %  Comparison  with  ideai  image 

end; 

C.  SIMULATION  OF  DIS  AND  ISAR  SIGNAL  PROCESSING 

The  DIS  and  ISAR  simulations  are  implemented  by  the  simhwchk_v5.m  file.  The 
input  parameters  for  the  execution  of  this  file  are  the  DRFM  phase  samples  from  the 
ISAR  intercepted  pulses  and  the  DIS  modulation  coefficients.  The  DIS  processes  the 
ISAR  pulses  sequentially.  Within  each  pulse,  the  phase  samples  are  also  processed  in  the 
same  manner.  The  simulation  model  emulates  the  DIS  using  a  loop  structure  that  is  re¬ 
peated  for  each  received  pulse.  Within  that  loop,  an  additional  looped  process  is  repeated 
for  each  intrapulse  sample.  The  DIS  hardware  processes  each  intrapulse  sample  simulta¬ 
neously  (in  parallel)  at  all  range  bins.  Similarly,  the  simulation  model  emulates  this 
process  using  vectors  with  one  element  at  each  range  tap.  The  samples  at  the  range-bin 
outputs  are  properly  delayed  and  added  to  build  the  DIS  output  modulated  pulses.  The 
ISAR  simulation  is  modeled  with  a  range  compression  of  the  DIS  output  pulses,  followed 
by  an  azimuth  compression.  The  end  result  is  a  matrix  with  intensity  values  representing 
the  target  range  Doppler  image.  A  flow  diagram  summarizing  the  processes  that  are 
modeled  by  the  simhwchk_v5.m  file  is  shown  in  Figure  24. 

I.  Loading  of  Modulation  Coefficients 

The  modulation  coefficients  are  loaded  from  the  paraMULTIq4Ship.txt  file.  This 
file  also  contains  the  values  for  the  number  of  samples  in  the  ISAR  chirp  pulse,  the  num¬ 
ber  of  pulses  integrated  by  the  radar  and  the  target  extent  in  range  bins.  The  number  of 
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pulse  samples  is  saved  under  the  variable  name  nRangeCell  and  is  used  to  control  the  to¬ 
tal  number  of  times  the  inner  loop  is  set  for  each  pulse.  The  number  of  pulses  integrated, 
which  also  corresponds  to  the  number  of  Doppler  resolution  cells,  is  saved  under  the 
variable  name  nDopplercell  and  is  used  to  control  the  number  of  times  the  outer  loop  is 
used.  The  target  extent  is  saved  under  the  variable  name  numtaps  and  determines  the 
length  of  the  vectors  that  are  used  to  simulate  the  parallel  process  at  the  range  bins. 


Load  DRFM  pulse 
and  Modulation  N- 
Coefficients 


Process  Sanple  at 
all  Range  Bins 


Store  Output  Sample 


Load  Output  Sanples 


Figure  24  Summarized  Flow  Diagram  of  the  simhwchk_v5.m  File. 

After  the  first  three  values  are  loaded  from  the paraMULTIq4Ship.txt  file,  the  gain 
coefficients  can  be  found  starting  from  the  fourth  value  contained  in  the  file.  For  the  case 
of  a  32-bin  target  and  128  integrated  pulses,  the  gain  coefficients  can  be  extracted  by 
reading  the  values  from  location  4  through  4099.  These  coefficients  are  stored  in  a  4096- 
element  vector  named  gainin.  The  elements  of  this  vector  are  rearranged  into  a  128  by 
32  element  matrix  labeled  as  gain.  The  rows  of  this  matrix  contain  coefficients  belong¬ 
ing  to  the  same  pulses,  while  the  columns  represent  same  range  bins.  The  phase  modula¬ 
tion  coefficients  can  be  read  immediately  after  the  last  gain  coefficient,  from  location 
5000  through  8195.  These  coefficients  are  stored  in  a  4096-element  vector  named  phi. 
Similarly  to  the  gain  coefficient  case,  the  phi  phase  coefficient  vector  is  rearranged  into  a 
128  by  32  element  matrix  labeled  as  phasinc,  where  coefficients  from  the  same  pulse  can 
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be  found  in  the  same  rows.  The  loading  of  modulation  coefficients  is  implemented  with 


the  following  code  lines: 

%  LOADING  OF  MODULATION  COEFFICIENTS 

fid  =  fopen('paraMULTIq4Ship.txt',V): 

tmp  =  fscanf(fid,'%f); 

nRangeCell  =  tnnp(1); 

nDopplerCell  =  tmp(2): 

numtaps  =  tnnp(3): 

gainin  =  tmp(4:4+numtaps*nDopplerCell-1): 
gain  =  reshape(gainln, numtaps, nDopplerCell); 
gain  =  gain'; 

phi  =  tmp(4+numtaps*nDopplerCell:end); 
phasinc  =  reshape(phi, numtaps, nDopplerCell); 
fclose(fid); 


%  Opens  file  w/DIS  coefficients 
%  Loads  file  contents 
%  number  of  pulse  samples 
%  number  of  pulses  integrated 
%  number  of  DIS  range  taps 
%  Extracts  gain  coefficients 
%  Arranges  them  in  a  matrix 
%  Transpose 

%  Extracts  gain  coefficients 
%  Arranges  them  in  a  matrix 


2.  Loading  of  DRFM  Phase  Signal 

The  simulated  DRFM  phase  samples  are  stored  in  the  rawint.txt  file.  The  file  has 


the  values  arranged  in  rows  and  separated  by  commas,  with  each  row  representing  a  dif¬ 
ferent  pulse.  Once  the  phase  samples  are  loaded,  they  are  stored  in  a  matrix  labeled  as 


raw,  where  each  row  also  represents  a  different  pulse.  The  loading  of  the  DRFM  signal 


is  obtained  with  the  following  code  lines: 

%  LOADING  OF  DRFM  SAMPLES 
raw  =  zeros(nDopplerCell, nRangeCell): 
fid  =  fopen('rawint.txt','r'); 
for  r  =  TnDopplerCell, 
for  i2  =  1  :nRangeCell-1 , 
raw(r,i2)=fscanf(fid,'%f,1): 
comma  =  fscanf(fid,'%c',1): 
end 

raw(r,nRangeCell)  =  fscanf(fid,'%f,1): 
end 

fclose(fid); 


%  Initialize  DRFM  phase  matrix 
%  Opens  file  w/DRFM  samples 
%  Loop  per  pulse 
%  Loop  per  sample 
%  Loads  DRFM  phase  sample 
%  Loads  comma  character 

%  Last  pulse  phase  sample 


In  order  to  reset  the  range-bin  inputs  after  the  last  DRFM  sample  has  been  fed  to 
the  DIS,  the  DRFM  phase  samples  are  appended  with  zeros  at  the  end  of  each  pulse. 
This  is  accomplished  by  concatenating  a  matrix  with  zeros  to  the  right  of  the  raw  matrix. 
The  resulting  matrix  is  labeled  as  phaDRFM.  The  code  lines  used  for  the  concatenation 
of  zeros  are  shown  below: 

[row, col]  =  size(raw): 

phaDRFM  =  [raw,zeros(row,numtaps-1 )];  %  Trailing  zeros  to  clear  taps 
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3.  Phase  Modulation 

The  two  loops  are  implemented  to  iterate  the  DIS  process  for  each  pulse  and  for 
each  sample  within  a  pulse.  The  first  step  that  takes  place  within  the  loops  is  the  phase 
addition  of  DRFM  phase  samples  with  phase  coefficients.  The  phase  samples  are  fed  in 
parallel  to  all  range  bins.  Thus,  it  is  necessary  to  transform  each  phase  value  into  a  vector 
of  equal  value  with  the  size  equal  to  the  total  number  of  range  bins.  The  resulting  vector 
is  labeled  as  tapin.  The  phase  coefficients  are  different  for  each  range  bin  and  might  be 
changed  each  time  the  arriving  samples  correspond  to  a  new  radar  pulse.  To  produce  the 
phase  modulation,  the  tapin  vector  containing  repetitions  of  the  DRFM  sample  is  added 
with  a  row  of  the  phasinc  phase  coefficient  matrix.  The  resulting  vector  is  labeled  phAd- 
dout  and  contains  the  modulated  phase  for  all  range  bins.  These  phase  values  are  later 
used  to  access  the  LUT.  The  code  lines  involved  in  the  DIS  phase  modulation  are  shown 
below: 

%  DIS  PHASE  MODULATION 

tapln(1:numtaps)  =  phaDRFM(batchCnt,inPlsCnt);  %  Loads  DRFM  Phase  in 

phAddOut=  tapln(1  :nunntaps)'+phaslnc(:,batchCnt);  %  Phase  adder  output 

4.  Look-up  Table  (LUT) 

The  output  of  the  phase  modulation  represents  the  phase  value  of  a  unitary  magni¬ 
tude  complex  expression.  The  sine  and  cosine  LUTs  convert  this  complex  expression 
into  in-phase  and  quadrature  phase  amplitudes.  For  the  simulation  model,  8-bit  resolu¬ 
tion  LUTs  are  stored  in  vector  form  using  two  text  fdes  labeled  as  cosine8.txt  and 
sine8.txt.  Additional  text  files  were  also  created  to  represent  LUTs  with  lower  bit  resolu¬ 
tions.  The  contents  of  the  cosine8.txt  and  sine8.txX  LUT  text  files  are  presented  in  Table 
8.  These  LUT  files  are  loaded  using  the  following  code  lines. 

load  -ascii  cosine8.txt  %  Loads  LUT  tables 

load  -ascii  sine8.txt 

Once  the  LUT  files  are  loaded,  their  contents  are  available  as  two  32-element  vec¬ 
tors  labeled  respectively  as  cosine8  and  sine8.  Each  element  of  these  vectors  represents 
one  possible  LUT  output.  The  LUT  input  phase  values  are  integer  numbers  from  0  to  3 1 . 
the  output  I  and  Q  values  are  contained  in  vectors  with  indices  from  1  to  32.  Therefore, 
the  downloading  of  a  specific  LUT  output  can  be  implemented  with  the  LUT  vector  name 
with  an  index  value  equivalent  to  the  phAddout  modulated  phase  added  by  one  as  in 
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Equation  (3.15).  The  output  value  is  calculated  simultaneously  for  all  range  bins  and  is 
stored  in  a  column  vector  labeled  as  lutOut.  The  calculation  of  the  LUT  output  is  ob¬ 
tained  using  the  following  code  lines: 

tmp  =  renn(phAddOut+depthLUT,depthLUT)+1 ;  %  Phase  limit  bet  0-31 

lutOut  =  cosine8(tmp)  +  j*sine8(tmp):  %  LUT  tabie  output 


cosine8.txt 

Cosine  ROM 

sine8.txt 

Sine  ROM 

1 

0.9922 

01111111 

0 

00000000 

2 

0.9766 

01111101 

0.1953 

00011001 

3 

0.9141 

01110101 

0.3828 

00110001 

4 

0.8281 

01101010 

0.5547 

01000111 

5 

0.7031 

01011010 

0.7031 

01011010 

6 

0.5547 

01000111 

0.8281 

01101010 

7 

0.3828 

00110001 

0.9141 

01110101 

8 

0.1953 

00011001 

0.9766 

01111101 

9 

0 

00000000 

0.9922 

01111111 

10 

-0.1953 

11100111 

0.9766 

01111101 

11 

-0.3828 

11001111 

0.9141 

01110101 

12 

-0.5547 

10111001 

0.8281 

01101010 

13 

-0.7031 

10100110 

0.7031 

01011010 

14 

-0.8281 

10010110 

0.5547 

01000111 

15 

-0.9141 

10001011 

0.3828 

00110001 

16 

-0.9766 

10000011 

0.1953 

00011001 

17 

-0.9922 

10000001 

0 

00000000 

18 

-0.9766 

10000011 

-0.1953 

11100111 

19 

-0.9141 

10001011 

-0.3828 

11001111 

20 

-0.8281 

10010110 

-0.5547 

10111001 

21 

-0.7031 

10100110 

-0.7031 

10100110 

22 

-0.5547 

10111001 

-0.8281 

10010110 

23 

-0.3828 

11001111 

-0.9141 

10001011 

24 

-0.1953 

11100111 

-0.9766 

10000011 

25 

0 

00000000 

-0.9922 

10000001 

26 

0.1953 

00011001 

-0.9766 

10000011 

27 

0.3828 

00110001 

-0.9141 

10001011 

28 

0.5547 

01000111 

-0.8281 

10010110 

29 

0.7031 

01011010 

-0.7031 

10100110 

30 

0.8281 

01101010 

-0.5547 

10111001 

31 

0.9141 

01110101 

-0.3828 

11001111 

32 

0.9766 

01111101 

-0.1953 

11100111 

Table  8  Contents  of  the  8-Bit  Resolution  LUT  Files  and  LUT  ROM. 


5.  Gain  Modulation 

The  I  and  Q  LUT  outputs  for  all  range  bins  are  stored  in  the  lutOut  complex  vec¬ 
tor.  The  gain  modulation  is  performed  by  multiplying  these  I  and  Q  values  by  the  gain 
coefficient  for  the  corresponding  range-bin/pulse  combination.  For  the  case  of  the  n* 
pulse,  the  gain  modulation  is  implemented  with  the  element-by-element  multiplication  of 
the  lutOut  vector  with  the  row  from  the  gain  matrix.  The  modulated  amplitudes  of  the 
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I  and  Q  channels  are  stored  in  a  complex  vector  labeled  as  gainOut.  The  following  code 
line  performs  the  gain  modulation: 

gainOut  =  gain(batchCnt,:)'.*lutOut:  %  Gain  modulation 

Once  the  intrapulse  sample  has  been  modulated  by  the  DIS,  the  range-bin  outputs  must 
maintain  an  output  equal  to  zero.  This  shutting  down  of  range-bin  taps  is  accomplished 
in  the  hardware  design  using  the  control  signals  described  in  Chapter  III.  The  simulation 
model  implements  this  action  by  setting  the  value  of  the  gainOut  vector  to  zero  after  the 
last  pulse  sample  has  been  processed.  The  following  code  lines  perform  the  shutting 
down  of  range-bin  taps: 

if  inPIsCnt  >  nRangeCell,  %  Shuts  down  taps 

gain0ut(1  :(inPlsCnt-nRangeCell))  =  0;  %  after  pulse  passes 

end 

6.  Output  Adder  and  Overflow  Analysis 

The  chain  of  output  adders  described  in  Figure  17  perform  a  cumulative  sum  of 
the  I  and  Q  range-bin  outputs,  applying  one  sample  delay  between  them.  The  I  and  Q 
range-bin  output  samples  represented  by  the  gainOut  complex  vector  are  collected  as 
rows  in  a  matrix  labeled  as  tapOut.  Hence,  the  rows  of  this  matrix  represent  outputs  from 
the  same  DRFM  samples  and  the  columns  represent  outputs  from  the  same  range  bins. 
To  implement  one  sample  delay  on  an  output  element,  the  row  location  in  the  tapOut  ma¬ 
trix  must  be  increased  by  one.  Similarly,  all  elements  loaded  in  the  tapOut  matrix  are 
delayed  by  shifting  a  specific  amount  of  rows  depending  on  their  range-bin  number.  In 
other  words,  the  range-bin  output  gainOut  vectors  are  loaded  in  the  tapOut  matrix  as  di¬ 
agonals. 

The  amplitude  of  the  I  and  Q  range-bin  outputs  are  represented  by  the  bit  format 
shown  in  Table  5.  This  bit  format  limits  the  minimum  resolution  to  the  value  of  the  LSB. 
The  DIS  simulation  emulates  this  resolution  limit  by  truncating  the  I  and  Q  channel  am¬ 
plitude  defined  in  (3.19)  as: 


S{r,m,n^ 

LSB(5'(r,m,n)) 


LSB  (^5' 


(4.5) 


where  the  LSB  function  is  equal  to  the  value  of  the  least  significant  bit  in  the  binary  for¬ 
mat  of  the  argument.  The  simulation  defines  an  I  and  Q  bit  format  LSB  as  0.25,  storing  it 
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under  the  variable  name  of  o_res.  The  code  used  to  implement  the  delay  between  range- 

bin  outputs  and  the  truncation  for  the  bit  format  resolution  is  listed  below: 

%  OUTPUT  ADDER  DELAY  AND  MINIMUM  RESOLUTION 
for  idx3=1  inumtaps, 

delay  =  idx3-1 ;  %  Out.adder  register  delay  bet/taps 

tapOut(inPlsCnt+delay,idx3)=fix(gainOut(idx3)/o_res)*o_res: 
end  %  Applies  minimum  register  resolution 

The  sum  of  range-bin  outputs  can  be  performed  by  adding  the  row  elements  from 
the  tapOut  matrix.  This  complex  value  result  is  equivalent  to  the  I  and  Q  components  of 
one  sample  of  one  DIS  output  pulse.  As  described  in  Figure  17,  the  sum  between  delayed 
range-bin  outputs  is  implemented  serially.  The  output  of  each  sum  is  represented  by  the 
bit  format  of  Table  6.  Thus,  it  is  only  able  to  represent  amplitude  values  limited  by: 


(4.6) 

=  2  MSB(/ (m,«))-LSB(/ (m,n) 

l)  (4.7) 

where  the  MSB  function  is  equal  to  the  value  of  the  most  significant  bit 

mat  of  the  argument. 

in  the  binary  for- 

Whenever  an  I  or  Q  sum  output  is  forced  to  represent  a  value  outside  these  limits 

there  is  an  overflow  in  the  adder  format.  The  main  effect  of  the  format  overflow  is  that 

the  MSB  produces  a  carry  that  toggles  the  next  bit,  thus  switching  the  sign  of  the  two’s 

complement  bit  format.  The  resulting  erroneous  outputs  for  cases  where  the  lower  or  up¬ 
per  limits  are  surpassed  are  given  by: 

I  {m,n)  =  I  +  2  MSB(/ (m,n)) 

(4.8) 

l{m,n)  =  l{m,n\^  „  -2  MSB(/(OT,n)). 

V  /  V  /  Upper_overjlow  \  \  // 

(4.9) 

The  simulation  model  performs  the  cumulative  sum  storing  the  I  and  Q  adder  out¬ 
puts  in  a  complex  variable  named  suml.  After  each  sum  the  real  and  imaginary  values  of 
the  cumulative  sum  are  stored  as  lout  and  Qout,  and  checked  in  case  they  fall  outside  the 
overflow  limits  in  (4.6)  and  (4.7).  For  the  16-bit  format  implemented  in  the  DIS  hard¬ 
ware  and  shown  in  Table  6,  the  limits  are  -8196  and  +8195.75.  Whenever  the  real  and 
imaginary  components  lout  and  Qout  fall  outside  these  limits,  the  cumulative  sum  repre- 
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sented  by  suml  is  corrected  to  apply  the  erroneous  overflow  outputs  given  in  (4.8)  and 
(4.9).  In  addition,  a  counter  represented  by  the  variable  name  overflow  is  incremented 
every  time  an  overflow  occurs.  Once  the  cumulative  sum  has  been  completed,  the  result¬ 
ing  DIS  output  sample  is  stored  in  a  vector  labeled  as  partial  Japsum.  This  vector  stores 
the  DIS  output  samples  that  belong  to  the  same  pulse.  Thus,  it  stores  one  complete  out¬ 
put  pulse.  The  code  lines  that  implemented  the  range-bin  output  cumulative  sum  and 
overflow  simulation  for  each  sample  are  listed  below: 

%  OUTPUT  ADDER  CUMULATIVE  SUM  AND  OVERFLOW 
sunn1=0: 

for  tt=min(inPlsCnt,numtaps):-1  :max(1  JnPIsCnt-nRangeCell), 
suml  =  sum1  +  tapOut(inPlsCnt,tt); 
lout=real(sum1); 

Qout=imag(sum1); 
if  lout  >=  limit, 

sum1  =  lout-2*limit  +  fOout; 
overflow  =  overflow  +  1 ; 
elseif  lout  <  -limit, 

sum1  =  lout+2*limit  +  j*Qout: 
overflow  =  overflow  +  1 ; 
end 

if  Qout  >=  limit, 

sum1  =  lout  +  j*(Qout-2*limit): 
overflow  =  overflow  +  1 ; 
elseif  Qout  <  -limit, 

sum1  =  lout  +  j*(Qout+2*limit); 
overflow  =  overflow  +  1 ; 
end 

peak=max(max(abs(lout),abs(Qout)),peak):  %  End  of  overflow  simulation. 
partial_tapsum(inPlsCnt,1)  =  suml; 
end 


7.  DIS  Output 

After  the  last  sample  of  the  pulse  is  processed,  the  output  pulse  represented  by  the 
samples  contained  in  the  partial  Japsum  vector  is  stored  as  a  row  in  a  matrix  labeled  as 
finalAdderOut.  The  resulting  matrix  contains  the  complete  DIS  output  complex  signal 
with  each  row  representing  one  of  the  modulated  pulses.  Figure  25  displays  a  plot  of  one 
DIS  output  pulse  obtained  from  the  simulation  model.  The  figure  is  obtained  by  plotting 
the  first  row  of  the  finalAdderOut  matrix  (first  output  pulse). 

The  DIS  output  I  and  Q  channels  represented  by  the  finalAdderOut  complex  ma¬ 
trix  are  fed  to  I  and  Q  DACs.  These  converters  should  ideally  have  the  same  bit  resolu¬ 
tion  as  the  DIS  outputs  to  avoid  information  losses.  In  addition,  the  DAC  speed  should 


%  Overflow  simulation  using  +/-limit 

%  Positive  overflow  (in-phase) 

%  Adder  output  w/overflow 
%  Increases  counter 
%  Negative  overflow  (quad-phase) 
%  Adder  output  w/overflow 
%  Increases  counter 

%  Positive  overflow  (in-phase) 

%  Adder  output  w/overflow 
%  Increases  counter 
%  Negative  overflow  (quad-phase) 
%  Adder  output  w/overflow 
%  Increases  counter 
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be  high  enough  to  keep  up  with  the  required  clock  frequency.  For  cases  where  a  DAC 
has  less  bits  that  the  16-bit  DIS  output,  only  the  most  significant  DIS  output  bits  are  con¬ 
nected.  The  exclusion  of  output  bits  produces  a  truncation  that  is  dependent  on  the  value 
of  the  DIS  output  bit  position  that  is  connected  to  the  DAC  LSB,  such  that: 


r{m,n) 


l{m,n) 

LSB(D^C) 


LSB(D^C). 


(4.10) 


The  value  of  the  DIS  bit  connected  to  the  DAC  LSB  depends  on  the  number  of  DAC  bits 
Nb  and  is  given  by: 

hiB(DAC)  =  r-"‘.  ,,,,, 


Figure  25  Plot  of  the  DIS  Output  Modulated  Chirp 

The  DIS  model  simulates  the  effect  of  feeding  the  DIS  output  channels  through 
DACs  of  different  bit  resolutions  using  Equation  (4.10)  and  defining  the  LSB  of  the  DAC 
as  the  variable  name  dac_res.  The  following  code  lines  show  the  implementation  of  the 
DAC  loss: 


61 


%  DIS  OUTPUT  GENERATION  AND  DAO 

finalAdderOut(batchCnt,:)=partial_tapsum.':  %  DIS  output  mod.  chirps 

dac_res  =  0.25;  %  DAC  resolution  0.25  for  16bit 

finalAdderOut  =  fix(finalAdderOut/dac_res)*dac_res: 
figure 

plot(  real(  finalAdderOut(1 ,:)  )),axis  tight  %  Plots  1  DIS  output  chirp  pulse 

title('Plot  of  the  DIS  output  chirp  pulse'); 

8.  ISAR  Range  Compression 

The  ISAR  range  compression  is  modeled  in  a  pulse-by-pulse  basis  as  the  cross¬ 
correlation  between  a  DIS  output  pulse  and  the  reference  of  the  transmitted  pulse.  This 
cross-correlation  is  performed  with  the  product  of  the  FFT  of  both  signals.  The  FFT  of 
the  transmitted  pulse  reference  for  each  pulse  is  obtained  from  the pc_ref.mat  Matlab  data 
file  that  is  created  during  the  execution  of  the  mathost_v5.m  file.  A  plot  of  this  FFT  sig¬ 
nal  for  one  pulse  is  shown  in  Figure  26.  The  FFT  of  a  DIS  output  pulse  is  obtained  by 
calculating  the  FFT  of  the  respective  row  of  the  finalAdderOut  matrix.  A  plot  of  this  FFT 
signal  is  displayed  in  Figure  27.  Both  FFT  output  vectors  are  multiplied  element  by  ele¬ 
ment  and  saved  under  a  temporary  vector  imp.  This  Fourier  product  signal  is  shown  in 
Figure  28. 


Figure  26  FFT  of  the  Chirp  Reference  Signal. 
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x10^ 


Figure  27  FFT  of  the  DIS  Output  Signal. 


x10® 


Figure  28  FFT  Product  Signal  at  the  Fast  Cross-correlation  Output. 

The  output  of  the  cross-correlation  for  one  pulse  is  obtained  by  calculating  the  in¬ 


verse  FFT  of  the  FFT  product,  as  in  (3.24).  The  inverse  FFT  of  the  temp  vector  is  saved 
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under  a  second  temporary  vector  tmp2.  This  vector  contains  a  one  dimensional  range 
profile  for  one  returned  pulse.  A  plot  of  the  tmp2  vector  is  presented  in  Figure  29,  dem¬ 
onstrating  that  the  target  range  profile  can  be  extracted  from  the  first  Nr  samples,  where 
Nr  is  the  target  extent  in  range  bins.  The  range  profile  vector  for  each  pulse  is  extracted 
and  saved  as  a  row  in  a  new  matrix  labeled  as  priRgMapShift. 


The  cross-correlation  process  implemented  as  an  FFT  product  followed  by  an  in¬ 


verse  FFT  is  repeated  for  every  returned  pulse  using  a  loop  structure.  The  end  result  is 
the  priRgMapShift  matrix  containing  a  collection  of  range  profiles  from  different  pulses. 


The  code  lines  used  in  the  range  compression  model  are  listed  below: 


%  ISAR  PROCESSING  -  RANGE  COMPRESSION 
load  pc_ref 

priRgMapShift  =  zeros(fix(nunntaps*9/8),nDopplerCell): 
%  Range  compression: 

pcRefMapShift  =  fft(finalAdderOut.',2*rg_pts-1).': 
for  idx  =  1  :nDopplerCell, 

tmp  =  cref(:,idx).*pcRefMapShift(idx,:).': 
tmp2  =  ifft(tmp); 

priRgMapShift(:,idx)  =  tmp2(  1:fix(numtaps*9/8) ); 
end 


%  Loads  FFT  of  ISAR  reference 


%  FFT  of  DIS  output 
%  For  all  pulses. 

%  Fast  correlat.(FFT  product) 

%  InvFFT  to  obtain  range  profile 
%  Collect  range  profile  as  a  row 
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Figure  30  shows  a  collection  of  16  range  profiles  extracted  from  the  simulation 
model.  They  were  obtained  by  plotting  the  amplitude  of  the  first  16  rows  from  the 
priRgMapShift  complex  matrix. 


w 

A 

A 

fi 

Ih 

Figure  30  Range  Profiles  from  16  Different  Pulses  (Vertical  Axis  Not  Scaled). 

The  magnitude  of  the  priRgMapShift  complex  matrix  can  be  plotted  in  order  to 
view  the  range  profiles  vertically  stacked.  Figure  3 1  presents  this  surface  plot  with  the 
intensities  mapped  in  a  color  scale.  Due  to  scaling  differences,  only  the  highest  magni¬ 
tude  range  profiles  stand  out  in  the  figure. 

9.  ISAR  Azimuth  Compression  and  Image  Generation 

The  priRgMapShift  matrix  contains  the  one-dimensional  profile  from  each  re¬ 
turned  pulse.  Thus,  this  matrix  can  be  interpreted  as  a  time  history  of  range  profiles, 
where  time  is  represented  by  the  pulse  index.  The  azimuth  compression  process  consists 
of  an  FFT  along  the  time  dimension  in  order  to  perform  a  conversion  to  the  frequency 
domain.  This  is  equivalent  to  an  FFT  of  the  priRgMapShift  matrix  calculated  along  the 
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column  dimension  (pulse  indices).  The  scatterers  at  each  range  bin  are  rearranged  in 
Doppler  resolution  cells  depending  on  their  frequency  component.  The  end  result  is  a 
matrix  whose  magnitude  distribution  matches  the  range  Doppler  image  of  the  false  target 
that  is  loaded  from  the  Doppler  template  fde.  The  magnitude  matrix  is  stored  under  the 
name  of  dpRgMapShiftMOD.  This  matrix  is  truncated  to  the  first  Nr  columns  (target  ex¬ 
tent  in  range  cells)  and  Np  rows  (number  of  integrated  pulses)  and  saved  under  the  matrix 
name  quantimage.  The  code  lines  involved  in  this  azimuth  compression  process  and 
truncation  of  the  final  image  matrix  are  listed  below: 


%  ISAR  PROCESSING  -  AZIMUTH  COMPRESSION 

dpRgMapShiftMOD  =  abs(fft(priRgMapShift.',nDopplerCell)): 

quantimage  =  dpRgMapShiftMOD(1  :nDopplerCell*100/prf,1  :fix(numtaps/8*9)): 

save  finresoldata  quantimage  %  Saves  Range  Doppler  profile 


x10^ 


Range  Cells. 


Figure  3 1  Surface  Plot  of  Range  Profile  Matrix. 


To  plot  the  target  range  Doppler  image  using  frequency  values  in  the  azimuth 
axis,  it  is  necessary  to  relate  the  Doppler  filter  number  n  with  the  Doppler  frequency  f. 
This  relation  is  given  by: 


PRF 

n 


(4.12) 
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Using  a  Doppler  filter  vector  from  1  to  128  and  applying  (4.12),  a  Doppler  frequency 
vector  is  obtained  and  saved  under  the  name  dopp.  Similarly,  a  range  axis  vector  is  de¬ 
fined  from  1  to  the  number  of  columns  of  the  quantimage  matrix.  This  vector  is  saved 
under  the  name  rang.  The  range  Doppler  image  for  a  target  of  32  range  cells  and  128 
integrated  pulses  is  presented  in  Figure  32  and  can  be  compared  with  the  intended  target 
profile  in  Figure  19.  The  code  lines  used  to  define  the  axes  and  plot  the  target  image  are 
listed  below: 

%  ISAR  PROCESSING  -  IMAGE  DISPLAY 

dopp  =  (1  :size(quantinnage,1  ))*prf/nDopplerCell:  %  Doppler  axis 

rang  =  1:size(quantimage,2);  %  Range  axis 

colormap(hot(100)) 

surf(rang, dopp, quantimage),  shading  interp,  lighting  phong,  view(0,90): 


5  10  15  20  25  30  35 

Down  Range  Cells 


Figure  32  Range  Doppler  Image  for  a  False  Target  of  32  Range  Cells. 

D.  IDEAL  INFINITE  RESOLUTION  MODEL 
I.  Definition 

The  ideal  model  is  arbitrarily  defined  as  the  DIS  model  that  has  been  modified  to 
include  the  following  assumptions: 

•  The  phase  coefficients  are  not  limited  to  integer  values  (no  quantization  is 
used) 
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•  The  phase  samples  generated  by  the  DRFM  are  not  limited  to  integer  val¬ 
ues  (no  quantization  is  used) 

•  The  LUT  stage  is  replaced  by  a  Matlab  calculation  of  Sine  and  Cosine 

•  The  gain  coefficients  are  linearly  assigned  assuming  that  the  DIS  uses  a 

gain  multiplier  instead  of  a  binary  shifter.  Thus,  coefficients  are  not  lim¬ 
ited  to  integer  values  (no  quantization  is  used). 

•  The  limiting  resolution  of  0.25  due  to  the  LSB  at  the  I  and  Q  channel  for¬ 
mats  is  non-existent 

•  There  is  no  overflow  at  the  output  adder 

•  The  DAC  used  at  the  output  channels  of  the  DIS  has  infinite  resolution 

It  must  be  pointed  out  that  the  term  infinite  resolution  is  idealistic  in  the  sense  that 
the  simulation  is  always  bounded  to  the  resolution  of  the  MATLAB  program  and  the 
computer. 

2.  Software  Implementation 

Copies  of  the  three  main  simulation  files  are  made  and  renamed  as  ex- 
tract_y5 _inf_res.m,  mathost_v5 _inf_res.m  and  simhwchk_y5 _inf_res.m.  These  files  are 
moved  to  a  new  directory.  This  subsection  describes  the  code  that  must  be  applied  to 
simulate  the  ideal  model. 

a.  Phase  Coefficients 

The  phase  coefficients  for  the  DIS  model  are  obtained  from  Equation  (3.5) 
resulting  in  integer  values  from  0  to  31.  In  contrast,  the  coefficients  used  in  the  ideal 
model  are  not  rounded.  Thus,  the  symbol  [  J  is  not  used.  The  resulting  phase  coefficient 

values  range  from  0  to  31.999.  To  generate  phase  coefficients  for  the  ideal  model,  the 
following  loop  structure  from  the  extract_v5 _inf_res.m  file  are  disabled: 
for  i8  =  1  :dp_pts, 

for  i9  =  1  :numtaps,  %  5-bit  quantization  in  32  levels 

phasecoeff(i9,i8)  =  fix(renn(phaseinc(i9,i8)+32,32)): 
end 
end 

and  substituted  by  the  code  line: 

phasecoeff  =  phaseinc;  %  Added  For  Infinite  Resolution 
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b.  DRFM  Phase 

Similarly  to  the  phase  coefficient  case,  the  DRFM  phase  samples  defined 
in  (3.14)  are  not  rounded  for  the  ideal  case.  The  resulting  DRFM  phase  sample  values 
also  range  from  0  to  31.999.  The  code  change  required  for  the  ideal  model  is  the  elimina¬ 
tion  of  the  fix  rounding  function  from  the  following  mathost_y5_inf_res.m  code  line: 

pri_rg_phaseq(idx1 ,:)  =fix(oldphase/p):  %  Phase  quantized  in  2"  leveis 

c.  Sine  and  Cosine  LUT 

The  sine  and  cosine  of  the  phase  modulation  output  in  the  ideal  model  is 
achieved  using  the  MATLAB  sin  and  cos  functions.  The  input  phase  values  need  to  be 
previously  converted  from  the  0  to  31.99  scale  to  radian  values.  The  ideal  model  LUT 
output  is  obtained  by  enabling  the  simhwchk_v5_inf_res.m  file  code  line: 

lutOut=cos(phAddOut*2*pi/32)+sqrt(-1)*sin(phAddOuf2*pi/32):  %  ideai  LUT 
and  disabling  the  code  line: 

tmp  =  renn(phAddOut+depthLUT,depthLUT)+1:  %  Phase  iimited  between  0-31 

lutOut  =  cosine8(tnnp)  +  j*sine8(tnnp):  %  8-bit  LUT  tabie  output 

d.  Gain  Coefficients 

The  gain  coefficients  in  the  ideal  model  are  linearly  assigned  by  multiply¬ 
ing  the  normalized  magnitude  by  the  maximum  shifter  gain  of  1024.  This  ideal  gain  as¬ 
signment  can  be  expressed  as: 

2g(''.«)  =  1024  r^(r,n)  (4.13) 

where  is  the  normalized  magnitude  of  the  superposition  of  scatterers  in  (3.6). 

Equation  (4.13)  is  implemented  in  the  simhwchk_v5 _inf_res.m  file  by  enabling  the  fol¬ 
lowing  code  line: 

gain  =  1024*amp/nnax(nnax(annp));  %  Added  for  infinite  resolution 

and  disabling  the  following  loop  structure: 

gain  =  ones(nunntaps,dp_pts):  %  Default  gain  coefficient  is  1 . 

for  i1  =  Lnumtaps, 
for  i2  =  1  :dp_pts, 

normAmp  =  amp(i1  ,i2)/nnax(nnax(amp)):  %  Magnitudes  are  normalized, 

for  s=1:10, 

if  normAmp  >  (2''-(1 1-s))*3/2, 
gain(i1,i2)  =  2''s: 

end 
end 
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%  Thresh:  0.001 5.. ..,0.375,0.75 

%Gain:  2,4,8,16, . 512,1024 

%  Assigns  1 1  power-of-2  gains 


end 

end 

e.  I  and  Q  Channel  Resolution  Phase  Coefficients 
The  resolution  of  the  I  and  Q  ehannels  in  the  ideal  model  is  not  restrieted 
by  a  bit  format.  Thus,  the  implementation  to  simulate  the  use  of  an  LSB  must  be  elimi¬ 
nated,  which  is  accomplished  in  the  simhwchk_y5 _inf_res.m  file  by  enabling  the  code 
line: 

tapOut(inPlsCnt+delay,idx3)  =  gainOut(idx3): 
and  disabling  the  code  line: 

tapOut(inPlsCnt+delay,idx3)  =  fix(gainOut(idx3)/o_res)*o_res: 


f  Output  Adder  Overflow 

The  output  adder  in  the  ideal  model  has  an  unlimited  output  amplitude. 
Thus,  the  overflow  simulation  included  in  the  simhwchk_v5 _inf_res.m  file  must  be  elimi¬ 
nated.  This  code  change  is  achieved  by  disabling  the  following  code  lines.- 


lout=real(sum1); 

%  Overflow  simulation  using  +/-limit 

Qout=innag(sunn1); 

if  lout  >=  limit, 

%  Positive  overflow  (in-phase) 

sum1  =  lout-2*limit  +  j*Qout: 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

elseif  lout  <  -limit, 

%  Negative  overflow  (quad-phase) 

sum1  =  lout+2*limit  +  j*Qout: 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

end 

if  Qout  >=  limit. 

%  Positive  overflow  (in-phase) 

sum1  =  lout  +  j*(Qout-2*limit): 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

elseif  Qout  <  -limit. 

%  Negative  overflow  (quad-phase) 

sum1  =  lout  +  j*(Qout+2*limit); 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

end 

peak=max(max(abs(lout),abs(Qout)),peak): 

%  End  of  overflow  simulation. 

g.  DIS  output DAC 

The  ideal  model  assumes  that  the  DIS  output  channels  are  passed  through 

ideal  DACs  with  infinite  bit  resolution.  This  model  can  be  implemented  by  disabling  the 

following  code  lines  from  the  simhwchk_v5 _inf_res.m  file: 

dac_res  =  0.25;  %  DAC  resolution  0.25  for  16bit 

finalAdderOut  =  fix(finalAdderOut/dac_res)*dac_res: 
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h.  Storage  of  the  Ideal  Model  Target  Image 

The  ideal  model  range  Doppler  image  that  is  obtained  by  the 
simhwchk_y5_inf_res.m  file  is  saved  under  the  variable  name  image  in  the  in- 
fresoldata.mat  data  file.  This  change  is  implemented  by  enabling  the  following  code 
lines  from  the  simhwchk_v5_inf_res.m  file: 

quantimage  =  dpRgMapShiftMOD(1  :nDopplerCell*100/prf,1  :fix(numtaps/8*9)): 

save  infresoldata  image  %  Saves  Range  Doppler  profile 

and  disabling  the  code  lines: 

quantimage  =  dpRgMapShiftMOD(1  :nDopplerCell*100/prf,1  :fix(numtaps/8*9)): 

save  finresoldata  quantimage  %  Saves  Range  Doppler  profile 

E.  CHAPTER  SUMMARY 

The  simulation  model  presented  in  this  chapter  is  able  to  generate  the  DIS  input 
signals,  simulate  the  signal  processing  that  takes  place  in  the  chip  and  produce  the  jam¬ 
ming  output  samples.  The  model  also  simulates  an  ISAR  processor  that  takes  the  DIS 
output  samples  to  produce  a  range  Doppler  image  of  the  target.  The  simulation  of  the 
DIS  signal  processing  incorporated  the  use  of  the  different  bit  formats,  in  order  to  prop¬ 
erly  recreate  the  quantization  and  resolution  errors  involved  in  each  processing  stage. 

A  second  version  of  the  DIS  model  was  created  by  using  full  resolution  signals 
and  ignoring  the  quantization  effects  derived  from  the  use  of  finite-bit  formats.  In  the 
next  chapter,  this  ideal  model  is  used  as  a  reference  to  compare  output  images  in  order  to 
study  the  DIS  hardware  resolution  and  quantization  effects. 
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V.  DIS  ARCHITECTURE  QUANTIZATION  AND  OVERFLOW 

ANALYSIS 


The  simulation  algorithm  described  in  Chapter  IV  can  be  used  to  generate  range- 
Doppler  images  from  a  variety  of  DIS  design  alternatives.  This  chapter  studies  the  per¬ 
formance  of  different  DIS  quantization  (different  hardware  implementations).  The  analy¬ 
sis  is  performed  by  comparing  the  output  image  with  the  ideal  model.  The  effects  of 
individual  DIS  hardware  stages  are  isolated  and  measured  in  terms  of  output  image  qual¬ 
ity.  As  a  result,  the  DIS  hardware  contributions  to  the  overall  image  error  are  identified. 

A.  IMAGE  ANALYSIS  TOOLS 

Two  processes  are  developed  to  compare  the  hardware  image  model  with  the 
ideal  image  model.  These  image  analysis  tools  are  the  noise  calculation  and  the  two- 
dimensional  correlation  coefficient.  Both  methods  identify  the  images  as  matrices  of  in¬ 
tensities  of  the  same  dimensions  and  produce  an  index  to  quantify  matrix  similarity. 

1.  Mean  Square  Error  (MSE)  Calculation 

The  MSE  method  compares  the  image  distribution  of  two  images  by  obtaining  the 
average  difference  between  squared  intensities  at  equivalent  locations.  The  first  step  be¬ 
fore  comparing  image  intensities  is  to  normalize  both  matrices.  This  is  accomplished  by 
dividing  each  matrix  by  the  highest  intensity  value  and  obtaining  new  values  ranging 
from  0  to  1.  Having  both  matrices  under  this  condition,  the  ideal  matrix  can  be  sub¬ 
tracted  from  the  real  matrix,  giving  the  following  error  voltage: 

(5,1) 

ideal)  ^^A^DIs) 

If  the  matrix  intensities  C^^^^^and  in  (5.1)  are  understood  as  voltage  values, 
the  subtraction  of  these  two  normalized  matrices  is  equivalent  to  a  voltage  noise  error. 
The  deviation  of  the  real  image  from  the  ideal  image  is  due  to  the  quantization 
and  overflow  of  the  bit  formats. 

The  error  voltage  in  (5.1)  is  squared  and  averaged  over  all  the  pixel  positions  of 
the  image.  Thus,  the  result  is  equivalent  to  the  average  noise  energy  and  is  given  by: 
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e  = 


1 


(5.2) 


m=0  ^=0 


where  Np  is  the  number  of  image  matrix  rows  (Doppler  cells)  and  Nr  is  the  number  of  im¬ 
age  matrix  columns  (range  bins).  The  average  error  energy  e  will  be  referred  as  the  MSB 
coefficient.  This  quantity  ranges  from  0  to  1 ,  with  values  very  close  to  0  for  the  current 
DIS  model  (low  quantization  noise). 

2.  Image  Correlation  Coefficient 

A  two-dimensional  correlation  coefficient  is  a  mathematical  tool  that  can  be  used 
as  a  measure  of  similarity  between  two-matrices  of  the  same  size.  As  with  the  MSB  coef¬ 
ficient  method,  both  the  ideal  and  test  images  can  be  interpreted  as  matrices  of  voltage 
intensities.  The  first  step  of  this  calculation  is  to  subtract  each  matrix’s  average  intensity. 
This  is  equivalent  to  taking  the  DC  voltage  component  out  of  the  signal,  and  leaving  both 
images  with  bipolar  intensities  and  a  zero  mean. 

A  new  matrix  Ri  is  formed  with  the  element-by-element  product  of  the  two  zero 
mean  matrices: 


where  the  over  bar  denotes  the  mean.  Two  additional  matrices  R2  and  are  formed  as: 


R^{m,k)  =  {m,k)-  f 

(5.4) 

R^{m,k)  =  {^C^js{m,k)-C^rs)  ■ 

(5.5) 

For  each  of  these  three  new  matrices,  a  sum  of  all  elements  can  be  obtained  as: 


h  = 


^2  = 


r,  = 


Nr-lNp-l 

^^R^{m,k) 

m=0  k=0 

(5.6) 

Np-lNp-l 

Zj  Zj^2{fn,k) 

m=0  k=0 

(5.7) 

Nr-lNp-l 

Z  Z,R3[m,k). 

m=0  ^=0 

(5.8) 
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The  ratio  of  the  ri  coefficient  in  (5.6)  can  be  divided  by  the  square  root  of  the  product  of 
r2  and  r-^  in  (5.7)  and  (5.8)  to  obtain  a  two-dimensional  correlation  coefficient  of  the  two 
original  image  matrices  and  The  resulting  two-dimensional  coefficient  can  be 
summarized  as  [19]  : 


r  = 


6 


IS  [c ideal  ) 


m=0  k=0 


^2  r, 


(5.9) 


ZZ  [Cicleal{m^k)-C,,^^,  ')  ZZ(c„(»..q-c») 


m=0  k=0 


m=0  k=0 


The  correlation  coefficient  in  (5.9)  can  be  implemented  in  MATLAB  with  the 
function  corr2(A,B),  where  A  and  B  are  matrices  of  the  same  size.  This  correlation  func¬ 
tion  can  be  used  for  many  template-matching  applications,  such  as  face  detection  or 
speech  recognition  [20].  The  two-dimensional  correlation  coefficient  values  range  from 
0  to  1,  with  values  very  close  to  1  for  the  DIS  models.  This  indicates  the  image  is  very 
similar  to  the  ideal  model  image. 

3.  Image  Analysis  Routine 

A  simple  image  analysis  routine  is  performed  during  the  DIS  simulation  immedi¬ 
ately  after  the  ISAR  range-Doppler  image  generation.  The  analysis  routine  compares  the 
recently  generated  image  with  the  ideal-case  image  obtained  from  a  previous  simulation. 
The  comparison  is  in  terms  of  the  MSB  coefficient  and  the  two-dimensional  correlation 
coefficient  analysis. 

The  image  analysis  routine  is  listed  in  Appendix  A  and  is  saved  as  the  im¬ 
age  _analy  sis. m  file.  This  routine  assumes  that  the  real  and  ideal  images  have  been  stored 
as  matrices  of  the  same  size.  The  real  image  is  loaded  from  the  quantimage  variable 
stored  in  the  finresoldata.dat  data  file  and  generated  by  the  simhwchk_v5.m  file.  The 
ideal  image  is  loaded  from  the  image  variable  stored  in  the  infresoldata.dat  data  file  and 
generated  by  the  simhwchk_v5 _inf_res.m  file.  Once  loaded,  both  image  matrices  are 
normalized.  The  following  code  lines  implement  this  intensity  normalization: 

normimage  =  image/nnax(nnax(image)):  %  Normalizes  real  image. 

normquantimage  =  quantimage/max(max(quantimage)):  %  Normalizes  ideal  image. 
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The  normalized  matrices  are  subtracted  and  the  resulting  matrix  is  squared  at  each  ele¬ 
ment  location.  The  matrix  containing  the  squared  errors  is  saved  as  noiseEnergy.  The 
squared  error  at  each  pixel  location  is  found  with  the  following  loop  structure: 

for  nn=1:dp_pts/2, 

for  n=1  :numtaps*1 .0625, 

noiseEnergy(m,n)  =  (normimage(m,n)-normquantimage(m,n))''2;  %  Pixel  square  error 
end; 

end; 

The  MSB  coefficient  is  then  found  as  the  average  of  all  the  elements  of  the  noiseEnergy 
matrix.  The  calculation  is  implemented  summing  all  the  elements  of  noiseEnergy  and 
dividing  the  sum  with  the  product  of  the  number  of  rows  and  the  number  of  columns. 
This  average  is  calculated  in  the  following  code  lines: 

[a,b]  =  size(noiseEnergy); 

noise_ave  =  sum(sum(noiseEnergy))/(a*b)  %  MSE  coefficient. 

The  two-dimensional  correlation  coefficient  is  implemented  using  the  Matlab 
function  corr2( ),  using  the  real  and  the  ideal  image  matrices  as  the  two  parameters.  This 
Matlab  function  implements  the  calculation  presented  in  Equation  (5.9).  As  shown  in  the 
following  code  line,  the  correlation  coefficient  is  stored  as  corel_coeff. 

corel_coeff  =  corr2(image,quantinnage)  %  2D  correlation  coeff. 

Finally,  the  software  routine  plots  the  image  distribution  of  the  squared  error 
{noiseEnergy  matrix),  with  the  purpose  of  identifying  possible  concentrations  of  errors  at 
specific  range  bins  or  Doppler  cells. 

B.  PHASE  QUANTIZATION  ANALYSIS 

This  section  analyzes  an  ideal  model  that  is  modified  to  implement  the  DIS  phase 
modulation  stage  using  5-bit  phase  formats.  The  three  files  containing  the  ideal  model 
are  copied  to  new  files  under  the  names  of  extract _para_ship_v5_isol_quant.m, 
mathost_v5_isol_quant.m  and  simhwchk_v5 _isol_quant.m.  This  last  file  is  modified  in 
order  to  store  the  non-ideal  output  image  as  the  variable  quantimage  in  the  finre- 
soldata.dat  data  file. 

I.  Identification  of  Error  Sources 

The  intercepted  ISAR  signal  is  fed  to  the  DIS  in  terms  of  DRFM  generated  phase 
samples  as  described  in  (3.14).  This  collection  of  phase  samples  is  represented  by  the 
unsigned  5-bit  format  in  Table  2.  The  use  of  this  format  is  a  source  of  quantization  error 
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induced  in  the  phase  information.  The  maximum  phase  error  due  to  the  5 -bit  format  is 
given  by: 

=360V2' =11.25°.  (5.10) 

The  5-bit  DRFM  phase  input  is  incremented  by  the  phase  modulation  coefficient 
of  the  same  bit  format  and  error.  Thus,  the  resulting  phase  sum  from  this  modulation  has 
a  resolution  error  ranging  from  0  to  22.5  degrees.  Assuming  the  DRFM  phase  samples 
and  the  phase  coefficients  have  a  uniform  probability  of  distribution,  the  error  of  the 
modulated  phase  should  have  a  mean  of  1 1.25  degrees.  The  incidence  of  the  phase  error 
is  reflected  in  the  sine  and  cosine  LUT  table.  At  this  next  stage,  the  phase  samples  carry¬ 
ing  errors  higher  than  the  1 1.25-degree  average  produce  an  LUT  output  off  by  one  level. 
As  shown  in  Table  8  from  Chapter  IV,  the  I  and  Q  amplitude  difference  between  adjacent 
8-bit  LUT  levels  ranges  from  0.015  to  0.195  (equivalent  to  0.7  to  10  percent  error). 

2.  Implementation  of  the  Software  Simulation 

The  extract _para_ship_v5_isol_quant.m  file  is  modified  to  add  the  use  of  5-bit  for¬ 
mats  in  the  DIS  phase  modulation  stage.  The  new  model  is  run  and  the  image  effects  are 
compared  with  the  ideal  model  image.  The  image  analysis  tools  described  in  Section  A 
are  used  to  quantify  the  isolated  effect  of  the  5 -bit  phase  quantization.  The  inclusion  of 
5 -bit  phase  modulation  coefficients  is  done  enabling  the  following  loop  structure: 

for  i8  =  1  :dp_pts, 

for  i9  =  1  :numtaps,  %  5-bit  quantization  in  32  levels 

phasecoeff(i9,i8)  =  fix(renn(phaseinc(i9,i8)+32,32)): 
end 

end 

and  disabling  the  following  line: 

phasecoeff=phaseinc:  %  Added  For  Infinite  Resolution 

The  inclusion  of  the  5-bit  format  in  the  DRFM  phase  is  done  by  adding  the  fix  rounding 
function  to  the  chirp  phase  generated  in  the  mathost_v5 _isol_quant.m  file.  This  can  be 
done  by  disabling  the  code  line: 

pri_rg_phaseq(idx1 ,:)  =oldphase/p:  %  Phase  quantized  in  2"  levels 

and  enabling  the  code  line: 

prLrg_phaseq(idx1 ,:)  =fix(oldphase/p):  %  Phase  quantized  in  2"  levels 
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3.  Simulation  Results 

The  MSB  coefficient  obtained  when  the  effects  of  the  5 -bit  phase  format  in  the 
phase  modulation  are  added  is  9.8486 -10"^  The  two-dimensional  correlation  coefficient 
obtained  for  the  same  case  is  0.99931.  The  small  value  in  the  MSB  coefficient  or  noise 
energy  indicator  and  the  correlation  coefficient  value  vary  in  the  vicinity  of  one,  indicat¬ 
ing  that  the  quantization  effects  of  the  5 -bit  phase  format  are  minimum.  In  fact,  the  dif¬ 
ferences  between  the  resulting  image  and  the  ideal  case  cannot  be  identified  with  the 
human  eye.  The  image  obtained  from  this  model  and  the  distribution  of  the  MSB  are  in¬ 
cluded  in  Figure  33. 
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Figure  33  a)  Image  Generated  by  an  Ideal  DIS  with  Phase  Quantization,  b)  Image  Dis¬ 
tribution  of  the  Phase  Quantization  MSB. 

C.  BUT  RESOLUTION  AND  SCALING  EFFECTS 
1.  Ideutificatiou  of  Error  Sources 

The  output  of  the  DIS  phase  modulation  circuit  is  converted  into  I  and  Q  ampli¬ 
tude  components  using  cosine  and  sine  calculations.  Due  to  the  high-speed  system  re¬ 
quirements,  an  LUT  table  with  8-bit  resolution  is  used.  The  BUT  is  a  5-bit  decoder, 
where  each  of  the  32  possible  phases  that  can  be  represented  by  the  format  in  Table  2  has 
an  equivalent  sine  and  cosine  value.  The  errors  induced  by  the  use  of  an  LUT  with  8  bits 
of  resolution  to  calculate  cosine  and  sine  can  be  classified  as  resolution  or  quantization 
errors. 
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The  LUT  resolution  errors  originate  from  the  output  format  of  the  I  and  Q  chan¬ 
nels.  These  LUT  outputs  have  an  8 -bit  two’s  complement  format  with  a  sign  bit  and  7 
fractional  bits.  Thus,  the  sine  and  cosine  calculations  have  a  maximum  resolution  error 
of  2"’.  This  value  in  the  output  scale  from  - 1  to  +1  represents  a  0.4  percent  error.  If  the 
8-bit  LUT  table  is  replaced  by  a  7,  6,  5,  4  or  3-bit  table,  the  maximum  resolution  error 
increases  to  2“*’ ,  2“^ ,  2“'* ,  2“^  and  2“^,  respectively.  Also,  these  values  represent  a  0.8, 
1.6,  3.2,  6.4  and  12.8  percent  error,  respectively. 

2.  Implementation  of  the  Software  Simulation 

The  LUT  quantization  effects  are  isolated  by  reversing  the  changes  done  in  the 
phase  quantization  analysis  to  return  to  the  ideal  model  and  then  adding  an  LUT  to  re¬ 
place  the  sine  and  cosine  calculations.  The  inclusion  of  the  LUT  table  is  implemented 
changing  the  code  in  the  simhwchk_v5_isol_quant.m  file  to  disable  line: 

lutOut=cos(phAddOut*2*pi/32)+sqrt(-1)*sin(phAddOuf2*pi/32):  %  Ideal  LUT 

and  to  enable  the  code  lines  listed  below: 

tmp  =  renn(phAddOut+depthLUT,depthLUT)+1;  %  Phase  limited  between  0-31 

lutOut  =  cosine8(tmp)  +  j*sine8(tmp):  %  8-bit  LUT  table  output 

where  cosineS  and  sineS  are  32-element  array  variables  storing  the  contents  of  the  8-bit 
LUT.  An  additional  code  line  must  be  added  to  round  up  the  infinite  resolution  phase 
(non-integers)  to  be  used  as  integer  addresses  for  the  LUT.  The  following  code  line  must 
be  added  immediately  before  the  line  defining  the  lutOut  vector: 

tmp  =  fix(tmp):  %  Needed  to  isolate  LUT  effects 

3.  Simulation  Results 

The  new  model  was  run  and  the  results  are  similar  to  the  phase  quantization  case. 
The  8-bit  LUT  quantization  showed  a  very  small  MSB  coefficient  (6. 1069  10”®)  and  a 
correlation  coefficient  very  close  to  one  (0.999916).  Thus,  the  effects  of  the  8-bit  LUT 
quantization  can  be  considered  minimal.  The  image  generated  by  the  model  and  the 
MSB  distribution  is  shown  in  Figure  35. 
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Figure  34  Plot  of  the  Sine  Curve  Using  32  Phase  Values  in  Different  LUTs. 
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Figure  35  a)  Image  Generated  by  an  Ideal  DIS  with  LUT  Quantization,  b)  Image  Distri¬ 
bution  of  the  LUT  Quantization  MSB. 

To  evaluate  the  effect  of  additional  LUT  configurations,  the  arrays  cosineS  and 
sineS  are  replaced  to  implement  LUT  options  from  7  bits  through  3  bits.  The  results  are 
summarized  in  Table  9.  The  DIS  images  for  the  different  LUT  cases  are  not  presented 
here  since  they  are  too  similar  to  visually  recognize  their  differences. 
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LUT  resolution 

MSE  coefficient 

2D  Correl.Coefficient 

8-Bit 

6. 1069. 10“*’ 

0.9999160 

7-Bit 

6. 1382- 10“*’ 

0.9999155 

6-Bit 

6.2004. 10“*’ 

0.9999149 

5-Bit 

6.4654.10“*’ 

0.9999112 

4-Bit 

6.7981.10“*’ 

0.9999057 

3-Bit 

1.1416.10“*’ 

0.9998420 

Table  9  Image  Quality  Statistics  for  Different  LUT  Resolutions. 

In  summary,  the  simulations  showed  that  the  quantization  effects  of  the  LUTs  are 
negligible.  Even  at  the  extreme  case  of  cosine/sine  calculations  using  a  3-bit  resolution 
LUT  where  only  seven  output  amplitude  levels  are  obtained,  the  effects  are  minimal. 

D.  GAIN  QUANTIZATION  ANALYSIS 

1.  Identification  of  Error  Sources 

In  the  gain  modulation  stage  of  the  DIS  pipeline,  the  I  and  Q  outputs  of  the  LUT 
are  multiplied  by  a  power-of-two  gain  (2^'^'”"*).  This  multiplication  is  implemented  using 
a  binary  shifter  using  the  control  codes  listed  in  Table  4  from  Chapter  III.  As  described 
in  that  chapter,  the  generation  of  these  coefficients  is  performed  for  each  range-bin/pulse 
combination.  The  default  quantization  to  obtain  the  gain  coefficients  follows  the  scheme 
shown  in  Table  1,  where  normalized  magnitudes  in  the  range  from  0  to  1  are  respectively 
assigned  gains  from  2^  to  2^*.  This  quantization  scheme  is  shown  as  a  staircase  shape 
curve  in  Figure  36.  Ideally,  the  gain  values  should  be  linearly  related  to  the  normalized 
magnitudes,  and  thus  following  the  diagonal  line  of  that  figure.  However,  the  exponen¬ 
tially  nature  of  the  gain  shifter  limits  the  gain  assignment  to  exponentially  distributed 
levels. 


2.  Implementation  of  the  Software  Simulation 

The  gain  quantization  effects  are  isolated  by  reverting  the  previous  changes  to  re¬ 
turn  to  the  ideal  model  and  adding  new  code  modifications  to  the  ex¬ 
tract _para_ship_v5Jsol_quant.m  file.  The  following  code  line  is  disabled  to  eliminate  the 
ideal  gain  quantization: 

gain  =  1024*amp/nnax(nnax(annp));  %  Added  for  infinite  resolution 
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and  the  following  loop  structure  is  enabled  to  return  to  the  default  gain  quantization 
scheme: 


gain  =  ones(nunntaps,dp_pts): 
for  i1  =  hnumtaps, 
for  i2  =  1  :dp_pts, 

normAmp  =  amp(i1,i2)/nnax(nnax(amp)): 
for  s=1:10, 

if  normAmp  >  (2'^-(1 1-s))*3/2, 
gain(i1,i2)  =  2^s: 
end 
end 
end 
end 


%  Default  gain  coefficient  is  1 . 


%  Magnitudes  are  normalized. 

%  Thresh:  0.001 5...., 0.375, 0.75 

%Gain:  2,4,8,16, . 512,1024 

%  Assigns  1 1  power-of-two  gains. 


Figure  36  Quantization  Scheme  for  Gain  Modulation  Coefficients. 


3.  Simulation  Results 

The  new  model  is  run  yielding  an  MSB  coefficient  of  5.220-10"^  and  a  correla¬ 
tion  coefficient  of  0.993 1 1 .  The  model  image  using  the  DIS  gain  quantization  and  the 
MSB  distribution  are  shown  in  Figure  37. 
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Figure  37  a)  Image  Generated  by  an  Ideal  DIS  with  Gain  Quantization,  b)  Image  Distri¬ 
bution  of  the  Gain  Quantization  MSB. 

E.  I  AND  Q  CHANNEL  RESOLUTION  ANALYSIS 

1.  Identification  of  Error  Sources 

The  format  of  the  I  and  Q  BUT  channels  after  truncation  of  the  5  LSBs  has  a  13- 
bit  two’s  complement  format  as  shown  in  Table  5.  The  LSB  of  that  format  is  2“^ .  Thus, 
all  the  possible  I  and  Q  magnitudes  are  multiples  of  0.25.  The  resolution  error  in  the  DIS 
performance  is  quantified  in  this  section. 

2.  Implementation  of  the  Software  Simulation 

After  reverting  the  changes  implemented  in  the  previous  model  to  return  to  the 
ideal  model,  new  code  modifications  are  performed  to  the  simhwchk_y5 _isol_quant.m  file  in 
order  to  implement  a  truncation  by  multiples  of  0.25.  That  truncation  is  included  by  ena¬ 
bling  the  following  code  line: 

tapOut(inPlsCnt+delay,idx3)  =  fix(gainOut(idx3)/o_res)*o_res: 
and  by  disabling  the  following  code  line: 

tapOut(inPlsCnt+delay,idx3)  =  gainOut(idx3); 
where  the  variable  o  res  is  the  bit  format  LSB  with  a  value  of  0.25. 
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3.  Simulation  Results 

The  model  is  ran  with  a  32-range-bin  target  and  128  radar  pulses.  The  image  gen¬ 
erated  had  an  MSB  coefficient  of  0.023378  •10'^  and  the  two-dimensional  correlation 
with  the  ideal  image  had  a  coefficient  of  0.9999998.  The  simulation  results  indicate  that 
the  resolution  effects  of  the  LSB  in  the  I  and  Q  channels  are  negligible.  The  image  ob¬ 
tained  in  this  simulation  is  shown  along  with  the  MSB  distribution  in  Figure  38. 
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Figure  38  a)  Image  Generated  by  an  Ideal  DIS  with  BSB  Resolution  Bffects  in  the  I  and 
Q  Channels,  b)  Image  Distribution  of  the  MSB. 

F.  OUTPUT  ADDER  OVERFLOW  ANALYSIS 
1.  Ideutificatiou  of  Error  Sources 

The  two’s  complement  output  format  of  the  16-bit  output  adder  can  represent  val¬ 
ues  from  -2'^  to  2'^  -2“^  (-8192  to  8191.75),  as  shown  in  Table  6.  Although  the  out¬ 
put  values  that  are  originated  in  each  range  bin  vary  from  - 1028  to  1027.75  (cases  for 
maximum  gain),  it  is  feasible  that  the  cumulative  sum  of  range  bins  exceeds  the  -8192 
and  8191.75  output  adder  limits.  An  overflow  occurs  in  the  output  adder  when  the  cumu¬ 
lative  sum  of  the  range-bin  values  goes  beyond  these  limits. 

The  main  effect  of  the  output  adder  overflow  is  that  the  sum  of  the  most  signifi¬ 
cant  bits  (2*^)  produces  a  carry  that  changes  the  sign  of  the  two’s  complement  format. 
This  sign  change  effect  was  studied  in  order  to  predict  the  output  value  for  overflow 
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cases.  It  was  found  that  the  result  of  an  overflow  due  to  the  sum  of  a  value  A  and  a  value 
B  can  be  expressed  as: 


(A+BL,^„={A+B)  +  2.8192  if(.4  +  S)<-8192  (5.11) 

(A+Bk,,„=(A+B)-2.8192  if  (J  +  i!)>+8191.75  (5.12) 


where  (5.11)  and  (5.12)  represent  a  negative  and  a  positive  overflow  respectively. 

2.  Implementation  of  the  Software  Simulation 

The  code  changes  performed  for  the  previous  model  are  reversed  to  return  to  the 
ideal  mode.  Next,  the  extract _para_ship_v5_isol_quant.m  file  code  is  changed  to  imple¬ 
ment  equations  (5.1 1)  and  (5.12)  by  enabling  the  following  code  lines: 


lout=real(sunn1): 

%  Overflow  simulation  using  +/-limit 

Qout=imag(sum1): 

if  lout  >=  limit, 

%  Positive  oflow  (in-phase) 

sum1  =  lout  -  2*limit+j*Qout; 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

elseif  lout  <  -limit, 

%  Negative  oflow  (quad-pha) 

sum1  =  lout  +  2*limit+j*Qout: 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

end 

if  Qout  >=  limit. 

%  Positive  oflow  (in-phase) 

sum1  =  lout  +  j*(Qout-2*limit): 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

elseif  Qout  <  -limit. 

%  Negative  oflow  (quad-pha) 

sum1  =  lout  +  j*(Qout+2*limit); 

%  Adder  output  w/overflow 

overflow  =  overflow  +  1 ; 

%  Increases  counter 

end 

peak=max(max(abs(lout),abs(Qout)),peak): 

%  End  of  oflow  simul. 

where  the  variable  limit  was  previously  set  to  8192  and  the  variable  peak  is  used  to  regis¬ 
ter  the  highest  overflow  value. 

3.  Simulation  Results 

The  simulation  of  this  model  is  run  using  a  32-range-bin  target  and  128  pulses. 
After  comparing  the  model  image  with  the  ideal  quantization  image,  the  MSB  coefficient 
obtained  was  0.43933  •  10'^  and  the  two-dimensional  correlation  coefficient  was  0.99407. 
For  this  particular  target,  the  number  of  overflows  is  22  with  a  maximum  value  of  8443.5. 
The  image  of  this  model  and  the  MSB  distribution  is  shown  in  Figure  39. 
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Figure  39  a)  Image  Generated  by  an  Ideal  DIS  with  Output  Adder  Overflow,  b)  Image 
Distribution  of  the  MSE  Due  to  Overflow. 

G.  DIS  HARDWARE  MODEL  SIMULATION 

1.  Identification  of  Resolution  Error  Sources 

This  model  simulates  the  DIS  hardware  performanee  including  all  possible  quan¬ 
tization  and  resolution  effects  contained  in  the  DIS  being  fabricated,  except  only  32 
range-bin  processors  are  being  evaluated  instead  of  512.  The  model  includes  the  use  of  a 
5-bit  phase  format,  an  8-bit  resolution  LUT,  gain  coefficient  quantization  from  Table  1, 
the  minimum  resolution  of  the  I  and  Q  channels  and  the  overflow  effects  of  the  output 
adders.  This  model  also  assumes  the  DIS  output  is  connected  to  16-bit  DACs.  Thus, 
truncation  does  not  occur  due  to  bit  format  mismatch.  The  simulation  is  performed  using 
the  defaults  simulation  files:  extract _para_ship_v5.m,  mathost_v5.m  and  simhwchk_v5.m. 

2.  Simulation  Results 

The  default  model  was  run  with  a  32-range-bin  target  and  integrating  128  pulses. 
The  simulation  yielded  an  MSE  of  0.074854  - 10'^  and  a  two-dimensional  correlation  co¬ 
efficient  of  0.990004.  The  number  of  overflow  counts  is  22  with  a  peak  value  of  8410. 

H.  RESOLUTION  EFFECTS  OF  OUTPUT  DAC  CONVERSION 

1.  Identification  of  Resolution  Error  Sources 

The  resolution  effects  of  the  I  and  Q  channels  (Section  E)  minimally  impacted  the 
DIS  performance.  Therefore,  it  may  be  expected  that  truncating  the  DIS  output  by  the 
least  significant  bits  (LSB)  would  have  the  same  negligible  effect.  The  DIS  output  chan- 
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nels  have  a  16-bit  format.  The  intention  here  is  to  simulate  the  hardware  model  of  Sec¬ 
tion  G  added  with  the  effects  of  coupling  the  DIS  output  with  a  14,  a  12  and  a  10-bit 
DAC,  and  thus  truncating  different  amounts  of  LSBs. 

2.  Implementation  of  the  Software  Simulation 

The  simhwchk_v5.m  file  is  changed  in  order  to  set  the  value  of  the  DAC  minimum 
resolution  (dac_res  variable)  according  to  Table  10.  The  implementation  of  this  code 
change  ensures  that  all  I  and  Q  amplitudes  in  the  DIS  output  are  exact  multiples  of  the 
resolution  value  dac  res. 


DAC  resolution 

Resolution  value 

16  bit 

0.25 

Mbit 

1 

12  bit 

4 

10  bit 

16 

8  bit 

64 

6  bit 

256 

Table  10  Minimum  DIS  Resolution  for  Different  DAC  Configurations. 

3.  Simulation  Results 

The  16-bit  DAC  model  is  equivalent  to  the  hardware  model  of  Section  G,  where 
the  DAC  does  not  produce  any  truncation  effect.  Changing  the  DAC  to  14-bit,  where  the 
2  DIS  output  LSBs  are  discarded,  yielded  an  MSB  coefficient  of  0.74799-10”^  and  a  two- 
dimensional  correlation  coefficient  of  0.990005.  For  12-bit  DAC,  4  output  LSBs  are  dis¬ 
carded  and  the  results  are  respectively  e  =  0.73572-10”^  and  r  =  0.990078.  For  10-bit 
DAC,  6  LSBs  are  discarded  obtaining  an  MSB  coefficient  of  e  =  0.72228-10”^  and  a  two- 
dimensional  coefficient  of  r  =  0.990082.  There  is  a  slight  apparent  image  quality  im¬ 
provement  when  reducing  the  DAC  resolution  for  some  cases.  However,  this  is  attribut¬ 
able  to  an  error  cancellation  effect  that  can  be  neglected  because  of  the  small  magnitude. 
For  an  8-bit  DAC  case,  e  =  1.200*10'^  and  r  =  0.98605.  As  shown  in  Figure  40,  only 
when  DAC  resolution  is  reduced  to  4  bits  is  there  a  visually  notable  image  deterioration 
with  r  =  0.9279.  These  results  clearly  show  that  most  of  the  vital  information  of  the  DIS 
output  is  included  in  the  4  MSBs  and  that  DACs  of  6  or  8  bits  can  be  used  without  affect¬ 
ing  the  system  performance.  As  a  final  note,  the  DAC  must  sample  at  the  rate  of  the  DIS 
clock  frequency,  which  dictates  the  range  resolution  of  the  device. 
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Figure  40  Target  Image  and  MSB  for  Different  Output  DACs:  (a)  16-Bit,  (b)  8-Bit,  (c) 
6-Bit  and  (d)  4-Bit. 


I.  CHAPTER  SUMMARY 

The  effeets  of  the  bit  format  quantization  and  resolution  on  the  DIS  output  images 
were  studied  for  different  models  in  terms  of  the  similarity  to  an  ideal  image.  The  iso- 
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lated  effects  from  the  5 -bit  phase  format,  the  sine/cosine  LUT,  the  4-bit  gain  control  code 
format,  the  13 -bit  I  and  Q  channel  format  and  the  output  adder  overflow  were  modeled 
independently.  Only  the  effects  of  gain  quantization  and  adder  overflow  errors  stood  out. 
The  rest  were  mathematically  negligible. 

The  hardware  architecture  is  modeled  considering  all  bit  format  effects  (Section 
G).  This  model  is  basically  the  true  simulation  of  the  DIS  performance.  The  effects  of 
the  combined  quantizations  on  the  image  errors  were  minimum  {r  =  0.990)  and  not  visu¬ 
ally  recognizable.  Finally,  the  output  truncating  effects  of  using  DACs  with  different  bit 
resolution  were  simulated.  From  these  results,  it  can  be  concluded  that  the  DIS  output 
information  is  concentrated  in  the  first  8  MSBs  from  the  16-bit  format.  To  support  this 
fact,  the  discarding  of  the  8  LSB  out  of  the  output  format  resulted  in  image  effects  that 
were  mathematically  negligible. 

A  summary  of  the  image  effects  of  the  different  simulated  models  is  included  in 
Table  11.  In  general,  the  DIS  simulation  models  proved  that  the  architecture  design  is 
robust  and  able  to  endure  the  effects  of  error  originated  from  the  use  of  values  truncated 
by  bit  formats. 


Effect  Simulated 

MSE  from 

Ideal  Image 

2D  Correlation 
Coefficient 

Phase  5 -bit  format 

0.0098486  10“^ 

0.999310 

Sine/Cosine  8-bit  LUT 

0.006 1069- 10“^ 

0.999916 

Sine/Cosine  7-bit  LUT 

0.006 13  82  10“^ 

0.999915 

Sine/Cosine  6-bit  LUT 

0.0062004  10-' 

0.999914 

Sine/Cosine  5 -bit  LUT 

0.0064654  10-' 

0.999911 

Sine/Cosine  4-bit  LUT 

0.006798  MO-' 

0.999905 

Sine/Cosine  3 -bit  LUT 

0.011416-10-' 

0.999842 

Gain  control  4-bit  format 

0.522000-10-' 

0.993110 

I  and  Q  13 -bit  format 

0.0000233-10-' 

0.999999 

Output  Adder  Overflow 

0.43933-10-' 

0.994074 

Hardware  model  with  1 6-bit  D AC 

0.74854-10-' 

0.990004 

Hardware  model  with  14-bit  DAC 

0.74799-10-' 

0.990005 

Hardware  model  with  12-bit  DAC 

0.73572-10-' 

0.990078 

Hardware  model  with  10-bit  DAC 

0.72228-10-' 

0.990082 

Hardware  model  with  8 -bit  DAC 

1.200-10-' 

0.986059 

Table  1 1  Quantization  and  Resolution  Error  Effects  on  the  DIS  Image. 
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VI.  DIS  PERFORMANCE  UNDER  DIFFERENT  ISAR 

PARAMETERS 


This  chapter  explains  the  ISAR  requirements  for  PRF  and  eoherent  integration 
time  and  relates  these  parameters  to  the  Doppler  resolution  obtained  at  the  ISAR  target 
image.  The  simulation  model  is  run  to  verify  the  expeeted  results  and  the  effeet  of  these 
parameters  under  a  different  set  of  values. 


A.  ISAR  PRF  AND  INTEGRATION  TIME  REQUIREMENTS 
1.  ISAR  PRF  Limits 

As  with  any  radar  design,  ISAR  PRF  must  be  low  enough  to  avoid  range  ambigui¬ 
ties  and  high  enough  to  avoid  Doppler  ambiguities.  The  limitation  for  maximum  PRF  is 
given  by  [2]  : 


2R 


FP 


(6.1) 


where  Rpp  is  the  range  extension  of  the  real  antenna  footprint.  For  ISARs  that  are  fo¬ 
cused  on  a  target  (spotlight),  the  range  window  is  relatively  short.  For  a  range  footprint 
of  1km,  the  maximum  PRF  is  150,000  pulses  per  seeond  (pps).  To  avoid  Doppler  ambi¬ 
guities,  the  PRF  must  exceed  the  maximum  Doppler  frequeney  spread.  For  airborne  ap- 
plieations,  the  minimum  ISAR  PRF  is  defined  as  [8]: 

(6-2) 

where  Vr  is  the  velocity  magnitude  of  the  aircraft  and  is  the  azimuth  antenna  real 

beamwidth.  Using  worst-ease  values  with  a  speed  of  205m/s  (400  knots),  a  beamwidth  of 
0.0175  radians  (1  degree)  and  a  wavelength  of  0.03m  (10  GHz),  the  resulting  minimum 
PRF  is  240  pps.  This  minimum  PRF  ean  be  understood  as  the  Nyquist  sampling  limit. 
Radar  pulses  are  equivalent  to  Doppler  samples  (complex).  Thus,  the  PRF,  aeting  as  a 
sampling  rate,  must  be  equal  to  the  highest  Doppler  frequeney  differenee. 
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2.  ISAR  Coherent  Integration  Time 

The  coherent  integration  time  is  the  time  period  that  the  radar  uses  to  collect  com¬ 
plex  samples  from  pulse  returns.  At  the  end  of  the  integration  time,  the  radar  samples  are 
processed  as  a  group.  In  Equation  (2.16),  it  is  demonstrated  that  the  Doppler  resolution  is 
the  inverse  of  the  coherent  integration  time. 

(6.3) 

Thus  the  azimuth  or  cross-range  resolution  improves  (decreases)  with  longer  integration 
times.  Since  the  increase  in  integration  time  is  directly  related  to  the  number  of  samples 
processed,  the  limit  is  set  by  the  signal  processing  capacity  of  the  ISAR  processor. 

B.  DIS  PERFORMANCE  UNDER  DIFFERENT  PRES  AND  INTEGRATION 

TIMES 

1.  DIS  Dependence  from  ISAR  Parameters 

The  ISAR  pulse  samples  processed  by  the  DIS  are  fed  from  a  DRFM.  Therefore, 
the  DRFM  controls  the  rate  at  which  these  samples  are  processed.  This  sampling  rate  is 
based  on  the  Nyquist  criteria  and  depends  on  the  maximum  frequency  value  of  the  analog 
signal  before  sampling.  The  PRF  defines  the  rate  at  which  the  ISAR  pulses  arrive  at  the 
DIS.  Each  radar  pulse  is  received  by  the  DIS  in  the  form  of  a  sequence  of  complex  sam¬ 
ples.  Therefore,  from  the  DIS  point  of  view,  the  ISAR  PRF  is  transparent.  In  contrast, 
the  number  of  pulses  integrated  by  the  ISAR  is  very  relevant  for  the  DIS,  since  it  deter¬ 
mines  the  number  of  modulation  coefficients  that  must  be  used  at  each  range  bin  to  syn¬ 
thesize  a  column  of  the  range-Doppler  image.  The  DIS  uses  different  modulation 
coefficients  for  every  pulse.  If  the  ISAR  receiver  is  designed  to  integrate  Np  pulses,  then 
the  DIS  must  repeat  the  set  of  modulation  coefficients  every  Np  pulses.  The  number  of 
integrated  ISAR  pulses  is  directly  related  to  the  PRF  and  integration  time  such  that: 

Np=PRFTj.  (6.4) 

Therefore,  DIS  dependence  on  the  ISAR  PRF  and  integration  time  are  tied  together  and 
cannot  be  analyzed  separately.  For  a  fixed  PRF,  it  can  be  expected  that  if  the  number  of 
integrated  pulses  is  increased,  the  integration  time  must  increase,  thus  the  Doppler  resolu- 
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tion  is  improved.  Similarly,  for  a  fixed  number  of  integrated  pulses,  if  the  PRF  is  in¬ 
creased,  the  Doppler  resolution  is  also  increased.  Table  12  shows  various  combinations 
of  PRFs  and  the  number  of  integrated  pulses,  along  with  their  equivalent  integration 
times  and  the  resulting  Doppler  resolution. 


PRF 

Number  of 
Integrated  Pulses 

Integration 

Time 

Doppler 

Resolution 

50  pps 

100  pps 

200  pps 

400  pps 

256  pulses 

128  pulses 

64  pulses 

32  pulses 

5.12  s 

0.20  Hz 

100  pps 

200  pps 

400  pps 

800  pps 

256  pulses 

128  pulses 

64  pulses 

32  pulses 

2.56  s 

0.39  Hz 

200  pps 
400  pps 
800  pps 
1600  pps 

256  pulses 

128  pulses 

64  pulses 

32  pulses 

1.28  s 

0.78  Hz 

400  pps 
800  pps 
1600  pps 
3200  pps 

256  pulses 

128  pulses 

64  pulses 

32  pulses 

0.64  s 

1.56  Hz 

800  pps 
1600  pps 
3200  pps 
6400  pps 

256  pulses 

128  pulses 

64  pulses 

32  pulses 

0.32  s 

3.12  Hz 

1600  pps 
3200  pps 
6400  pps 
12800pps 

256  pulses 

128  pulses 

64  pulses 

32  pulses 

0.16  s 

6.24  Hz 

Table  12  Doppler  Resolution  for  Various  ISAR  Parameters. 


2.  Simulation  Runs 

The  DIS  performance  is  modeled  for  each  of  the  six  resolution  cases  listed  in 
Table  12.  To  run  the  model,  the  number  of  integrated  pulses  is  set  to  128  pulses  while 
the  PRF  is  increased  from  100  pps  to  800  pps.  In  addition,  the  pulsewidth  is  adapted  for 
each  case  in  order  to  keep  a  constant  duty  cycle.  Figure  41  presents  the  target  range 
Doppler  images  for  each  resolution  case.  For  this  specific  target,  a  maximum  Doppler 
resolution  close  to  2  Hz  is  required  for  a  good  target  profile. 
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Figure  41  DIS  Image  for  Different  Doppler  Resolutions:  (a)  0.2  Hz,  (b)  0.39  Hz,  (e)  0.78 
Hz,  (d)  1.56  Hz,  (e)  3.12  Hz  and  (f)  6.24  Hz. 


3.  PRI  Multiplicity 

During  the  simulation  runs,  it  was  observed  that  the  use  of  a  duty  cycle,  in  which 
the  PRI  was  not  an  exact  multiple  of  the  pulsewidth,  did  not  properly  synthesize  the  tar¬ 
get  image.  The  default  values  for  PRF  and  PW  used  in  the  previous  chapters  were  200 
pps  and  0.5  us  respectively.  If  the  PRF  is  increased  from  the  default  value  to  201  pps,  the 
chirp  wave  is  not  properly  constructed  by  the  routine  in  the  mathost_v5.m  file.  Figure  42 
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shows  a  plot  of  the  chirp  wave  when  the  ISAR  PRF  is  201  pps.  This  chirp  pulse  suffers  a 
phase  reversal  near  three  quarters  of  the  pulsewidth.  After  this  ISAR  pulse  is  processed 
by  the  DIS,  the  resulting  output  modulated  pulse  returned  to  the  radar  is  that  presented  in 
Figure  43.  The  resulting  target  range  Doppler  image  is  not  properly  synthesized  by  the 
ISAR  routine  as  can  be  seen  in  Figure  44.  This  image  presents  ambiguities  in  the  Dop¬ 
pler  dimension,  resulting  in  the  display  of  three  strong  targets 


Pulse  samples 

Figure  42  Erroneous  Generation  of  a  Chirp  Pulse  of  0.5us  PW  and  201pps  PRF. 


The  range  Doppler  image  obtained  with  the  default  PRF  of  200  pps  and  128 
pulses  is  that  shown  in  Figure  41b.  The  PRI I  PW  ratio  for  this  case  is  10000,  equivalent 
to  a  duty  cycle  of  110“^.  For  the  case  where  the  PRF  is  201  pps,  the  PRI !  PW  ratio  is 
9950.2  and  the  duty  cycle  is  slightly  increased  to  1.05-10“''. 

To  better  understand  the  phenomena  of  Doppler  ambiguities  when  the  duty  cycle 
is  changed,  different  target  images  are  generated  by  increasing  the  PRF  from  200  pps  and 
keeping  the  pulsewidth  fixed  at  0.5  us.  These  images  are  displayed  in  Figure  45  for  dif¬ 
ferent  cases.  The  quality  of  the  resulting  target  images  is  summarized  in  Table  13.  Of 
special  interest  is  the  image  obtained  in  Figure  45b  with  a  PRF  of  220  pps.  This  plot 
shows  a  target  image  with  a  better  range  resolution  than  the  default  plot  used  throughout 
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this  thesis.  From  the  simulations  results  listed  in  that  table,  it  is  not  clear  when  the  Dop¬ 
pler  ambiguity  is  produced.  However,  this  problem  does  not  occur  if  the  PRI/PW  ratio 
has  an  integer  value.  Such  is  the  case  with  a  coherent  transmitter. 
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Figure  43  DIS  Output  Pulse  Corresponding  to  the  Chirp  Wave  in  Figure  42. 
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Down  Range  Cells 

Figure  44  DIS  Image  for  an  ISAR  with  0.5us  PW  and  20Ipps  PRF. 
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PRF 

PRI  /  PW 

Duty  Cycle 

Doppler  Resolvability 

210 

9523.2 

1.05-1  O'" 

Ambiguous 

220 

9090.9 

l.MO^ 

Good 

230 

8695.7 

1.15-10^ 

Ambiguous 

240 

8333.3 

1.2-10' 

Good 

250 

8000 

1.25-10^ 

Good 

260 

7.6923 

1.3-10" 

Ambiguous 

Table  13  Doppler  Resolvability  for  Different  ISAR  Duty  Cycles 


Figure  45  DIS  Image  for  Different  PRF  Cases  with  a  Fixed  0.5us  PW:  (a)  210  pps,  (b) 
220  pps,  (c)  230  pps,  (d)  240  pps,  (e)  250  pps  and  (f)  260  pps. 
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C.  CHAPTER  SUMMARY 

The  simulation  model  was  used  to  test  the  DIS  performance  and  the  ISAR  image 
generation  under  different  PRFs  and  coherent  integration  times.  Two  main  conclusions 
can  be  extracted  from  this  study.  The  first  is  that  Doppler  resolution  of  the  image  is  in¬ 
versely  proportional  to  the  integration  time,  which  is  equal  to  the  product  of  the  radar 
PRF  and  the  number  of  pulses  integrated.  Since  both  parameters  are  radar  dependent,  the 
Doppler  resolution  of  the  target  image  is  independent  from  the  DIS  hardware.  The  sec¬ 
ond  is  that  an  ISAR  with  a  PRI  that  is  not  an  exact  multiple  of  the  pulse  width  is  subject 
to  phase  incoherence,  which  in  turn,  induces  Doppler  ambiguities  in  the  target  profile.  In 
the  next  chapter,  the  target  image  matrix  is  analyzed  in  terms  of  the  dynamic  range  inten¬ 
sity  in  order  to  relate  it  to  the  DIS  output. 
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VII.  DIS  OUTPUT  AND  IMAGE  INTENSITY  DYNAMIC  RANGE 


The  DIS  output  is  a  dual  channel  I  and  Q  signal,  the  magnitude  of  which  is  de¬ 
pendent  on  the  cumulative  effect  of  the  phase  and  amplitude  modulations  exerted  at  dif¬ 
ferent  range  bins.  The  purpose  of  this  chapter  is  to  study  the  dynamic  range  of  the  output 
signal  amplitude  and  the  effect  on  the  ISAR  image  dynamic  range.  The  study  is  sup¬ 
ported  by  simulation  results  of  the  synthesized  image  for  a  particular  target  using  a  fixed 
number  of  integrated  radar  pulses.  This  target  is  the  same  32-range-bin  target  used  in  the 
study  of  quantization  and  resolution  effects  in  Chapter  V. 

A.  DIS  OUTPUT  DYNAMIC  RANGE 

The  I  and  Q  DIS  output  channels  are  generated  from  the  16-bit  output  adder. 
However,  since  the  image  is  an  array  of  scatterers  arranged  by  range  and  Doppler,  the 
intensity  distribution  cannot  be  related  to  the  output  of  the  16-bit  adder,  where  the 
contribution  of  different  range  bins  have  already  been  superimposed.  Therefore,  the  DIS 
signal  must  be  analyzed  at  the  output  of  each  range  bin. 

The  I  and  Q  range-bin  outputs  are  represented  by  a  13 -bit  format  given  in  Table  5. 
Since  the  format  uses  a  two’s  complement  arrangement,  the  magnitudes  of  the  I  and  Q 
amplitudes  are  represented  by  12  bits.  The  amplitude  dynamic  range,  based  on  the  num¬ 
ber  of  unsigned  bits  x,  is  given  by: 

DR  =  20\og,,(r).  (7.1) 

Thus,  for  the  range-bin  output  x  =  12,  the  dynamic  range  is  72.2  dB.  It  is  expected  that 
the  cumulative  sum  that  takes  place  in  the  follow-on  stage  produces  a  superposition  of 
phasors  that  can  be  successfully  unwrapped  in  the  ISAR  pulse  compression. 

Using  the  32-bin  test  target  with  the  integration  of  128  pulses,  the  DIS  model  is 
run  to  extract  the  output  signal.  The  amplitude  of  two  range-bin  channels  is  captured 
from  the  simulation  model  and  displayed  in  Figure  46.  The  two  selected  signals  were 
respectively  modulated  by  gain  coefficients  of  1024  and  1.  Thus,  they  contain  samples 
with  the  maximum  and  minimum  signal  magnitudes  present  at  any  range-bin  output. 
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Figure  46  Range-bin  Outputs  Signals  for  a  Gain  of  1024  (a.)  and  for  a  Gain  of  1  (b.). 

From  Table  3  of  Chapter  III,  it  was  shown  that  the  non-zero  magnitudes  of  the 
LUT  output  samples  range  from  2’’  to  1  -2’’.  For  the  current  range-bin/pulse  combina¬ 
tion  the  gain  coefficient  is  1024  and  the  binary  samples  are  shifted  10  binary  places.  Fol- 

o  n 

lowing  gain  modulation,  the  samples  are  truncated,  discarding  the  5  LSBs  from  2'  to  2'  . 
The  resulting  range-bin  output  non-zero  magnitudes  vary  from  2“’2“’,  equivalent  to 
eight,  to  a  value  of  (^1-2"’)  2'°,  equivalent  to  1016.  These  magnitude  limits  can  be  seen 

in  the  range-bin  output  plots  shown  in  Figure  46a.  For  the  case  of  a  gain  of  1,  no  shift  of 
binary  places  in  the  range  bin  occurs.  Therefore,  after  a  truncation  of  the  5  LSBs  the 
non-zero  magnitude  varies  from  2'  to  1  -  2'  .  The  extracted  range-bin  output  signal  in 
Figure  46b  is  bounded  by  these  limits. 

The  amplitude  dynamic  range  defined  in  (7.1)  can  also  be  found  graphically.  The 
magnitude  of  the  range-bin  output  signals  can  be  logarithmically  plotted  as  shown  in 
Figure  47.  The  output  dynamic  range  is  then  found  from  the  difference  between  the  peak 
logarithmic  magnitude  for  the  highest  gain  case  (1024)  and  the  lowest  logarithmic  magni¬ 
tude  for  the  lowest  gain  case  (1).  In  terms  of  the  bounding  limits  that  can  be  extracted 
from  the  plots  in  Figure  47,  the  range-bin  output  dynamic  range  is  given  by: 

=  (60.2)-{-12)  =  72.2  dB  (7.2) 

which  matches  the  theoretical  values  obtained  with  Equation  (7.1). 
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Figure  47  Logarithmic  Plot  of  the  Magnitudes  of  the  Range-bin  Signals  in  Figure  46. 


B.  ISAR  PROCESSING  DYNAMIC  RANGE 

The  output  signals  of  different  range  bins  are  superimposed  by  the  chain  of  output 
adders  from  Figure  17  in  Chapter  III.  The  resulting  cumulative  sum  is  then  converted  to 
an  analog  signal  and  transmitted  to  the  ISAR.  The  radar  down  converts  the  signal  and 
performs  range  compression.  To  simplify  this  analysis,  it  is  assumed  that  both  the  at¬ 
mospheric  attenuation  and  the  receiver  processing  do  not  affect  the  DIS  signal  dynamic 
range.  The  last  assumption  is  very  practicable,  since  a  dynamic  range  greater  than  70  dB 
is  typical  in  radar  receivers. 

ISAR  range  compression  in  the  simulation  model  is  implemented  using  cross¬ 
correlation  (Chapter  II,  Figure  6).  The  DIS  signal  is  cross-correlated  with  the  ISAR  ref¬ 
erence  by  multiplying  them  in  the  Fourier  domain,  previously  complex  conjugating  one 
of  them.  Since  both  signals  represent  voltage  values,  the  range  compression  output  signal 
should  represent  a  power  signal.  The  output  of  the  range  compression  is  equivalent  to  a 
one  dimensional  target  profile  on  a  pulse-by-pulse  basis.  The  range  compression  outputs 
for  128  pulses  were  extracted  from  the  simulation  model.  These  128  target  range  profiles 
are  superimposed  in  one  plot  in  Figure  48.  Based  on  the  maximum  and  minimum  mag¬ 
nitudes  from  the  range  profiles,  the  dynamic  range  at  the  range  compression  output  can 
be  obtained  as: 

=  10  l°g..  (4.9  ■  10=  )-10  log,,  (670)  =  28.6  dB.  (7.3) 
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Figure  48  Plot  of  Superimposed  Target  Range  Profiles  for  128  Pulses. 


Figure  49  Logarithmic  Plot  of  the  Target  Range  Profiles. 

The  superimposed  range  profile  is  plotted  logarithmically  in  Figure  49,  where  the  28.6 
dB  of  dynamic  range  can  be  observed  from  the  difference  between  the  maximum  and 
minimum  magnitude.  It  must  be  pointed  out  that  the  intensities  at  a  single  range  profile 
of  the  range  compression  output  are  superimposed  in  terms  of  Doppler.  That  is,  the  com¬ 
plex  value  at  each  location  in  the  one-dimensional  profile  is  the  result  of  the  addition  of 
scatterers  with  different  Doppler  frequencies.  Since  the  Range  Doppler  image  does  not 
superimpose  scatterers  with  different  Doppler  components,  the  dynamic  range  obtained 
in  (7.3)  cannot  be  used  to  derive  the  image  dynamic  range. 
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The  collection  of  range  compression  outputs  form  a  matrix  to  represent  the  time 
history  of  the  target  range  profile.  The  azimuth  compression  produces  a  Fourier  trans¬ 
form  of  this  matrix  along  the  time  dimension.  Thus,  the  scatterers  are  rearranged  accord¬ 
ing  to  their  Doppler  components.  The  result  is  the  target  range-Doppler  image.  Figure 
50  shows  the  range-Doppler  image  along  with  one-dimensional  views  from  the  two  axes. 


Down  Range  Cells 


Figure  50  Range  Doppler  Image:  a)  2D  View,  b)  ID  View  Along  the  Doppler  Axis, 
c)  ID  View  Along  the  Range  Axis. 

From  Figure  50c  the  maximum  intensity  is  8  10'’.  The  minimum  value  extracted 
from  the  simulation  model  is  761.7.  Therefore,  the  dynamic  range  for  the  range  Doppler 
image  can  be  obtained  as: 

=10log,„(8-10‘)-101og,.(761)  =  40.2  dB  (7.4) 

The  ISAR  image  dynamic  range  in  (7.4)  is  higher  than  that  calculated  for  the  DIS 
range-bin  outputs  in  (7.2).  This  can  be  explained  since  the  range-bin  output  samples  are 
not  resolved  in  Doppler  as  the  ISAR  image  pixels.  It  can  be  concluded  that  the  DIS  dy- 
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namic  range  cannot  be  related  to  the  image  dynamic  range,  unless  there  is  a  unique  Dop¬ 
pler  component  per  range  bin.  This  can  be  simulated  using  a  point  target. 

C.  IMAGE  DYNAMIC  RANGE  OF  A  POINT  TARGET 

A  single  scatterer  with  a  unitary  radar  cross-section  and  a  Doppler  component  of 
50  Hz  is  defined  as  the  false  target.  The  model  is  run  integrating  128  pulses.  Both  the 
target  template  and  the  induced  range  Doppler  image  at  the  ISAR  are  shown  in  Figure  51. 
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Figure  5 1  Unitary  RCS  Point  Target:  a)  Target  Template,  b)  Range-Doppler  Image. 

The  generated  point  target  appeared  in  the  ISAR  range  Doppler  profile  at  the 
Doppler  filter,  corresponding  to  50  Hz.  A  smearing  effect  in  the  range  dimension  is  ap¬ 
parent.  This  effect  can  be  observed  in  detail  by  extracting  the  magnitudes  across  the  33'^'* 
Doppler  filter  as  shown  in  Figure  52. 
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Figure  52  Magnitude  Plot  of  the  Doppler  Filter  Containing  the  Point  Target 
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The  smearing  effect  across  the  range  dimension  is  due  to  the  superposition  of 
complex  scatterer  samples  from  different  range  bins  at  the  DIS  output  adders.  Ideally, 
this  range  superposition  is  reversed  at  the  ISAR  range  compression  but  imperfections  in 
the  idealistic  approximation  of  the  inverse  Fourier  transform  shifts  the  range  location  of  a 
few  samples  following  the  Sine  distribution  noted  in  Figure  52.  This  effect  is  identified 
in  Equation  (2.26)  as  the  PSF  of  source  deconvolution.  The  smearing  effect  is  noted  bet¬ 
ter  in  the  three  plots  shown  in  Figure  53.  While  the  range  dimension  in  Figure  53c  shows 
a  wide  spike  at  the  scatterer  location,  the  Doppler  dimension  shows  narrow  spike  in 
Figure  53b. 
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Target  Size  =  1  Range  Cell 
Pulses  integrated  =  128. 

Pulse  Rq)etition  Freq.=  20^ps 
Doppler  Resolution  =  1 .56  Hz 


Figure  53  Point  Target  Range  Doppler  Image,  a)  2D.View.  b)  ID.  View  Along  the 
Doppler  Axis,  c)  ID.View  Along  the  Range  Axis. 


A  logarithmic  plot  of  the  point  target  image  intensities  is  presented  in  Figure  54. 
From  the  maximum  and  minimum  intensity  values,  the  Dynamic  range  for  the  point  tar¬ 
get  image  can  be  obtained  as: 

=64-25  =  39  dB  (7.5) 
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which  is  only  1  dB  smaller  from  the  dynamic  range  image  of  a  ship-like  target  in  (7.4).  It 
can  be  observed  that  the  target  scatter  distribution  barely  affects  the  image  intensity  dy¬ 
namic  range. 


Down  Range  Cells 

Figure  54  One  Dimensional  Logarithmic  Plot  of  the  Point  Target  Image. 

D.  CHAPTER  SUMMARY 

This  chapter  studied  the  amplitude  range  of  the  DIS  samples  in  order  to  relate 
them  to  the  dynamic  ratio  of  the  ISAR  azimuth  compression  output,  which  in  turn  deter¬ 
mines  the  intensity  contrast  in  the  range-Doppler  image.  The  analysis  was  done  for  both 
a  test  ship  target  and  a  point  source  target.  The  results  showed  that  the  dynamic  range  is 
reduced  from  72  dB  in  the  DIS  samples  at  the  output  of  the  DIS  range  bins,  to  40.4  dB  in 
the  Range  Doppler  image.  Similar  results  are  obtained  for  both  the  ship  target  and  the 
point  source  target. 

The  simulation  studies  presented  in  Chapter  V,  VI  and  VII  used  a  32-range-bin 
target  to  reduce  the  long  simulation  times  that  would  involve  larger  targets.  In  the  next 
chapter,  the  current  version  512-range-bin  chip  is  simulated. 
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VIII.  SYNTHESIS  OF  A  512-RANGE  CELL  TARGET 


The  existing  DIS  simulation  models  used  previous  to  this  thesis  do  not  run  prop¬ 
erly  with  target  extents  having  more  than  256  range  eells.  This  ehapter  explains  the  code 
changes  that  are  implemented  to  extend  the  simulation  capabilities  to  any  target  size  and 
describes  a  simulation  run  for  a  512-range  cell  target. 

A.  SOFTWARE  MODEL  MODIFICATION  REQUIREMENTS 

The  original  software  code  of  the  mathost.m  file  created  a  chirp  wave  where  the 
first  sample  is  repeated  a  specific  number  of  times  depending  on  the  number  of  range 
bins  of  the  DIS  array  (target  extent).  The  repetition  of  the  first  phase  originates  from  an 
earlier  DIS  model  in  which  the  DRFM  samples  were  fed  serially,  modeling  the  sample 
delay  between  taps  at  their  inputs.  The  DIS  model  was  eventually  changed  to  the  current 
model  where  the  DRFM  samples  are  fed  in  parallel  to  all  range  bins  and  the  delay  be¬ 
tween  them  is  implemented  in  the  output  adder. 

The  default  values  for  the  pulsewidth  and  sampling  rate  are  0.5  us  and  1  GHz  re¬ 
spectively.  The  resulting  number  of  infra  pulse  samples  is  their  product,  equivalent  to 
500  samples.  For  the  case  of  a  32-range-bin  target,  the  old  model  repeated  the  first  phase 
sample  32  times.  Thus,  the  chirp  that  is  processed  by  the  DIS  routine,  and  later  by  the 
ISAR  compression  routines  had  an  overhead  of  32  samples  (6%).  Nonetheless,  the 
model  worked  correctly  and  synthesized  a  clear  image  of  the  intended  target.  However 
for  a  target  of  512  range  cells,  there  are  more  overhead  samples  than  the  actual  chirp 
phase  samples  (50.6%).  The  initial  solution  to  this  problem  was  to  increase  the  number 
of  infra  pulse  samples  by  increasing  the  sampling  rate.  This  solution  demanded  an  exces¬ 
sive  processing  load  with  long  simulation  time.  The  definite  solution  for  this  problem  is 
to  eliminate  the  overhead  of  repeated  samples  and  apply  code  changes  in  the 
simhwchk_y5.m  file  to  ensure  that  each  range  bin  has  a  zero  output  until  the  time  previous 
range  bin  has  already  generated  a  returned  scattering  sample. 

B.  SIMULATION  RUN 

A  Doppler  template  for  a  512-cell  target  is  defined  in  the  targets  12. mat  data  file. 
Once  loaded,  the  simulation  model  displayed  the  scatterer  arrangement,  as  shown  in 
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Figure  55.  The  scatterer  complex  values  are  processed  using  Equations  (3.2)  through 
(3.6)  to  obtain  magnitude  and  phase  modulation  coefficients. 


Figure  55  Scatterer  Distribution  for  a  5 12-Range-bin  Target  Using  128  Doppler  Cells. 

The  distribution  of  magnitude  coefficients  is  arranged  in  a  range  vs.  pulse  index 
map  in  Figure  56.  This  plot  does  not  show  the  symmetry  condition  observed  for  the  gain 
coefficient  distribution  for  a  32-range  cell  target. 


Down  Range  Cells 

Figure  56  Gain  Coefficient  Distribution  on  a  Range  Cell-Pulse  Index  Map. 
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The  range  vs.  pulse  index  distribution  of  the  phase  modulation  coefficients  is  pre¬ 
sented  in  Figure  57. 
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Figure  57  Phase  Coefficient  Distribution  on  a  Range  Cell-Pulse  Index  Map. 
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The  mathost_v5.m  file  generates  a  chirp-pulse  train,  with  each  pulse  having  the 
waveform  as  displayed  in  Figure  58. 
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Figure  58  Baseband  Plot  of  the  Intercepted  ISAR  Pulse  as  Generated  by  a  DRFM. 
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The  simhwchk_v5.m  file  processes  the  128  pulses  of  the  pulse  train  applying  the 
gain  and  phase  modulation  coefficients  depending  on  the  range  tap  of  the  DIS  array  and 
the  index  of  the  radar  pulse  being  processed.  The  resulting  DIS  output  signal  is  a  pulse 
train  with  modulated  waveforms.  Figure  59  presents  the  DIS  output  waveform  for  the 
first  pulse. 


Plot  of  the  DIS  output  chirp  pulse 


Figure  59  Plot  of  the  DIS  Output  Modulated  Chirp 

The  ISAR  routine  of  the  simhwchk_y5.m  file  processes  each  pulse  returned  to  the 
DIS  with  a  range  compression  process.  For  the  case  of  the  simulation  model,  this  process 
is  implemented  with  a  fast  cross-correlation  as  described  in  Figure  6.  The  FFT  of  the 
first  returned  signal  is  presented  in  Figure  60. 

The  FFT  of  the  chirp  reference  signal  is  presented  in  Figure  61.  After  complex 
conjugation  of  the  latter,  both  FFT  signals  are  multiplied  element  by  element.  The  result¬ 
ing  Fourier  domain  signal  is  presented  in  Figure  62. 
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Figure  60  FFT  of  the  DIS  Output  Signal. 


ISAR  Reference  Intrapulse  Samples  (Fourier  Domain) 
Figure  61  FFT  of  the  Chirp  Reference  Signal. 
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Figure  62  FFT  Product  Signal  at  the  Pulse  Fast  Correlation  Output. 


Figure  63  Fast  Correlation  Output  Showing  the  Target  Range  Profile. 
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Figure  64  Surface  Plot  of  Range  Profile  Matrix. 
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Figure  65  Range  Doppler  Image  for  a  False  Target  of  5 12  Range  Cells. 
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C.  CHAPTER  SUMMARY 

This  chapter  presented  the  simulation  model  for  the  current  DIS  chip  design  with 
512  range-bin  processors.  From  the  values  selected  for  DRFM  sampling  frequency  (1 
GHz)  and  radar  PW  (5us),  the  resulting  number  of  samples  at  each  pulse  is  500.  Previous 
DIS  software  models  were  unable  to  simulate  a  scenario  where  the  target  extent  had  as 
many  range  bins  as  the  number  of  samples  in  a  radar  pulse.  The  software  improvements 
applied  to  this  model  made  possible  a  simulation  of  a  target  whose  extent  in  range  cells 
(512)  was  superior  to  the  number  of  radar  pulse  samples  (500).  It  is  expected  that  the 
model  would  perform  similarly  for  larger  targets. 

The  new  version  of  the  DIS  software  model  is  able  to  represent  the  current  DIS 
hardware  design.  The  simulation  results  at  this  level  showed  that  the  hardware  design 
works  properly.  The  results  must  be  supported  by  lower  level  DIS  hardware  simulations 
using  VHDL. 
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IX.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  investigated  the  theoretical  operation  of  the  DIS  integrated  circuit  as 
an  ISAR  countermeasure  and  developed  a  software  simulation  model  to  test  the  perform¬ 
ance.  The  software  model  simulates  the  main  signal  processing  steps  that  take  place  in 
the  DIS  and  later  in  the  ISAR  receiver  and  processor.  The  extraction  and  plotting  of  sig¬ 
nals  processed  by  the  simulation  model  provided  a  good  tool  to  study  the  DIS  design  in 
depth. 

A.  CONCLUSIONS 

•  The  DIS  software  simulation  model  provided  with  this  thesis  is  better  than 
the  predecessors  in  being  closer  in  the  resemblance  to  the  actual  hardware, 
faster  to  execute  and  shorter  in  code  extension.  Thus,  the  DIS  design  team 
has  a  more  efficient  and  easier  to  understand  simulation  tool. 

•  A  method  to  evaluate  the  DIS  performance  based  on  image  comparison 
from  an  ideal  model  was  integrated.  The  performance  under  the  different 
simulation  scenarios  could  be  evaluated  based  on  the  quality  of  the  output 
image  and  the  similarity  to  the  ideal  model  image.  The  image  comparison 
tool  made  it  possible  to  isolate  individual  processing  stages  in  order  to  find 
the  major  contributors  to  the  overall  image  error.  This  tool  can  also  be 
used  to  measure  the  impact  of  future  design  changes  at  any  stage  of  the 
DIS  process. 

•  The  output  adder  overflow  was  identified  as  the  main  source  of  image  er¬ 
ror.  This  thesis  suggests  that  the  overflow  can  be  managed  by  solely  con¬ 
trolling  the  thresholding  values  that  assign  the  highest  DIS  gain  coefficient 
value  (1024),  thus  preventing  any  hardware  design  change. 

•  Software  improvements  in  the  simulation  code  enabled  the  modeling  of  a 
target  with  an  extension  of  512  range  cells.  These  code  changes  elimi¬ 
nated  a  theoretical  limit  for  the  simulation  even  bigger  targets. 
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•  The  simulation  model  was  efficiently  used  to  produce  input  and  output 
vectors  for  VHDL  simulations  for  DIS  arrays  of  different  sizes.  Output 
results  of  both  the  MATLAB  simulation  model  and  the  VHDL  model 
matched  perfectly. 

•  The  software  model  has  been  described  thoroughly  in  order  to  make  it  eas¬ 
ier  to  understand  when  implementing  future  code  changes. 

B.  RECOMMENDATIONS 

•  The  software  code  can  still  be  improved  by  substituting  for  some  of  the  it¬ 
erative  loop  structures  with  vector  and  matrix  operations.  These  changes 
would  increase  the  percentage  of  parallel  processing  operations,  thus 
greatly  improving  the  simulation  time.  Additional  time  saving  changes 
include  the  reduction  of  intermediate  variables,  the  use  of  column  opera¬ 
tions  (vice  rows)  and  the  reduction  of  conditional  (if-then)  operations. 
Care  must  be  taken  to  keep  the  software  code  easy  to  understand  if  exces¬ 
sive  reduction  in  code  lines  is  achieved. 

•  The  performance  of  the  DIS  against  stepped  chirp  pulses  can  be  added  to 
this  study.  From  the  quantization  analysis  obtained  from  chirp-pulse 
cases,  it  can  be  predicted  that  the  DIS  should  be  able  to  deceive  stepped 
chirp  pulses  that  resemble  a  discrete  linear  FM  chirp  and  scrambled  fre¬ 
quency  stepping  pulses  (discrete  chirp  with  scrambling  of  pulse  samples). 
In  the  particular  case  of  interpulse  frequency  stepping  where  each  pulse 
has  a  different  constant  frequency,  the  DIS  should  not  be  effective. 
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APPENDIX  A.  MATLAB  CODES 


1.  DIS  SIMULATION  FILES  -  VERSION  5 
a.  extract_para_ship_v5.m 

(^************************************************************************************************************ 

%  extract_para_ship_v5.nn  * 

% 

%  -  Extracts  modulation  coefficients  (phase  and  magnitude)  from  a  false  * 

%  target  to  be  used  by  DIS  simulation  files.  * 

%  -  Runs  'doppler_XX_target.m'  to  create  a  Range-Doppler  target  profile.  * 

%  -  Saves  modulation  coefficients  in  'paraMULTIq4Ship.txt'.  * 

% 

%  Modified  by  F.LeDantec  7/7/02  based  on  previous  version  from  P.E.Pace  * 

% 

(^************************************************************************************************************ 


clear  all; 

global  dp_pts  rg_pts  rag_pts  numtaps  tO  k  prf  nbitsdop 


numtaps  =  32 

%  #  of  taps  (target  extent) 

dp_pts  =  128 

%  #  of  Doppler  bins  (number  of  pulses  integrated) 

bw2  =  500e6 

%  Bandwidth  of  the  chirp  signal,  delta 

fs  =2*  bw2 

%  Sampling  frequency 

prf  =  200; 

%  Pulse  repetition  frequency 

nbitsdop  =  5; 

%  #  of  bits  used  for  phase  quantization 

pwc  =  1/(bw2); 

%  Compressed  pulsewidth 

pw  =  0.5e-6*200/prf: 

%  Uncompressed  pulsewidth 

pri  =  l/prf; 

%  Pulse  repetition  interval 

k  =  bw2/pw: 

%  Pulse  compression  rate  (deltaBW/uncompr.pWidth) 

Ts  =  1/fs: 

%  Sampling  time 

to  =  Ts:Ts:pw: 

%  Interpulse  sample  index 

rag_pts  =  length(tO); 

%  Number  of  samples  per  chirp  pulse  (range  cells) 

rg_pts  =  max(rag_pts,numtaps*10);  %  #  of  range  FFT  pts  including  padding 

p  =  2*pi/(2''nbitsdop): 

%  Quantization  factor 

tstop  =  pri*dp_pts; 

%  Pulse  train  duration  (integration  time) 

ti  =  linspace(0, tstop, dp_pts): 

%  Vector  with  pulse  leading  edge  times 

switch  numtaps 

case  2, 

doppler_02: 

%  Creates  and  saves  a  2  range-bin  profile 

load  target02 

%  Loads  profile 

case  32, 

doppler_32: 

%  Creates  and  saves  a  32  range-bin  profile 

load  target32 

%  Loads  profile 

case  128, 

doppler_1 28; 

%  Creates  and  saves  a  32  range-bin  profile 

load  target128 

%  Loads  profile 

case  512, 

doppler_512: 

%  Creates  and  saves  a  512  range-bin  profile 

load  target512 

%  Loads  profile 

end 
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%  Doppler  profile  is  reduced  (slower  rotation) 


if  dp_pts  <  1024, 
freq  =  freq/IO; 
end; 

%  PLOT  OF  THE  FALSE  TARGET  RANGE-DOPPLER  PROFILE  : 
figure; 

colordef  biack 
plot(freq,'y*'): 
nt  =  nunn2str(nunntaps): 
nd  =  nunn2str(dp_pts): 

title([  nt  '-range  bin  Faise  Target  Profile  - '  nd  '  Radar  Pulses'  ]) 

xlabel('Range  bins'); 

yiabei('Doppier  frequency  [Hz]'); 

axis([0  fix(numtaps/8*9)  0  128]); 

axis  square 


%  CREATION  OF  A  COMBINED  SIGNAL  AT  EACH  RANGE  GATE,  BY  SUPERPOSITION  OF 
%  COMPLEX  RETURNS,  TO  OBTAIN  A  1  DIM  RANGE  PROFILE  FOR  EACH  PULSE  : 


%  Scatterer  Amplitudes  (uniform) 


%  At  each  range  bin 


A  =  ones(numtaps,dp_pts); 
targetSum  =  zeros(numtaps,dp_pts); 
for  i1  =  1:numtaps, 

target  =  zeros(dp_pts,length(ti)); 

for  i2  =  1  :numfreq(i1 ),  %  At  each  scatterer 

for  i3  =  1  :length(ti),  %  At  each  integrated  pulse 

target(i2,i3)  =  A(i1,i2)  *  exp(-j*2*pi*freq(i1,i2)*ti(1,i3)); 
end  %  Creates  complex  scatterer 

end 

for  i4  =  1:numfreq(i1)  %  Combined  signal  p/rbin  p/pulse 

targetSum(i1,:)  =  targetSum(i1,:)  +  target(i4,:); 
end 
end 


%  CREATION  OF  GAIN  MODULATION  COEFFICIENTS  : 


amp  =  abs(targetSum); 
gain  =  ones(numtaps,dp_pts); 
fq  =  1.5; 

for  i1  =  1:numtaps, 
for  i2  =  1  :dp_pts, 

normAmp  =  amp(i1,i2)/max(max(amp)); 
for  s=1:10, 

if  normAmp  >  (2'^-(1 1-s))*fq, 
gain(i1,i2)  =  2^s; 
end 
end 
end 
end 


%  Magnitude  of  combined  signal. 
%  Default  gain  coefficient  is  1 . 

%  Arbitrary  threshold  factor 


%  Magnitudes  are  normalized. 

%  Thresh:  0.001 5....,0.375,0.75 

%Gain:  2,4,8,16, . 512,1024 

%  Assigns  1 1  power-of-two  gains 


%  CREATION  OF  PHASE  MODULATION  COEFFICIENTS  : 


%  Phase  of  combined  signal 
%  Adjusts  values  from  0  to  31 .99 
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vphasel  =  angle(targetSum); 
vphase2  =  vphasel  *2''nbitsdop/(2*pi); 


phaseinc  =  zeros(nunntaps,dp_pts); 

for  i6  =  1  inumtaps,  %  Turn  phases  into  phase  increments 

for  i7  =  1  :dp_pts, 
if  i7==1, 

phaseinc(i6,i7)  =  phaseinc(i6,i7)+vphase2(i6,i7): 
else  %  Finds  phase  increments  bet/pulses 

phaseinc(i6,i7)  =  phaseinc(i6,i7-1  )+vphase2(i6,i7-1  )-vphase2(i6,i7): 
end 
end 
end 

for  i8  =  1  :dp_pts, 

for  i9  =  1  :numtaps,  %  5-bit  quantization  in  32  levels. 

phasecoeff(i9,i8)  =  fix(rem(phaseinc(i9,i8)+32,32)): 
end 
end 


%  MODULATION  COEFFICIENTS  ARE  STORED  SERIALLY,  SO  THEY  CAN  BE  LOADED  BY 
%  THE  DIS  SIMULATION  PROGRAM  : 


%  Stores  modulation  coefficients 
%  Store  #  of  chirp  pulse  samples 
%  Store  #  of  Doppler  bins 
%  Store  #  of  taplines  (target  extent) 


%  Stores  gain  coefficients 


f4  =  fopen('paraMULTIq4Ship.txt',W): 
fprintf(f4,'%d\r\n',rag_pts): 
fprintf(f4,'%d\r\n',dp_pts): 
fprintf(f4,'%d\r\n',numtaps): 
for  aa=1:dp_pts, 
for  bb=1:numtaps, 
fprintf(f4,'%d\r\n',gain(bb,aa)): 
end 
end 

for  aa  =  1  :dp_pts, 
for  bb  =  1  :numtaps, 

fprintf(f4,'%d\r\n',phasecoeff(bb,aa));  %  Stores  phase  coefficients 
end 
end 

fclose(f4): 


%  PLOT  OF  THE  GAIN  QUANTIZATION  SCHEME 
figure 

plot(amp/max(max(amp)),gain,'*'):  %  Plots  gain  coefficient  distrib. 


b.  mathost_v5.m 

0/************************************************************************************************************ 

/O 

%  mathost_v5.m  * 

% 

%  -  Creates  and  stores  an  intercepted  chirp  pulse  as  generated  by  a  * 

%  phase  sampling  DRFM.  * 

%  -  Creates  the  reference  signal  required  for  ISAR  range  compression.  * 

%  -  Runs  'simhwchk_v5.m'  to  execute  the  DIS  and  ISAR  simulation.  * 

%  -  Runs  'image_analysis.m'  to  compare  image  w/infinite  resolution  case.  * 

% 

%  Modified  by  F.LeDantec  7/7/02  based  on  previous  version  from  P.E.Pace  * 

% 

(^  ************************************************************************************************************ 


global  dp_pts  rg_pts  numtaps  rag_pts  k  tO  prf  nbitsdop 
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%  GENERATION  OF  DRFM  SIGNAL  CONTAINING  PHASE  SAMPLES  FOR  ALL  PULSES  : 
pri_rg_phaseq  =  zeros(dp_pts,rag_pts): 
p  =  2*pi/(2''nbitsdop): 
pri  =  1/prf: 
for  idx1  =  1  :dp_pts, 
t1  =  to  +  idx1*pri: 
oldphase  =  2*pi*(k*t1  .*t1  )/2: 
oldphase  =  rem(oldphase,2*pi): 
prLi'g_phaseq(idx1,:)  =  fix(oldphase/p): 
crefc(:,idx1)  =  exp(j*oldphase.'): 
end 

%  STORAGE  OF  DRFM  PHASE  MATRIX  : 
f2  =  fopen('rawlnt.txt','w'): 
for  i  =  1  :dp_pts, 

int_raw  =  pri_rg_phaseq(i,1:rag_pts-1): 
fprintf(f2,'%d,',int_raw): 
int_raw  =  pri_rg_phaseq(i,rag_pts): 
fprintf(f2,'%d\r\n',int_raw); 
end; 

fclose(f2): 

%  GENERATION  OF  ISAR  REFERENCE  FOR  RANGE  COMPRESSION  (FAST  CORRELAT.) : 

cref  =  conj(fft(crefc,2*rg_pts-1 ));  %  FFT  of  ISAR  reference 

save  pc_ref  cref  t1  %  Saves  FFT  of  ISAR  reference. 

%  DIS  SIMULATION,  IMAGE  GENERATION  AN  ANALYSIS  : 

sinnhwchk_v5:  %  DIS  simulation  &  ISAR  procesing 

if  numtaps==32  &  nDopplerCell==128  &  prf==200, 

lmage_Analysis:  %  Comparison  with  ideal  image 

end; 


%  Initializes  DRFM  phase  matrix 
%  Phase  quantization  factor 

%  Loop  per  pulse 
%  PRI  time  shift  for  pulse  samples 
%  Chirp  pulse  phase  samples 
%  limit  phase  from  0  to  2pi 
%  Phase  quantized  in  2''nbit  levels 
%  Store  phase  for  ISAR  reference 


%  Opens  file  for  DRFM  ph.  storage 
%  Loop  for  all  pulses. 

%  Stores  pulse  phase  samples 

%  Stores  last  pulse  phase  sample 


c.  simhwchk_v5.m 

0/************************************************************************************************************ 

/O 

%  simhwchk_v5.m  * 

% 

%  Simulates  the  DIS  architecture  by:  * 

%  -  Loading  its  inputs  from  files  'paraMULTIq4Ship.txt'  and  'rawint.txt'  * 

%  -  Processing  the  modulations  of  the  tap  array.  * 

%  -  Generating  the  output  chirp.  * 

%  Simulates  the  ISAR  processing  and  plots  the  range-Doppler  profile  of  * 

%  the  false  target  as  seen  in  the  radar  screen.  * 

% 

%  Modified  by  F.LeDantec  7/7/02  based  on  previous  version  from  P.E.Pace  * 

% 

0/************************************************************************************************************ 

/O 


clear 

global  rg_pts  prf 

fid  =  fopen('paraMULTIq4Ship.txt','r'):  %  Opens  file  w/DIS  coefficients 

tmp  =  fscanf(fid,'%f):  %  Loads  file  contents 

nRangeCell  =  tmp(1 );  %  number  of  pulse  samples 
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nDopplerCell  =  tmp(2): 
numtaps  =  tnnp(3): 

gainin  =  tmp(4:4+numtaps*nDopplerCell-1): 
gain  =  reshape(gainln, numtaps, nDopplerCell); 
gain  =  gain'; 

phi  =  tmp(4+numtaps*nDopplerCell:end); 
phasinc  =  reshape(phi, numtaps, nDopplerCell); 
fclose(fid); 

raw  =  zeros(nDopplerCell,nRangeCell); 
fid  =  fopen('rawint.txt','r'); 
for  r  =  1:nDopplerCell, 
for  i2  =  1  :nRangeCell-1 , 
raw(r,i2)=fscanf(fid,'%f  ,1 ); 
comma  =  fscanf(fid,'%c',1); 
end 

raw(r,nRangeCell)  =  fscanf(fid,'%f,1); 
end 

fclose(fid); 

[row, col]  =  size(raw); 

phaDRFM  =  [raw,zeros(row,numtaps-1)]; 


%  number  of  pulses  integrated 
%  number  of  DIS  range  taps 
%  Extracts  gain  coefficients 
%  Arranges  them  in  a  matrix 
%  Transpose 

%  Extracts  gain  coefficients 
%  Arranges  them  in  a  matrix 

%  Initializes  DRFM  phase  matrix 
%  Opens  file  w/DRFM  samples 
%  Loop  per  pulse 
%  Loop  per  sample 
%  Loads  DRFM  phase  sample 
%  Loads  comma  character 

%  Last  pulse  phase  sample 


%  Trailing  zeros  to  clear  taps 


colordef  white 

set(0,'defaultAxesFontSize',12);  %  Font  size  for  graph  labels 

figure  %  Plots  one  intercepted  pulse 

plot(  real(  exp(j*raw(  1 ,1  mRangeCell)  /32*2*pi ) ) ) 
title('Plot  of  an  intercepted  ISAR  chirp  as  generated  by  a  DRFM'); 


depthLUT  =  32; 
load  -ascii  cosine8.txt 
load  -ascii  sine8.txt 
lutOut  =  zeros(numtaps,1); 

tapOut  =  zeros(nRangeCell  +  (numtaps-1),numtaps); 

o_res  =0.25; 

overflow=  0; 

peak  =  0; 

limit  =8192; 

partial_tapsum  =  zeros(nRangeCell+numtaps-1,1); 


%  Size  of  the  LUT  table 
%  Loads  LUT  tables 

%  Initializes  LUT  output 

%  Minimum  output  bit  resolution 
%  Initializes  overflow  counter 
%  Initializes  maximum  amplitude 
%  Adder  overflow  limit  2'^13 


%  Loop  per  ISAR  pulse 


%  Per  out  pise  sample 
%  Loads  DRFM  Pha  in// 


for  batchCnt  =  ImDopplerCell, 

disp(['Processing  Pulse  ',num2str(batchCnt)]); 
for  inPIsCnt  =  1;(nRangeCell+numtaps-1), 

tapln(1  :numtaps)  =  phaDRFM(batchCnt,inPlsCnt); 
phAddOut=  tapln(1:numtaps)'+phaslnc(;, batchCnt);  %  Pha  adder  output 
tmp  =  rem(phAddOut+depthLUT,depthLUT)+1 ;  %  Phas  limit  bet  0-31 

lutOut  =  cosine8(tmp)  +  j*sine8(tmp);  %  LUT  table  output 

gainOut  =  gain(batchCnt,:)'.*lutOut;  %  Gain  modulation 

if  inPIsCnt  >  nRangeCell,  %  Shuts  down  taps 

gain0ut(1  :(inPlsCnt-nRangeCell))  =  0;  %  after  pulse  passes 

end 

for  idx3=1:numtaps, 

delay  =  idx3-1 ;  %  Out.adder  register  delay  bet/taps 

tapOut(inPlsCnt+delay,idx3)=fix(gainOut(idx3)/o_res)*o_res; 
end  %  Applies  minimum  register  resolution 

sum1=0; 

for  tt=min(inPlsCnt,numtaps):-1  :max(1  ,inPlsCnt-nRangeCell), 
sumi  =  sumi  +  tapOut(inPlsCnt,tt); 
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lout=real(sum1): 

Qout=imag(sum1): 
if  lout  >=  limit, 

sum1  =  iout-2*iimit  +  j*Qout: 
overflow  =  overflow  +  1 ; 
elseif  lout  <  -limit, 

sum1  =  lout+2*limit  +  j*Qout: 
overflow  =  overflow  +  1 ; 
end 

if  Qout  >=  limit, 
sum1  =  lout  +  j*(Qout-2*limit): 
overflow  =  overflow  +  1 ; 
elseif  Qout  <  -limit, 

sum1  =  lout  +  j*(Qout+2*limit): 
overflow  =  overflow  +  1 ; 
end 

peak=max(max(abs(lout),abs(Qout)),peak): 
partial_tapsum(inPlsCnt,1 )  =  sum1 ; 
end 
end 

finalAdderOut(batchCnt,:)=partial_tapsum.': 

end 

overflow 

peak 

dac_res  =  0.25; 

finalAdderOut  =  fix(finalAdderOut/dac_res)*dac_res: 
figure 

plot(  real(  finalAdderOut(5,:)  )),axis  tight 
title('Plot  of  the  DIS  output  chirp  pulse'); 


%  Overflow  simulation  using  +/-limit 

%  Positive  overflow  (in-phase) 

%  Adder  output  w/overflow 
%  Increases  counter 
%  Negative  overflow  (quad-pha) 

%  Adder  output  w/overflow 
%  Increases  counter 

%  Positive  overflow  (in-phase) 

%  Adder  output  w/overflow 
%  Increases  counter 
%  Negative  overflow  (quad-pha) 

%  Adder  output  w/overflow 
%  Increases  counter 

%  End  of  overflow  simulation. 


%  DIS  output  mod.  chirps 

%  Number  of  DIS  overflows 
%  Maximum  DIS  amplitude 
%  DAC  resol  0.25  for  16bit 


%  Plots  1  DIS  output  chirp  pulse 


%  ISAR  PROCESSING  : 


load  pc_ref  % 

priRgMapShift  =  zeros(fix(numtaps*9/8),nDopplerCell): 

%  Range  compression: 

pcRefMapShift  =  fft(finalAdderOut.',2*rg_pts-1 ).';  % 

for  idx  =  1  :nDopplerCell,  % 

tmp  =  cref(:,idx).*pcRefMapShift(idx,:).':  % 

tmp2  =  ifft(tmp):  % 

priRgMapShift(:,idx)  =  tmp2(  1  :fix(numtaps*9/8) );  % 

end 

%  Azimuth  compression: 

%  FFT  of  matrix  of  range  profiles  in  the  pulse  index  dimension  (time) 

%  gives  the  target  range  Doppler  profile. 
dpRgMapShiftMOD  =  abs(fft(priRgMapShift.',nDopplerCell)); 
quantimage  =  dpRgMapShiftMOD(1  :nDopplerCell*100/prf,1  :fix(numtaps/8*9)): 
save  finresoldata  quantimage  %  Saves  Range  Doppler  profile 

%  ISAR  IMAGE  PLOTS  (Range  Doppler  profile): 


Loads  FFT  of  ISAR  reference 


FFT  of  DIS  output 
For  all  pulses. 

Fast  correlat.(FFT  product) 
InvFFT  to  obtain  range  profile 
Combined  range  profile  per  pulse 


dopp  =  (1:size(quantimage,1))*prf/nDopplerCell: 
rang  =  1:size(quantimage,2); 
colordef  black 
figure; 

Ncontours  =  120; 

contour(rang,dopp,quantimage,Ncontours); 


%  Doppler  axis 
%  Range  axis 
%  Background  color 

%  #  of  plot  contour  lines 
%  Contour  Plot 
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%title('Range-Doppler  profile  of  a  false  target'); 
xlabel('Down  Range  Cells'): 
ylabel('Doppler  frequency  [Hertz]'); 
axis  tight 
figure: 

colormap(hot(100)) 

surf(rang,dopp,quantinnage),  shading  interp,  lighting  phong,  view(0,90): 

%title('Range-Doppler  profile  of  a  false  target'); 

xlabel('Down  Range  Cells'): 

ylabel('Doppler  frequency  [Hertz]'); 

axis  tight 


d.  image_analysis.m 


(^***************************************************************************************************** 

%  innage_analysis.nn  * 

% 

%  -  Load  the  matrix  representing  the  real  DIS  range-Doppler  image.  * 

%  -  Load  the  matrix  representing  the  ideal  DIS  image  range-Doppler.  * 

%  -  Calculates  the  Mean  Squared  Error  coefficient  (average  noise  energy)  * 

%  -  Calculates  the  Two-dimensional  correlation  coefficient.  * 

%  -  Plots  the  noise  energy  image  distribution.  * 

% 

%  Fernando  LeDantec  7/7/02.  * 


% 

(^***************************************************************************************************(^ 


dp_pts  =  128; 
numtaps  =  32; 
load  finresoldata 
load  infresoldata 


%  Number  of  Doppler  cells 
%  Number  of  range  cells 
%  Loads  real  DIS  image 
%  Loads  ideal  DIS  image 


normimage  =  image/max(max(image)):  %  Normalizes  real  image, 

normquantimage  =  quantimage/max(max(quantimage)):  %  Normalizes  ideal  image. 


for  m=1:dp_pts/2, 

for  n=1  :numtaps*1 .0625, 

noiseEnergy(m,n)  =  (normimage(m,n)-normquantimage(m,n))''2;  %  Pixel  square  error, 
end; 
end; 


[a,b]  =  size(noiseEnergy): 

noise_ave  =  sum(sum(noiseEnergy))/(a*b)  %  MSE  coefficient. 

corel_coeff  =  corr2(image,quantimage)  %  2D  correlation  coeff. 


figure:  % 

dopp  =  (1  :a)*prf/nDopplerCell;  % 

rang  =  1  :b;  % 

colormap(hot(100)) 

surf(rang, dopp, noiseEnergy), lighting  phong,  view(0,90): 
title('Range-Doppler  Map  of  the  Quantization  and  Overflow 
xlabel('Range  Cells'); 
ylabel('Doppler  Cells'): 
axis  tight; 


Plots  Noise  energy. 
Doppler  axis 
Range  axis 


Noise'): 
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APPENDIX  B.  DIGITAL  RF  MEMORY 


1.  DRFM  THEORY 

Most  military  radars  currently  operate  using  coherent  processing,  where  the  signal 
phase  is  preserved  for  detection  and  performance  and  low  vulnerability  to  repeater  jam¬ 
ming.  Among  EA  target-like  jamming  techniques  are  those  that  radiate  a  signal  similar 
to  the  signal  scattered  by  the  target.  Selected  signal  parameters  are  then  changed  to  de¬ 
ceive  or  confuse  the  radar  operator.  The  advent  of  coherent  radar  processing  proved  to 
be  very  powerful  to  sustain  non-coherent  deception  and  confusion  jamming  techniques 
with  older  repeaters.  Jamming  technology  was  forced  to  use  coherent  techniques.  The 
development  of  the  Digital  Radio  Frequency  Memory  (DRFM)  provided  a  technological 
breakthrough  to  achieve  this  effort. 

The  DRFM  concept  was  developed  to  store  a  signal  and  replicate  it  coherently. 
Thus,  a  DRFM  can  be  used  for  coherent  jamming  techniques.  Figure  66  shows  a  block 
diagram  representation  of  a  DRFM.  The  first  stage  is  the  frequency  down-conversion  of 
the  signal  into  an  intermediate  frequency  (IF).  This  is  implemented  using  two  very  stable 
local  oscillators  (LO)  in  quadrature  that  generate  two  channels  I  and  Q.  High  stability  of 
these  oscillators  is  essential  in  order  to  preserve  the  phase  information  of  the  signal.  Both 
channels  are  fed  into  ADCs  for  sampling  at  intervals  of  a  few  microseconds.  The  sam¬ 
ples  are  quantized  into  digital  formats  and  stored  in  digital  Random  Access  Memory 
(RAM).  Once  the  signal  is  in  memory,  it  is  regenerated  after  any  user-defined  delay. 
The  signal  I  and  Q  samples  are  loaded  from  RAM,  fed  to  DACs  and  up-converted  in  fre¬ 
quency  using  the  same  LO  signal. 

The  signal  samples  stored  in  RAM  can  be  managed  for  a  wide  variety  of  ECM 
application.  The  DIS  described  in  this  thesis  certainly  exploits  this  capacity.  Other  op¬ 
tional  processing  techniques  include  the  generation  of  noise,  synthesis  of  a  different  in- 
band  frequency,  and  modulation  of  false  targets. 

The  conversion  of  microwave  signals  into  stream  of  bits  is  limited  by  the  speed  of 
ADCs  and  semiconductor  memories.  Hence,  the  use  of  an  intermediate  frequency  is  es¬ 
sential  to  match  the  converter  sampling  and  RAM  access  speed.  For  example,  a  10-GHz 
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signal  can  be  down-converted  to  a  1  GHz  intermediate  frequency.  According  to  the  Ny- 
quist  criterion,  the  sampling  frequency  should  be  twice  the  signal  frequency.  Therefore, 
the  IF  signal  should  be  sampled  at  2  GHz.  However,  since  the  sampling  is  done  simulta¬ 
neously  in  the  I  and  Q  channels,  each  of  them  can  sample  at  half  of  the  required  sampling 
frequency  of  1  GHz.  This  sampling  frequency  matches  current  RAM  access  speeds. 


Figure  66  Functional  Diagram  for  a  DRFM. 
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APPENDIX  C.  THE  AN/APS-137  RADAR 


The  APS- 137  is  a  maritime  surveillanee  radar  used  by  the  U.S.  Navy  onboard  the 
S-3B  anti-submarine  warfare  aircraft  and  will  soon  be  installed  onboard  the  P-3C  mari¬ 
time  patrol  aircraft.  The  system,  manufactured  by  Raytheon,  provides  long  range  detec¬ 
tion  and  ship  classification  capability  using  ISAR  imaging  among  other  operational 
modes.  The  radar  operating  modes  include  surface  surveillance  (land,  sea  and  littoral), 
navigation,  periscope  detection,  strip-map  SAR,  spotlight  SAR  and  ISAR.  The  processed 
signals  are  displayed  in  a  variety  of  options:  PPI,  B-scan,  A-scan  range  profding,  SAR 
map  and  ISAR  imagery. 

Periscope  detection  is  achieved  using  a  high  PRF  and  a  high  anterma  scan  rate. 
The  high  resolution  obtained  provides  periscope  detections  to  a  range  up  to  26  nm  while 
flying  at  1500  feet.  The  ISAR  processor  differentiates  this  radar  from  its  predecessor,  the 
APS- 116,  used  on  S-3A  aircraft.  The  ISAR  processor,  coupled  with  onboard  classifica¬ 
tion  aids,  is  a  powerful  tool  to  identify  rolling  or  pitching  surface  targets  before  an  anti¬ 
surface  strike.  ISAR  requires  a  tight  beam  be  held  on  the  target  ship  to  measure  the  Dop¬ 
pler  shifts  as  the  target  rotates. 

The  APS- 137  radar  uses  two  antennas,  a  narrow  beam  parabolic  antenna  for  ISAR 
applications  in  a  search  light  mode,  and  a  flat  plate  array  for  sector  scan  and  normal  sur¬ 
veillance. 


Figure  67  The  S-3B  Viking. 
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APS-137  Parameters  [21]: 


Frequency: 

9.3-10.1  GHz. 

Pulsewidth: 

5  us  (periscope  detection) 

12  us  (SAR  modes) 

10  us  (all  other  modes) 

Pulsewidth  compression: 

2.5ns  (periscope  detection) 

6  to  13  ns  (ISAR  mode) 

14  ns  (search  modes) 

240  ns  (navigation  mode) 

PRF: 

388  pps  (search/navigation) 

500  to  1000  (SAR-ISAR) 

1854  pps  (periscope  detection) 

Best  resolution: 

0.9  m  (ISAR) 

3  m  (SAR) 

Azimuth  coverage: 

360° 

Average  power: 

200  W  (search/navigation). 

500  W  (ISAR). 

350  W  (SAR). 

460  W  (periscope  detection). 
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APS-137  Antennas: 


Antenna  rotation: 
Parabolic  antenna: 

Planar  array  antenna: 


6,60  and  300  rpm. 

1 .3°  H  beamwidth. 
4.5°  V  beamwidth. 
35  dB  Gain 
1.05°  H  beamwidth, 
15°  V  beamwidth. 
32.5  dB  Gain. 
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